{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Varying initial concentrations and parameters" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "A version of this notebook may be run online via Google Colab at https://tinyurl.com/rxd-initial-and-parameters\n", " (make a copy or open in playground mode)." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:04.934614Z", "iopub.status.busy": "2025-08-18T03:36:04.934463Z", "iopub.status.idle": "2025-08-18T03:36:05.616406Z", "shell.execute_reply": "2025-08-18T03:36:05.615944Z" } }, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Often we will want to see how the choice of initial conditions affects the dynamics. We can do this by setting the initial attribute of an `rxd.Species` and rerunning." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For example, suppose at a single point we have the bistable dynamics introduced in the first part of this tutorial. That is, $u'=-u(1-u)(\\alpha - u)$. (Here we use $u$ instead of a specific molecule name to indicate that we are not describing any particular molecule's kinetics, but rather modeling a class of phenomena.)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This time, we'll use an `rxd.Parameter` for $\\alpha$ instead of a constant. This offers two advantages: (1) this allows $\\alpha$ to vary spatially, and (2) this allows us to change the parameter values and rerun without changing the reaction/rate specification." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We start by defining the model and setting up the recordings:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:05.619817Z", "iopub.status.busy": "2025-08-18T03:36:05.619587Z", "iopub.status.idle": "2025-08-18T03:36:06.487859Z", "shell.execute_reply": "2025-08-18T03:36:06.487426Z" } }, "outputs": [], "source": [ "from neuron import n, rxd\n", "from neuron.units import mV, ms, mM\n", "from matplotlib import pyplot as plt\n", "\n", "n.load_file(\"stdrun.hoc\")\n", "\n", "soma = n.Section(name=\"soma\")\n", "cyt = rxd.Region([soma], name=\"cyt\", nrn_region=\"i\")\n", "u = rxd.Species(cyt, name=\"u\")\n", "α = rxd.Parameter(cyt, value=0.3 * mM)\n", "\n", "rate = rxd.Rate(u, -u * (1 * mM - u) * (α - u))\n", "\n", "n.finitialize(-65 * mV)\n", "\n", "t = n.Vector().record(n._ref_t)\n", "y = n.Vector().record(soma(0.5)._ref_ui)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now we define a function that sets the initial value for ca, runs the simulation, and plots the curves:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:06.490769Z", "iopub.status.busy": "2025-08-18T03:36:06.490126Z", "iopub.status.idle": "2025-08-18T03:36:06.493467Z", "shell.execute_reply": "2025-08-18T03:36:06.493113Z" } }, "outputs": [], "source": [ "def plot_it(ca_init):\n", " u.initial = ca_init\n", " n.finitialize(-65 * mV)\n", " n.continuerun(10 * ms)\n", " plt.plot(t, y)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And now let's run it a few times and plot it:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:06.495725Z", "iopub.status.busy": "2025-08-18T03:36:06.495102Z", "iopub.status.idle": "2025-08-18T03:36:10.476558Z", "shell.execute_reply": "2025-08-18T03:36:10.474481Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAksAAAGwCAYAAAC5ACFFAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAeLZJREFUeJzt3Xd8HOWBN/DfzGzvq95lyZZt3LuxnQoOpoQ3BHKBhISWkDc5CBBfCuQCJCHBAZK8PkogyYWSS8Akd4EUOB+cqQbbuGCMm2zZsmSrt+19Z94/ZneltSVZtiXtrvT7fj7zmZlnnp19Vra1Pz/zzDOCoigKiIiIiGhQYqYbQERERJTNGJaIiIiIhsGwRERERDQMhiUiIiKiYTAsEREREQ2DYYmIiIhoGAxLRERERMPQZLoBE4Esy2htbYXVaoUgCJluDhEREY2Aoijwer0oKyuDKA7df8SwNApaW1tRWVmZ6WYQERHRWTh+/DgqKiqGPM6wNAqsVisA9Ydts9ky3BoiIiIaCY/Hg8rKytT3+FAYlkZB8tKbzWZjWCIiIsoxpxtCwwHeRERERMNgWCIiIiIaBsMSERER0TAYloiIiIiGwbBERERENAyGJSIiIqJhMCwRERERDYNhiYiIiGgYDEtEREREw2BYIiIiIhoGwxIRERHRMBiWiIiIiIbBsJTFAnEZW12+TDeDiIhoUmNYylLHQxEsencfvvDBEbiisUw3h4iIaNJiWMpSFXotyg1aBGUFz7b1Zro5REREkxbDUpYSBAFfqSgEADzZ0oW4omS4RURERJMTw1IW+2yRE3laCSdCUbzS7c50c4iIiCYlhqUsZpREXFuaDwD49xPdGW4NERHR5MSwlOWuLy+AJADvuHzY7wtmujlERESTDsNSlqsw6HBpgQMA8FhzZ2YbQ0RENAkxLOWAW6uLAAAvdPThWDCc4dYQERFNLgxLOWC+1YRP5lkhA3i0ib1LRERE44lhKUfcUV0MAHi+vRetoUiGW0NERDR5MCzliOUOC863mxFVFDzKsUtERETjhmEph3y7pgQA8B+tPWji2CUiIqJxwbCUQz7itOLjTiuiioIHG9sz3RwiIqJJgWEpx/zr1FIAwF86+rCP8y4RERGNOYalHDPPasJnihxQAPzkSCsUPjOOiIhoTDEs5aA7a0qhFQS83uvFqz2eTDeHiIhoQmNYykE1Jj3+b2UhAODuwy0IxeUMt4iIiGjiYljKUd+qLkaJToumUASPH+dUAkRERGOFYSlHmTUSfjitDADwcFMHH4NCREQ0RhiWcthnihxY5bAgKCv4l4PHOdibiIhoDDAs5TBBEPCLmZUwigLecfnwh7aeTDeJiIhowmFYynFTjHp8r0ade+lHDa1o4XPjiIiIRhXD0gRwc2UhFtlM8MVl3HagGXFejiMiIho1DEsTgCQIeOS8KhhFEe+4fPgVH7RLREQ0ahiWJoipJgN+WlcOAHigsQ3vewIZbhEREdHEkFNh6a233sLll1+OsrIyCIKAF198cdj6f/nLX/CpT30KhYWFsNlsWLFiBf7nf/4nrc4Pf/hDCIKQtsycOXMMP8XY+UJpHj5daEdMAb627xh6o7FMN4mIiCjn5VRY8vv9mD9/Ph577LER1X/rrbfwqU99Ci+//DJ27tyJT37yk7j88svx/vvvp9WbPXs22traUsvmzZvHovljThAE/GJGJaYYdTgeiuCW/U0cv0RERHSONJluwJm45JJLcMkll4y4/vr169P277//fvz1r3/F3//+dyxcuDBVrtFoUFJSMlrNzCi7VoPfzanBp3cewuu9Xvy8sR3fqy3NdLOIiIhyVk71LJ0rWZbh9XqRl5eXVn748GGUlZWhtrYW1157LZqbm4c9TzgchsfjSVuyyWyLEQ/NqAQA/L+mDvy5vTfDLSIiIspdkyos/fznP4fP58PnP//5VNny5cvx9NNPY+PGjXj88cfR2NiIj370o/B6vUOeZ926dbDb7amlsrJyPJp/Rj5XkodbqooAAGsPHsfmvqE/DxEREQ1NUHL0GRmCIOCFF17AFVdcMaL6zz77LG6++Wb89a9/xerVq4es53K5UF1djV/+8pf4yle+MmidcDiMcLj/WWwejweVlZVwu92w2Wxn9DnGkqwo+Mb+Jvy10wWbRsTfF03HDLMh080iIiLKCh6PB3a7/bTf35OiZ2nDhg346le/ij/96U/DBiUAcDgcmD59OhoaGoaso9frYbPZ0pZsJAoC/m1mFZbZzfDEZHzxgyNoC3OGbyIiojMx4cPSc889hxtvvBHPPfccLrvsstPW9/l8OHLkCEpLJ8agaIMk4um5Nag16tESjuKq94+gIxzNdLOIiIhyRk6FJZ/Ph927d2P37t0AgMbGRuzevTs1IPuuu+7Cddddl6r/7LPP4rrrrsMvfvELLF++HO3t7Whvb4fb7U7V+fa3v40333wTx44dw7vvvovPfvazkCQJX/jCF8b1s42lPK0Gzy+YigqDFkeDYXxudwO6IgxMREREI5FTYWnHjh1YuHBh6rb/tWvXYuHChbjnnnsAAG1tbWl3sv3mN79BLBbDLbfcgtLS0tRy++23p+qcOHECX/jCFzBjxgx8/vOfR35+PrZu3YrCwsLx/XBjrNKgw38tmIYyvRaHA2H80+4j6I5w0koiIqLTydkB3tlkpAPEssHRQBifff8wOiIx1Jn02DB/KsoNukw3i4iIaNxxgPcEEI+HIcuj2/tTa9Ljvxb29zBdvusw6v2hUX0PIiKiiYRhKUspioIDB76LDz64CdGo+/QvOAPTTAb8bVEd6kx6tIajuGLXYexw+0f1PYiIiCYKhqUsFQgcRXfPa+jtewc7dl6FQKBxVM9fYdDhr4vqsMhmQl8sjs/tbsALHX2j+h5EREQTAcNSljKbp2Lxoj9Bry9FINCI7TuuRG/vO6P6HnlaDf68YCpW59sQktUJLH96pJUP3yUiIhqAYSmLWa3nYenSF2G3LUQs5sHuD25EU/NvoSjyqL2HWZLwzNwa3Jp4NMojzZ24/sNGuKO8U46IiAhgWMp6el0BFi36I0pKPgtFiaOh4WfYs+f/IhodvUtmkiDgB1PL8KtZ1TCIAv63x4MLd9RjJ8cxERERMSzlAlHUY9Z5D2HGjPsgijp097yGbe9dDpdrx6i+z5XFTvx1UR2qDTqcCEXxmfcP45GmDsi8LEdERJMYw1KOEAQBFeVfxJLF/wWjsRrhcBt27voCGhoehCyHT3+CEZpvNeHVpTPwmSIHYgrw06Nt+PzuI2gOjt57EBER5RKGpRxjtc7CsqV/RUnJFQBkNDX/Gu9t/ww8ng9H7T1sGglPzKrGL2dUwigK2Ozy4ZPb6/FMSzd7mYiIaNLhDN6jIFMzeHd2/Q8OHvwBotFeCIKEyoobUFNzGzQay6i9x9FAGN862IxtifFLqxwW/Gx6BerMhlF7DyIiokwY6fc3w9IoyOTjTiKRHtQf+iE6O18GAOh1xair+1cUFV0KQRBG5T1kRcHvTnTj/qOtCMoKtIKAr1UW4lvVxbBopFF5DyIiovHGsDSOsuHZcN09b+DQoR8hGFQfJJznXIXp0++B2Txt1N7jWDCMHxxuwf/2eAAAJTotfjitDJ8pcoxaMCMiIhovDEvjKBvCEqA+S66p+ddoanocshwBIKKs9HOoqb0dBn3JqL3PK91u3H24BU2hCABgqc2M708txQrH6F3+IyIiGmsMS+MoW8JSUiDQhMMN96O7+38BqFMPVFZcj+rqr0OrtY/Ke4TiMh4/3omHmzoQlNW/Qp/Ms+L7taWYazWNynsQERGNJYalcZRtYSnJ5d6JIw0PweXeDgDQaKyoqLgOlRU3QKfLG5X3aA9H8f+OteOPbT2IJf4mXV7owO3VRZjD0ERERFmMYWkcZWtYAgBFUdDT8wYajjwIv/8QAEAUjSgvvwZVlV+BwVA6Ku9zLBjGQ43t+EtHH5J/oT6ZZ8U3q4qxwmHmmCYiIso6DEvjKJvDUpKiyOjqehXHmn4Fr3cvAEAQtCgpvhwVFdfBZps7Ku9zwBfEw00d+GunC8kn2C22mfD1yiJcXGCHVmRoIiKi7MCwNI5yISwlKYqC3t7NONb0OFyubalyu30RKsq/jKKiiyGKunN+n6ZgGL9q7sSG9l6EE2OaSnRaXFeejy+V5qNIrz3n9yAiIjoXDEvjKJfC0kBu9/s4fuL36Oz8byhKFACg0xWitOSzKC29alSmHegMR/FkSzf+0NqD7mgMAKAVBFxaaMc1JXn4WJ4VEi/RERFRBjAsjaNcDUtJ4XAnWlo3oKXlWUQiXalym20BSkuvQnHRp6HVntvnCssyXupy48kTXdjhCaTKS3RaXFXixOdL8jCDs4ITEdE4YlgaR7kelpJkOYLuntfR1vYX9PS8DkWJAwBEUYeC/AtRVHQJCgo+CUk6t7vcPvQG8FxbL17o6ENfLJ4qn2814nPFebi00I5yw7lfCiQiIhoOw9I4mihhaaBwpBsd7X9Fa9t/pu6iAwBRNCA//xMoKroYBfkXQKMxn/17yDL+t8eDP7X3YlOPJzX1AAAssplwWaEDny60o9qoP5ePQkRENCiGpXE0EcNSkqIo8Pr2obPjZXR2/jeCoebUMVHUw+lcgfz8T6Ag/5MwGivO+n26IlH8tdOFf3S6sM3tx8C/lHMtRnyqwIYL82xYYDNxjBMREY0KhqVxNJHD0kCKosDn24+Ozv9GZ+fLCAab0o6bzXXIz/848vM+Drt9MSTp7HqEOsJRvNztxkudLrzr8qWmIACAPK2ET+TZcGGeFR/Ps6FApzmHT0RERJMZw9I4mixhaSBFUeD3H0JPzxvo7nkDbvfO1BgnQB3nZLctgtN5PpzOFbDZ5p3VlATdkRhe7XHjtR4v3uzzwBOT046fZzZgpcOCVU4LzndYkKdleCIiopFhWBpHkzEsnSwadaO3dzO6e15Hb+87iEQ6046LohEOxxI4HEthty2EzTYPGs2ZPXg3KivY6fFjU48Hr/V6sM8XSjsuAJhlUcPTcrsFi+0mlOo5UJyIiAbHsDSOGJbSKYqCQKARfa6t6Ovbgr6+rYhGe0+qJcJimQ6bbQHs9oWw2xbCZKqBIIgjfp+uSBRbXX686/LhnT4fDgVCp9Qp02uxyGbCYpsZi20mzLWaYJRG/h5ERDRxMSyNI4al4SmKDL//MPr6tsDl3gWP+32Ewq2n1JMkCyyWmbBaZ8FqmQ2rdRbM5mkjvnzXFYniXZcP7/b5sNMTwH5fEPJJdTQCMMNswGyLEXMtJsy2GDHbYoCdl++IiCYdhqVxxLB05sLhDrjdu+H27ILbvRte74eQ5fAp9QRBB4u5DhbrLFjM02E2T4PZPA16felpH87rj8XxgTeInR4/dnkC2OnxozMSG7RulUGHORYjZluMmG42oM6sR61RD53IXigioomKYWkcMSydO1mOIhA4Cq93H7y+/fB698Pn249YzDtofUmywGyeCrNpWipAGY1TYDRWDNkTpSgKWsJR7PMF8aE3qK59AZwIRQd/DwGYYtCjzqxHncmAOrNBXZv0sGikUfvsRESUGQxL44hhaWwoioJQ6AS83v3w+vbD72+A39+AYPAYFGXwHiJAhMFQBqOxCiZjNYymanVtrIbRWAVJMp7yClc0hn0+NTzt9QVx2B/G4UAIvvjJF/H6Feo0qDboMMWoR5VRXSf3i3Sa0/Z6ERFR5jEsjSOGpfEly1EEg02J8HRYXQeOIBhsQjweGPa1Wm0eDIYyddGXQW8oTW0bDGXQ6QogCCIURUF7JIrD/jAOBUI47A/hcEANUV1DXMpLMooCqox6lOu1KDfoUKbXokyvQ7lBXZfqtRxkTkSUBRiWxhHDUnZQFAWRSDeCwSYEg00IBJsQDDQhGGxGIHgMsZjntOcQBC0M+lLo9EXQ6wqh0xdCryuCTlcIvb4QOl0hglIBWmMmNIdiaAqG0RSM4FgwjKZQBC2hyCmDygeTp5VQrtehzKBFsU6LIp0WRTpNal2o16JQq4GBoYqIaMwwLI0jhqXcEI26EQq1IhRuRTjUltoOhdQlHO4ARhR1AEGQoNXmJwJUAbRaJ7TaPAiaPPSgEO1KPrpkKzrjRnTGtGiLCmgNx9ASiiIoj+w9AMCukVCk06AgEaQKtBrkaTVwaiXkJbbztBKciW32WBERjdxIv795vzRNGlqtHVqtHVbreYMel+UYIpFONThFOhEJdyIc6U6sOxGJdCMc7kQ02gtFiSMS6Txl8s0ka2KpTZUI0Gjs0OqcCGtK0CeWoVcoQg/y4VLscClm9Mom9MT16Ilr0B0TEVUAdywOdyyOw4FT7xQcjFEUEmGqP0TZNRKsGim1tkkibBopbbFrJJglkWOtiIgGwbBElCCKmtR4puHIchSRaA8i4S5EIl2IRHoRjfYiGu1DJNqHaLQP0WhvorwPsZgbgIJYzIVYzAWgEXYAdgA1Q7yHAiAAM1xwwAUnPEIBvGIJfGIefLDDCyu8sMCrmOBRDHDLesQhIiird/y1hAe/w2/Yzw+oYUojwaYRYZUkmCQRZkkNUv2LBLNGTBxT902iCLMmvY5JFKERGb6IKPflVFh666238NBDD2Hnzp1oa2vDCy+8gCuuuGLY17zxxhtYu3Yt9u3bh8rKSvzgBz/ADTfckFbnsccew0MPPYT29nbMnz8fjzzyCJYtWzZ2H4RymihqYdCXwKAvGVF9WY4hFnOpQSrSh2isD7GoF7GYB7GYB9GYB/GYF9GYB7FYf7k25oE51oJytKjpKZ5YBqEACMIIH2xqiIINPljhgxUBmBCA+aS1CUHBDD8sCMCEOCTI6O/JGi1aATCIgEEUoRcFGEURBkmEQZRglEQYUvuCui0mtqX0bWPi9VpBgE4UoRMF6AQBWlFIbIvqtiCkHdMKAnvLiOic5VRY8vv9mD9/Pm666SZceeWVp63f2NiIyy67DF//+tfxxz/+EZs2bcJXv/pVlJaWYs2aNQCA559/HmvXrsUTTzyB5cuXY/369VizZg3q6+tRVFQ01h+JJgFR1ECnK4BOVwCYz+y1ihJHLOZPBCg3YjE/4nE/4vEA4nE/YqntAOKJY7HEMXXpRDwWSNWT5eCAk/evItCdEqpCMCIEQ2oJJ5aTy9RtI8LQp7ZlQZ2HKqoA0TjgTU3DMHpBbKS0kKEVFGgEBTpBgTax6AQ1zGlFQCeo82ppBAGaAdv9ZULqmEYQB5SrvWepbUGARhQhCSI0AqAVpbRyjSBBIwhquShCgAiNKEIUREiCABFIrUVBgCictJ88LiBVJgmABAHCgLqnnAOJ1zA4Ep2VnB3gLQjCaXuWvve97+Gll17C3r17U2XXXHMNXC4XNm7cCABYvnw5li5dikcffRQAIMsyKisr8c1vfhN33nnniNoyFgO8FUVBNHrml1KIhqMocchyCHE5DDkegiyHIctByHIYcTkEWU6UpY6F0svlEOR4eMB+GLISgSJHISsRyHIUshxBRFYQUASEZAFhJbmIiECPKHSIQIcotIgktiOpMv2gx5KviUODWGKJQosYNKmyKDSIQQvlDJ4vOBmJkCFChgAFIpTEWk5sJ+oICpDYFxOpOllf3Vb3BWHAdqJ+/z4Sx5X0spPXwiCvG3guYZCyk14LRW1z8niKorYx2WYo6cfTtoXkcQWAkFifWjd1nkT9tNefXA+AoAzYUQYeT75P+vtjwM/41HX/G556PPlGQ7xeGebz4KSfW1Laa4QB73Oq9J/7MOccrJ6gpD578uwClPS3V4D7Pnol9Dr96c56RjjAG8CWLVuwevXqtLI1a9bgjjvuAABEIhHs3LkTd911V+q4KIpYvXo1tmzZMuR5w+EwwuH+Abcez+lvST9T0WgU999//6ifl+jcaBLLGXaRAQAUCIIMUZQhCDI0YhxaQYZVjEMQ4hBFPwTRA1GQIYpx9YtQkAdZD1Im9m8rAhAXRciCAFkUERfExH5yLSEmiIljkrqIiQghJNdCajsuSP3lA+rEBSmxHnhMLYsh+ZpEnQFrddFAhpCIKP1RRR50WxxQN/01aWUjDInJ153mj2qkf6RE526ov0cnJds7fT7o80Y3LI3UhA5L7e3tKC4uTisrLi6Gx+NBMBhEX18f4vH4oHUOHjw45HnXrVuHH/3oR2PSZqKJS4CiSIjHs+tRMWJiGbnTDCAb8jVKItANvg1BSfXWQJBTPQlCsisDif+BJ88nDPhffqJcAaAkunQUAIooQEluJ7qBZKjX8BQIUAQMOC6mukGUZNeJ0H9OIVGuCOolv9R7JVs6oIsk+Xol0dOhpLYT7Uh9ikR56nzp/ST9x4TEMWVA31OiXqL9yXrJ16d6KYSB+8KAvpdU0we8bgBhiPJUu/q7jE79rk//mfT3+QhpVYZ83VD7wkiOK/0/65PqpbV/mP1BXjLkz+JM66TalVZVOPnooOcSxaoRnX8sTOiwNFbuuusurF27NrXv8XhQWVk5qu+h1Wrx/e9/f1TPSZQrFFmBEpWhRONQwnHIkTiUiAwlktgOq9vJ42p5oiySeF1MVs+RXA/YzuoeEXWwEQR1kFH/WjpprQ5oSq8nicPXFwUIItQv1YHbAtTjp5SrwSq1nVYvua0eT5YrAiDHo4hFIohFI4jGwohFIohGw4hHo4jFIojH1GOxaEStFwsjHokimtyPhhGLRBGNhBCPRhCNhBGNhBGPhFMXatQRJEr/PhRA6b+MoySPKcnY1F93TP/4JA0kjQRRI0HSaCGIUqJMA1ESIWm0EDVqmSiKEEVRrSOKgCQmyiQIoghRkiBKEgRhYJmYOiYMqCtIIkQhWTZ0veT51HMKEERR/TMX1KlDBEEAxOS2+vcpfbv/GAR1SIwwWH1RSqyTdfvrIVlXwJDvgeSxAeFaq9WO+Z/fUCZ0WCopKUFHR0daWUdHB2w2G4xGIyRJgiRJg9YpKRn6Tie9Xg+9fmy7AgVBgE43+ANhibKVoiRCTigOORSDHIqltgcrk0NxKIkyOdwffJTIyCfuPJ30cSwCgPSeLUErDlgkCBoxvUwjAhp1LWgECJIIJNaCJKjHpJOOaQYeS+6LargZUJbcT4WaZOjIkHgsirDfj5Dfh5DPh7Dfh1DAj7DPh3DAj0goiEgwiEgwoK5DJ20HAoiGQlCU0fvzG4ooaaDR6aDRJhadFpJGC2nAtkang6QduK2DRquFpNVCo9Ul1snX6CBpNKn6okYDSZIgarSJdTLwqCFIHKQsFQZowpnQYWnFihV4+eWX08peffVVrFixAgCg0+mwePFibNq0KTVQXJZlbNq0Cbfeeut4N5coayiKAiUiQw5GIftjkANRyMHE2h/r3w4MWAejkENxID6K/3sXAEEnQdBJEPUSBJ04YFtK39YnjumkVMhJBSDtqQFI0EpqsJmAX27RSBhBjxtBj0ddez0I+rz9ASixhP3+tP1YeGSTn46IIEBnMEJnMqlroxE6gwEanR5avQEavR5afWJbp4fWYIBWp0+UGxLHEvuJ45rUWgdRzK7LuTSx5VRY8vl8aGhoSO03NjZi9+7dyMvLQ1VVFe666y60tLTg97//PQDg61//Oh599FF897vfxU033YTXXnsNf/rTn/DSSy+lzrF27Vpcf/31WLJkCZYtW4b169fD7/fjxhtvHPfPRzSWlKiMuC8C2RcdsI5C9kXUtT+xnQhH5xR6BEDQSxANGogGDQRDcluCYNT0byfWokEDQS9B0Gsg6kR1Oxl6JmCYOVPRcAj+vj74XX0IeAeGIHU74PWoZV43Ah73OYcendEEg8UCvdkCg8kMvdkCvckMnckIncGkBh/jgHUyDKXKjNDqDfyzowkjp8LSjh078MlPfjK1nxw3dP311+Ppp59GW1sbmpubU8dramrw0ksv4Vvf+hb+7d/+DRUVFfj3f//31BxLAHD11Vejq6sL99xzD9rb27FgwQJs3LjxlEHfRNlIURQo4Tji7jDinkhiCSPujiDuVQNRMgwp4bOY40gSIJo0EE1aiMbE2jRwrYGU2BaMiTqGRNDh7N3DUhQFIZ8X/r5e+F0u+F298PX1wu/qSyy9iYDUi0gwePoTnkSUNDDabDBZbTDa7DBYbTCYzTCYEyEobZ0ot1igN5nYa0NZIRqJw+8KI+AOI+SLoXZhYcbakrPzLGUTPkiXxoIiK5B9EcRcYTUMudUwJHvSg9EZje+RBEgWLUSLTl2btRCtOkgD1+b+QCTo2LNzNhRFQdDrgbe7C96ebnh7kuv+bV9vL+R4bMTn1Oj0MDscMNrsMNnsMFrtMNpsMFpt6r7NBqPVntrWGU38s6OsFYvE4esLw+cKw9cXgq93wHafug77+/99CALw9Uc/AXGUHxbOeZaIspwSV9ReoL4wYn0hxF3963hfCDFXeMSXwgSDBpJdB8mmg2TTJ9Y6iBYtJEv/WjBI/AIdBbIch6+3B+6Odrg62+Hp7EgFIU93F3w9PYhFIyM6l8FihdnhVBdnXtq2xeGEyeGExZnH8EM5JRKMwdMThKc7BE+3uvb2hlLBKOQf2aTLGp0Is0MPs12PaDgOvSkzk84yLBGNITkUQ6wnhFhPUF26Q4j1BhHvCyPuCQOn6xQSAMmuVxfbgDCUCEZiIhiJOl42GW0hvw/ujna4O9vh6miHp6sDrsS+p6trRL1CJrsD1vxCWPMLYC0o6N/OL4Q1Lx8mhxOaDN4OTXS24jEZ3p5EEOrpD0TqfjCtV2goGr0Eq1MPs0MPS54BFoceFqceFqchsdZDZ9RkxX8SGJaIzpEcjCHWnQxDwQHhKAT5dP97kgRIDj00ToO6dughOQ3QOPWQHAZIdp16yzmNiWgohL72VvS2nkBfWwv6WlvQ19YCV3sbQn7fsK8VJQm2wiLYi0pgLyqGraAoEYTUMGTJL2AQopymyAp8rjBcHYH+pVNde3tCON0gHoNFC1u+AbYCI2wFBljzjbA49bDmGWB26KE3ZUcQGgmGJaIRUBQFsieCaGcAsc4Aol1Bdd0ZgOwbPhCJFi00+UZo8g2pdTIQiRYdB0KPMUWW4enuRG9rC/paT6jrthb0trXA19M97GtNdgfsRcWwF5XAUVwCW1ExHEUlsBeXwJKXz4HQNCFEQjH0tQfQ1+ZHX0cA7mQo6gwiHh26+1ujExNByJgKRdZ8A+yF6lpnmDgRY+J8EqJRoCgK4n1hRNv9/cGoM4BYV3DYu8lEm64/EBWkByNRz39m40FRFPj6etBzvBndx5vQfbwJPSea0XO8GdFwaMjXGaw25JWWw1laDmdZOfLKyuEoKYOjqARag2EcPwHR2IqG4+hr96O3NbG0qWtv79D/PkRJgL3QCEexCY4ik7ouNsJeZILJpsuZnqFzxd/iNGnJoZgaitr9iLb5EW0PINruHzoUiVADUKEJ2qLk2gRNkZGBaJyFfD50NR1NC0Xdx5sQ9vsHrS9KGjhLy5BXVgFnaRmcZRXIK1MDktHKO1hpYpFlBa6OALpPeNFzwpcKRp6e0JCP+jHadMgrNcNZkgxEajiy5ulH/Q60XMTf8DThKYqCeG8IkRbfgGDkR7xviIn7JAHaQhM0JSZ1XZQIR/lG9dEUNG4URYGvtwedx46gs/Gouj52FJ6uzkHrC6IIZ0kZ8iurUFBZjYLKauRXVMNRUgpJw193NPFEI3H0tPjQfdyH7uNedJ/woeeED7EhLp8ZrVrklZrVpSyxlFpgsHB83XD424MmlIHBKNLiQ7TFh8gJH5TQ4HdmSHYdtCVmaEvN6rrEDE2hkYOqM0CRZfS1t6KzUQ1EnceOorPxCIJez6D1bYXFKKyegvyKqkQoqkJeWQU0fKYiTVDhQBSdTV50NXtTwcjVERh0oLVGJ6KgwoL8Civyy/rDkdHKfx9ng2GJclZyfFHkhFcNRYlFCQ4SjCQB2lIzdKUWaEtMqXAkmvi/qUzx9fag7cghtDccQntDPdqPNCASDJxSTxBF5JdXomhKLYpqpqKwuhZFU2phsFgy0Gqi8RGLxtF9wofOYx50HPOg85gXro5T/30A6iW0wkoLCiqsKKi0oLDSCluhESJvHhk1DEuUM+RwHJETXkSaPYg0eRE57oE82FweyWBUboGu3AptuQXaYhMvoWVQJBhAx9EGtDWo4ajtyKFB70TT6PQorJqCoppaFE2ZiqIptcivqoZWp89Aq4nGhywrcLUHEqFIDUc9LT7Ig0xKayswoKjahoJKCwoqrSiosMBs57+PscawRFlJURTEekKINHkQOe5FpMmDaLv/1MGJkgBtiRm6Cgu0yXDEYJRRiqLA3dGOlvr9aDm4D22H69Fz4jgUJX0MhSCIyK+sQum06SiZNh0lU6ejoLIaosTb8Wlii4Ri6DjmQfsRN9qOuNFx1I1I6NQbS4xWLYqm2FBUbUPxFBuKplhhtPAyWiYwLFFWUGIyIi0+hBvdiBzzDNlrJNn10FVboau0qesyC4NRhsnxODqPHUXLwf1oqd+HloP7EXC7TqlnLShE6VQ1GJVOm4Gi2qnQGYzj32CicebrC6PtiCsVjrpP+KDI6f/z0+hEFFXbEuHIiuIpNljzDZPm1vxsx7BEGSFH4mqPUaMb4aNuRI57oZx894ZGgK7cCl2lFbpqK/RVNkjsbs64SCiItkP1qWDUdrj+lHmMJI0GxVOno3zmLJRNPw+l06bD7HBmqMVE40dRFPS1B9B62IW2Bhfajrjh7Tl1HiOLU4/SqXaUTnOgZKod+WVm3qKfxRiWaFzI4RgixzwIN3rU3qMT3lMeEiuatdDX2KCbYoeuir1G2SIaDqG1/iCa932A4/v2oP3IYShyerDVm80onzELZTNmoXzmLJTU1vGuNJoUFEWBuzOIE/V9aD3Uh5ZDLgQ86Q9RFgQgv8KC0qkOlE61o2SqHdY8TniaSxiWaEwo0TjCxzwIH3Eh1OBCtMV3yngj0aaDvsYOfa0d+hq7ess+u5wzLhaNou3wQRzftwfH932ItsMHEY+lXxK1FRahfOZslM84D+UzZiG/ogqCyGBLE5+iKPB0B9FS70oFJL87PRxJWhEltXaUTbOjdKoDxbW2CfXoj8mIf3o0KhRZQbTVh1CDC+EGF8LH3EAsPR1JeQY1HNXYoK+xQ8rj9fhsIMfjaD9yGMf37UHz3g/QWn8AsWj6L39LfgGqZs9D5ex5qJo9D7bCogy1lmj8+d1hHD/QixMH+9BS3wffSRPaihoBJTV2lM9womKGA8VT7JC0/M/DRMKwRGclebdauKEP4cMuhI66T5nfSLTpYJjmgH6qumgcHG+ULdydHWja8z6OfbALzXs/QDiQ/pgQk92hBqM5akByFJcy2NKkEYvE0drgwvH9vTh+oBc9Len/PkRJQHGNDeXTnSif4URJjQ0aHe/inMgYlmjE5FAM4QYXQof6EDrUh7gr/X9Xgl6CfqpDDUjTHLyslkUioSCO7/sQxz7YhaY976OvrSXtuMFsQeXseaico/Yc5ZVX8s+OJg1FUdDb6kdzIhy1HnYhPvCGEwEoqrKiYmYeKmY4UTLVDq2e4WgyYViiISmKgmh7AKH6XoTq+xBp8gADb3eVBOirbdAnwpGu3ApB4hdsNlBkGZ3HjuLYnvfR9MEutNQfgBzv7/kTRBFl02eiet5CTJm/CMW10yCK/OVPk0fIH0Xzvh41IO3vPWVQttmhR9WsPFTOykPFTCfnN5rkGJYojRyKIXS4D6H6PoQP9SF+0i8QTYERhhlOGKY7oauxQ2TXc9aIBANo+nA3ju7ajqO7tp8y15G9qBhT5i9C9fxFqJo9D3qTOTMNJcoARVHQ0+JH095uNH3Yg/aj7rRnqmm0IsqmO1MByVliYu8qpTAsTXKKoiDWFURwfw9CB3sRafYAA3uftaJ6aS0RkDT5nEQwm7g62nF013s4ums7Tuz/MO2uNa3egMo58zBl/iJMmbcQjpIy/vKnSSUaiePEwT40fdiNpr09pwzMzi83o2p2Pqpm5aF0qoODsmlIDEuTkBJXED7mRuhAL4IHehA/acI0TaERhhl5MEx3Ql9jh8BfIFkjHouh9dCBVO9Rb8vxtOOO4lLULlqK2kXLUH7ebGi0fFAwTS6e7iCa9vbg2Ic9aKnvQzzW/78/SSuiYqYTU+bko3puAec6ohFjWJok5FAMofo+BA/0IFTfl37nmiTAMM0Bw8w8GGbkQcNfIFkl5Pehcdd2HNm1Hcc+2Imwv//OHEEUUTFzthqQFi+Ds7ScvUc0qSiygs5mLxp3d6FxTzd6W9PvXLPk6TFlbgGq5+SjYoaTd63RWWFYmsBivSE1HB3oRfioO21wtmjWwDAzH8bz8qCvc0LknR1ZxdvbjSPbt6Fhx1Yc37cHcrz/IZsGqw21CxajdvEyVM9bCIPZksGWEo2/eFxGa70LRz/oQuMH3fAPuDNXEAWU1NpSASmvzMz/QNA5Y1iaQJJ3rwX3diO0rxvR9kDacU2hEYZZakDSVdkgiPwFkk16Wo6jYftWNGzfgvaGQ2nH8iuqMHXJctQuWobSuum8c40mnUgohuZ9vWj8oAvHPuxBZEDvuEYvoXp2PmrmqwHJYOblZxpdDEs5TlEURE/4ENjbjdDebsQGjj8SAV21HcZZeTCclw9tAQdnZxNFltF+9LAakN7bgt7WE2nHS6fPxLQl56Nu2Qo4S8sz1EqizAl4Ijj2YTeO7u7CiQPp44+MVi1q5hWgZkEhKmY6odHyPxA0dhiWcpAiK4gccyO4twfBfd2ID3wukUaAoc4J45wCGGbmQeL/sLJKPBbDif17cXj7FhzZsRW+3p7UMVHSoGrOPExbugJTlyyHxZmXwZYSZYa3N4QjuzpxdHcX2o64054paSs0onZBIWrnF6C41g6RveM0ThiWcoQSkxE+4kJwXw+C+3og+6OpY4JOgmFmIiDNyOP4oywTj8XQvPcDHNq6GQ3btyLk86aOaQ1G1CxYjGnLVqB24RLOfUSTkqc7iCO7utCwqxOdxzxpxwqrrKhdUICa+YUcf0QZw7CUxeRIHOHDfWoP0oEeKKH+Qb6CUQPjrHwYZ+fDUOfk7f1ZJh6LovnDD1C/dTOObN+KkN+XOma02TF18XLULVuBqjnzodFxZmCafDzdQTTs7MSRXZ3obOr/DwQEoGyaA1MXFaJmfiFv76eswLCUpeKeMNof2gFlwPOJRIsWxtn5MM4pgL7WDkFiQMomqYC0ZTOO7EgPSCa7A3XLVmD6+R9BxXlzIErs/aPJx90VSASkLnQ19wckQQDK6hyYuqgItQsLYbbzoduUXRiWspRk00Ny6KFEZTUgzS3gHWxZKB6LounD3Ti05R007NiSNgeSGpBWqgFp1mzewUaTkqsjgIZdag9S9/H+/0AIAlA23Ylpi4tQu6AQJht7WCl7MSxlscKvzYNo0fIafZaJx6Jo2rNbHYO0Y+upAWn5Ksw4fxXKz2NAosnJ1aH2IDXs6kTPiQEBSRRQPt2RCkhGKwMS5QaGpSwm8RdJ1pDjcTTv/QAH330LDe9tQTjQH5DMDifqlqs9SOUzZzEg0aTk6wvh8I5OHN7ekX6JTRRQMdOJaYuKULOgAEYLf69R7mFYIhqCIstoPXQQB999E4e2voOA25U6ZnbmoW7ZSsw4/yMom3keAxJNSkFfBEd2deHw9g60NrhSt/kLooDKmU5MXVyE2vmFMFg4hQnlNoYlogEURUFXUyMOvvMmDr77FrzdXaljBqsNM85fhZkrP47ymbMgiBxgT5NPJBRD4+4uHNreiRMHeiEPeIxS6TQ76pYUY9riIl5iowmFYYkIQF9bCw6+8xYOvvNm2kzaWoMRdUvPx8xVH0fV3AWQNPwnQ5NPLBJH094eHN7egWN7exAfcJduYZVVDUhLinibP01YOfeb/7HHHsNDDz2E9vZ2zJ8/H4888giWLVs2aN1PfOITePPNN08pv/TSS/HSSy8BAG644QY888wzacfXrFmDjRs3jn7jKat4e7pR/+5bOPjuW+g42pAql7Ra1C5cipmrPoaaRUuh1fE2Zpp84nEZJw724fD2Dhzd3YXogHneHMUm1C0tRt2SIjhLOJEqTXw5FZaef/55rF27Fk888QSWL1+O9evXY82aNaivr0dRUdEp9f/yl78gEul/FEhPTw/mz5+Pf/qnf0qrd/HFF+Opp55K7ev1/HKcqAIeNw5vewcH33kLJw7uAxT1EoIgiqietxAzV34M05aez5m0aVJSZAVtR9w4vL0DDbs6EfL1PynA4tSjbkkx6pYWo6DSwrt0aVLJqbD0y1/+EjfffDNuvPFGAMATTzyBl156CU8++STuvPPOU+rn5aU/W2vDhg0wmUynhCW9Xo+SkpIRtyMcDiMcDqf2PR7PMLUp08KBABq2b8HBd99C0573ocj9lxDKZ87GzFUfx/TzV8Fks2ewlUSZoSgKupq9akDa2QlfX//vNqNVi6mLilC3tBiltXbO80aTVs6EpUgkgp07d+Kuu+5KlYmiiNWrV2PLli0jOsfvfvc7XHPNNTCb03sN3njjDRQVFcHpdOKCCy7AT37yE+Tn5w95nnXr1uFHP/rR2X0QGhfRSBiN7+/AwXfeROOuHYhF+3sYi2qmYuaqj2PGio/CVlCYwVYSZU5fux+HtnegYUcnXB2BVLnOIKF2YSHqlhajYoYTIp8UQJQ7Yam7uxvxeBzFxcVp5cXFxTh48OBpX//ee+9h7969+N3vfpdWfvHFF+PKK69ETU0Njhw5gu9///u45JJLsGXLFkhDPJLirrvuwtq1a1P7Ho8HlZWVZ/GpaDQlH1h78J030bB9CyLBYOpYXlmFGpBWfgx5ZeUZbCVR5nh7Qzi8vQOHd3SkzaYtaUVMmVuA6UuLUTUnDxotp8IgGihnwtK5+t3vfoe5c+eeMhj8mmuuSW3PnTsX8+bNw9SpU/HGG2/gwgsvHPRcer2e45qyhCLLaKnfj4PvvIVDWzcj6O2/JGotKMTMlR/DzFUfR2F1DcdY0KQU8ERwZJc6WWTbEXeqXBQFVM7KQ93SYtTML4DOMGm+DojOWM786ygoKIAkSejo6Egr7+joOO14I7/fjw0bNuDHP/7xad+ntrYWBQUFaGhoGDIsUWYpioLOxiM4+O5bqH/3bXh7+udCMtrsmLHiI5i58uMomz6TcyHRpBQOxnD0/S4c3tGBEwf7oCTnQhKAsmkO1C0txtRFhZxNm2iEciYs6XQ6LF68GJs2bcIVV1wBAJBlGZs2bcKtt9467Gv//Oc/IxwO40tf+tJp3+fEiRPo6elBaWnpaDSbRlFv6wl1ssh33kJfW0uqXGc0oW7ZSsxc9TFUzZkPcYjLp0QTWTQSx7E93Ti8vQNN+3ogx/oniyyqtqJuaTGmLS6GxclecaIzlTNhCQDWrl2L66+/HkuWLMGyZcuwfv16+P3+1N1x1113HcrLy7Fu3bq01/3ud7/DFVdcccqgbZ/Phx/96Ee46qqrUFJSgiNHjuC73/0upk2bhjVr1ozb56Khebq71LmQ3nkLnceOpMo1Oj1qFy/DzJUfRc2CJdDo+D9kmnziMRnHD/Ti8PYONH7QjWi4fy4kZ6kZ05cWYdqSYjiKTBlsJVHuy6mwdPXVV6Orqwv33HMP2tvbsWDBAmzcuDE16Lu5uRniSZdd6uvrsXnzZrzyyiunnE+SJOzZswfPPPMMXC4XysrKcNFFF+G+++7jmKQMCnjcOLRlMw6++yZaDu5PlYuSpM6FtOrjmLZkOXRGfgHQ5CPLCloPu3B4eweOvN+JsD+WOmbNN6TmQsovN3OcHtEoERRFUU5fjYbj8Xhgt9vhdrths9ky3ZyclJoL6Z030fTh7v65kAQBFefNxsyV6lxIRit/vjT5KIqCjmMeNGzvxOGdHQi4+6fCMNp0qFuszoVUXGNjQCI6AyP9/s6pniWaWKKRMBp3bcfBd97C0fe3Ix7tny24uLYOM1d9DDNWfhTWvIIMtpIoc3pafKlb/T3doVS53qRJzYVUPt0JkZNFEo0phiUaV/FYDM0f7lbnQtqxNX0upPJKzFz1Mcxc+TE4SzkXEk1O7q5gKiD1tvpT5RqdiJr5akCqmpUHScM7PYnGC8MSjTlFltFycD8OvvsmDm19J20uJFthEWasVAMS50KiycrvCqNhZycObe9A57H+fx+iRkD17HzULSnGlHkF0Op5pydRJjAs0ZhIzoV04J03Ub/lbfh6ulPHUnMhrfoEyupmcC4kmpRCviiOvK9OFtly2AUkp0ISgPIZTtQtLUbtgkIYzNqMtpOIGJZolHU3H0P91s2of/dtzoVEdJJIKIbGD9S5kI7v74Us999fU1JrT00WabbzblyibMKwROes58Rx1G95G4e2bkbPieZUuUarU+dCWvUxzoVEk1Y0HMexD7txZGcnju3tQTwqp44VVFpQt6QY05YUwZZvzGAriWg4DEt0VnpbT6gBactmdB9vSpVLGg2q5y/CjBUf5VxINGlFw3E07e1Bw84ONH3Yg9iAgGQvMqJuaTHqlhQjr9ScwVYS0UgxLNGI9bW34tCWzajf8ja6mhpT5aKkwZT5CzH9/I9g6pLlMJgtGWwlUWZEI3E0fdiDI7s6cezDbsQi/QHJVmDAtMXFmLa4CAWVFt7IQJRjGJZoWK6O9lQP0sDHjYiShKq5CxI9SOfDYGFAosknFomjaV8PGnZ24tiewQJSEaYtLmZAIspxDEt0CndnBw5t3Yz6LZvRcfRwqlwQRVTNma8GpKXnczZtmpRikTia9/WiYWcHGj/sQWzA89is+cmAVITCKisDEtEEwbBEAIDe1hYcfu9dHN72bnpAEkRUzp6LGSs/imlLV8Bks2ewlUSZEYvE0by/N9WDNPCBtdY8A6YmAlJRNQMS0UR0VmGpubkZTU1NCAQCKCwsxOzZs/ng2RyjKAq6jzfh8LZ3cHjbu2mDtAVBRMWsOZix4iOoW7YSJrsjcw0lypBIMIZje7tx9P0uNO3tSbvEZnHqU5fYiqYwIBFNdCMOS8eOHcPjjz+ODRs24MSJExj4/F2dToePfvSj+NrXvoarrroKIicZzEqKoqDjyGEceu9dHN72DlztbaljoiShcvY81C1biWlLz4fZ4cxgS4kyI+iNoHGPGpCOH+yFHOv/PWfJ02PqQrUHqXiKDQKfx0Y0aQjKwNQzhNtuuw3PPPMM1qxZg8svvxzLli1DWVkZjEYjent7sXfvXrz99tvYsGEDJEnCU089haVLl45H+7PCSJ9anAmyHEdr/QEc3vYuDr+3Bd6ertQxSavFlPmLULdsJWoXL4PRYs1gS4kyw9cXwtHd3Ti6uxOth1wY+BvRUWzC1IWFqF1YyDFIRBPQSL+/R9SzZDabcfToUeTn559yrKioCBdccAEuuOAC3Hvvvdi4cSOOHz8+qcJStonHoji+fy8Ob3sHDdu3IuB2pY5p9QbULFqK6ctXombBYs6DRJOSqzOAo+934ejuLnQ0etKOFVRaEgGpiPMgERGAEfYs0fCyoWcp5PehcfdOHNmxDY3v70AkGEgd05vNmLp4OeqWr0L1vAXQ6ji+jCYXRVHQ1exF455uNO7uQk+Lv/+gAJTW2lG7sBC1CwphK+BM2kSTxaj2LFF28nR1omHHNhzZuQ0n9n8IOd5/h47J7sC0JeejbvlKVM6eB0nDP2qaXGLROFrqXWjc041je7rhd4VTx0RRQPkMB2oXFqFmfgGfxUZEwxrxN+gFF1wwonqvvfbaWTeGhqcoCjobj6gBacfWtFm0ASC/ogpTlyzH1MXLUTptOgQOtKdJJuiLoOnDHjTu6Ubz/t60OZA0eglVs/JQM68AU+YVwGDWZrClRJRLRhyW3njjDVRXV+Oyyy6DVstfMuMlHovi+L4PUz1Ivp7u1DFBEFE+c5YakJYsh7OkLIMtJcoMV0cAjR90o3FPF9qPuNMGaJvtOkyZX4iaeQUon+GARitlrqFElLNGHJYeeOABPPXUU/jzn/+Ma6+9FjfddBPmzJkzlm2b1AIeNzY9+QSO7d6BSDCYKtfqDZgyfxGmLlmOmoVLOEkkTTpyXEZ7owfHPuhG455uuDoCaccLKi2YMq8ANfMKeAcbEY2KMx7gvWXLFjz55JP405/+hBkzZuCmm27CF7/4xay7ZX48jcUAbzkex+M3X4uQ3wezw4mpi9Xeo6o586HR6UblPYhyRcATQfP+HjTv7UHz/l6EA7HUMVESUD7Dmbq8Zs0zZLClRJRLRvr9fdZ3wwUCAfz5z3/GY489hv3796O1tXXSBqaxuhuufsvbsBUUoWRqHccf0aQiywo6mzxo3tuDpr096Gz2AgN+U+lNGlTPzUfNvEJUzcqDzsgbGIjozI353XC7du3Cm2++iQMHDmDOnDkcxzQGZqz4aKabQDRuQr4omvf3oGlfD5r39SLki6YdL6i0oHpOPqrnFKB4ihWixP9AENH4OKOw1NraiqeffhpPP/00PB4PvvSlL2Hbtm2YNWvWWLWPiCYoRVbQddyL5n1q71FHoydtcLbOIKFyVh6q5+SjanY+b+8noowZcVi69NJL8frrr+Oiiy7CQw89hMsuuwwazt1DRGfA2xvC8QO9OL6/FycO9iHkT+89yi83J3qP8lFca4fE3iMiygIjHrMkiiJKS0tRVFQ07N0lu3btGrXG5YpsmMGbKBtFgjG0HOrD8QN9OH6g95Q717QGCRUznKneIw7OJqLxNOpjlu69995RaRgRTVxyXEZnk1ftPTrQi46jHshy///HBAEorrGh4rw8VJ2Xh6IaG3uPiCjr8dlwo4A9SzRZKYqCvrYAWg714UR9H1rq+9Ju6wcAW6ERlYlwVD7DAb2JN4MQUXbgs+GIaNQpigJ3Z1ANRof60HLIhaAnklZHb9KgYoYTFeflofK8PNgL+WBaIsptZxyWenp6cM899+D1119HZ2cnZFlOO97b2ztqjSOizFIUBZ7ukBqM6tVwNPCBtAAgaUWU1NpRMcOBivPyUFRtgyhy1mwimjjOOCx9+ctfRkNDA77yla+guLiYjxIgmmC8vSE1GNX34cShPvh608ORqBFQUmNH+XQHymc4UVJjh6TluCMimrjOOCy9/fbb2Lx5M+bPnz8W7SGicaTICnrb/Gg74kbrYRfaGlzw9Z0UjkQBRVNsKJ+hhqPSWjs0Oj6QlogmjzMOSzNnzkRwwINdiSh3xGMyupq9aG1woa3BjbYjLoT96QOyBVFAYZUVFTMcKJ/uRMlUO3QGDm8kosnrjH8D/upXv8Kdd96Je+65Z9DHnPBuMKLsEQnF0H7UrQajBhc6Gj2IRdPHGWp0Iopr7CibZkdpnQPFU2wMR0REA5zxb0SHwwGPx4MLLrggrVxRFAiCgHg8PmqNI6KRS96p1t7oRvtRDzoa3eg54cPJk4MYLFqUTrWjdJoDZdMcKKiycK4jIqJhnHFYuvbaa6HVavHss89mZID3Y489hoceegjt7e2YP38+HnnkESxbtmzQuk8//TRuvPHGtDK9Xo9QKJTaVxQF9957L37729/C5XJh1apVePzxx1FXVzemn4PoXEWCMXQcU0NR+1EP2hvdp1xSAwBrvgFl0xwonaYGJGeJiTdmEBGdgTMOS3v37sX777+PGTNmjEV7hvX8889j7dq1eOKJJ7B8+XKsX78ea9asQX19PYqKigZ9jc1mQ319fWr/5C+JBx98EA8//DCeeeYZ1NTU4O6778aaNWuwf/9+GAx89AJlB0VW0NceQHujGx2NHrQfdaO3zQ+c1GskaUQUVVtRXGNDcY0dJbV2WJx8AC0R0bk447C0ZMkSHD9+PCNh6Ze//CVuvvnmVG/RE088gZdeeglPPvkk7rzzzkFfIwgCSkpKBj2mKArWr1+PH/zgB/jMZz4DAPj973+P4uJivPjii7jmmmvG5oMQDUNRFHh7Quhs8qKzyYPOJi+6mr2IBE/tNbIVGFBcY0dxjQ0ltXYUVFggaXhJjYhoNJ1xWPrmN7+J22+/Hd/5zncwd+7cUwZ4z5s3b9QaN1AkEsHOnTtx1113pcpEUcTq1auxZcuWIV/n8/lQXV0NWZaxaNEi3H///Zg9ezYAoLGxEe3t7Vi9enWqvt1ux/Lly7Fly5Yhw1I4HEY43H97tcfjOdePR5OUoijwu8KpYNTV5EVnkxchf/SUuhqdiKJqG0pqbamAZLaz14iIaKydcVi6+uqrAQA33XRTqkwQhDEf4N3d3Y14PI7i4uK08uLiYhw8eHDQ18yYMQNPPvkk5s2bB7fbjZ///OdYuXIl9u3bh4qKCrS3t6fOcfI5k8cGs27dOvzoRz86x09Ek1HAE+nvLUqsAyc9LgQARElAQYUFhdU2FFVbUVRtRV6pGSIHYhMRjbszDkuNjY1j0Y4xsWLFCqxYsSK1v3LlSpx33nn49a9/jfvuu++sz3vXXXdh7dq1qX2Px4PKyspzaitNLIqswN0dRPdxH7pPeNF9wofu475THhUCqPMa5ZWaUTTFiqJEOMovs3BWbCKiLHHGYam6unos2nFaBQUFkCQJHR0daeUdHR1Djkk6mVarxcKFC9HQ0AAAqdd1dHSgtLQ07ZwLFiwY8jx6vR56PS9/kCoWjaO31a8Go+OJYHTCh2h4kF5WAXAWm1BUbUNhtRXFU2zIr7BAyxmxiYiy1ojC0tatW3H++eeP6ISBQACNjY2pcUGjRafTYfHixdi0aROuuOIKAIAsy9i0aRNuvfXWEZ0jHo/jww8/xKWXXgoAqKmpQUlJCTZt2pQKRx6PB9u2bcM3vvGNUW0/TQwBTwQ9rb7+HqPjPvS1B6DIyil1JY2I/HIzCiosKKi0oqDCgvwKCyd8JCLKMSP6rf3lL38ZtbW1+OpXv4pLL70UZrP5lDr79+/HH/7wBzz11FN44IEHRj0sAcDatWtx/fXXY8mSJVi2bBnWr18Pv9+fujvuuuuuQ3l5OdatWwcA+PGPf4zzzz8f06ZNg8vlwkMPPYSmpiZ89atfBaCOtbrjjjvwk5/8BHV1dampA8rKylKBjCanSDCG3jY/elp86Gn1o7fVh95WP4LeUwdeA4DBrEVBZX8oKqi0wFls4hgjIqIJYERhaf/+/Xj88cfxgx/8AF/84hcxffp0lJWVwWAwoK+vDwcPHoTP58NnP/tZvPLKK5g7d+6YNPbqq69GV1cX7rnnHrS3t2PBggXYuHFjaoB2c3MzRLH/y6mvrw8333wz2tvb4XQ6sXjxYrz77ruYNWtWqs53v/td+P1+fO1rX4PL5cJHPvIRbNy4kXMsTRKxaBx97QH0tqrBqLfVj55WH3y9p44tAgAIgL3AmApEajiywuzQcaJHIqIJSlCUkx+GMLwdO3Zg8+bNaGpqQjAYREFBARYuXIhPfvKTyMvLG6t2ZjWPxwO73Q63281n42WpSCgGV0cAfe0B9LX7UwHJ3Rk45XEgSWaHHvllZuSVW9R1mRnOUjPHFxERTRAj/f4+q0kplyxZck6NIxoLiqIg4InANSAQJde+viF6igDoTRrkl1uQV2ZOrfNKzTCYtUO+hoiIJg+ONKWcE4/J8HQHT+kpcnUEEA6cOst1ktGqhbPEDGeJSV2XmpBfboHJxktoREQ0NIYlykrxuAxvdwiuzgDcnUG4OwNwdwXh6gzA2xMa8tKZIADWAmN/IBqwZk8RERGdDYYlyhg5LsPbG4IrEYYGrr09oUFvx0/S6ETYi0zIKzHBWWqGo9iEvFIz7EVGaLQcU0RERKOHYYnGVMgfhbcnBE93EJ7uxLonCHdXEN7uEOThApFWhL3ICHuRCY4iI+yFJtiLjHAUmWCy89IZERGND4YlOiexaDwRhpKBKAhPIhx5e0LDjiECAEkrwl6oBiB7oTEVhuxFJpjtOggiAxEREWXWiMLSwOegjdQPfvCDSTuVwEQSDcfh7Q3B1xtS131heHrUXiFPdxB+96kPgT2Z0aqFrcCoLvkG2ArVtb3IBItDz0BERERZbUTzLImiiBUrVkCn043opJs3b0Z9fT1qa2vPuYG5IFfnWVJk9VZ7byIIqaEonAhF6n7YP3zPEABo9RJsBQZY842wFxhhLTAkwpEBtnwjtHqOISIiouwz6vMsvfDCCygqKhpRXavVOtLT0hhRZAVBXxR+Vxg+V1hd96WHIV9fGHL89HOSag0SrHkGWPMMsOQZYM3TJ3qJ1EBksGg5foiIiCasEYWlp556Cna7fcQn/fWvf516BAmNvlgkngpAA8PQwP2AOzKiICQI6kzV1nwDLE5DIhTpE6FIDUd6I4e2ERHR5HXGjzuhU43FZbh4XMbBd9v6g1Bffyg63aDpFAEwWnWwOPQwJxZrnn5AD5EBZruOD3slIqJJacwed0LjQxQEvLXh0JC9QxqtmApAZoc+LRBZnOraZNdBYhAiIiI6J2cclkRRHHZ8SjweP6cGkUoQBUxfXgJJEgYNRHqThuOEiIiIxsEZh6UXXnghbT8ajeL999/HM888gx/96Eej1jACLrzuvEw3gYiIaNIbtTFLzz77LJ5//nn89a9/HY3T5ZRcnTqAiIhoMhvp9/eoDWg5//zzsWnTptE6HREREVFWGJWwFAwG8fDDD6O8vHw0TkdERESUNc54zJLT6UwbWKwoCrxeL0wmE/7whz+MauOIiIiIMu2Mw9L69evT9kVRRGFhIZYvXw6n0zla7SIiIiLKCmcclq6//vqxaAcRERFRVhrRmKU9e/ZAluURn3Tfvn2IxUY4yzQRERFRFhtRWFq4cCF6enpGfNIVK1agubn5rBtFRERElC1GdBlOURTcfffdMJlMIzppJBI5p0YRERERZYsRhaWPfexjqK+vH/FJV6xYAaPReNaNIiIiIsoWIwpLb7zxxhg3g4iIiCg78ZH0RERERMNgWCIiIiIaBsMSERER0TAYloiIiIiGwbBERERENAyGJSIiIqJhMCwRERERDYNhiYiIiGgYDEtEREREw2BYIiIiIhoGwxIRERHRMHIuLD322GOYMmUKDAYDli9fjvfee2/Iur/97W/x0Y9+FE6nE06nE6tXrz6l/g033ABBENKWiy++eKw/BhEREeWInApLzz//PNauXYt7770Xu3btwvz587FmzRp0dnYOWv+NN97AF77wBbz++uvYsmULKisrcdFFF6GlpSWt3sUXX4y2trbU8txzz43HxyEiIqIcICiKomS6ESO1fPlyLF26FI8++igAQJZlVFZW4pvf/CbuvPPO074+Ho/D6XTi0UcfxXXXXQdA7VlyuVx48cUXz7pdHo8HdrsdbrcbNpvtrM9DRERE42ek398507MUiUSwc+dOrF69OlUmiiJWr16NLVu2jOgcgUAA0WgUeXl5aeVvvPEGioqKMGPGDHzjG99AT0/PsOcJh8PweDxpCxEREU1MOROWuru7EY/HUVxcnFZeXFyM9vb2EZ3je9/7HsrKytIC18UXX4zf//732LRpEx544AG8+eabuOSSSxCPx4c8z7p162C321NLZWXl2X0oIiIiynqaTDdgvPzsZz/Dhg0b8MYbb8BgMKTKr7nmmtT23LlzMW/ePEydOhVvvPEGLrzwwkHPddddd2Ht2rWpfY/Hw8BEREQ0QeVMz1JBQQEkSUJHR0daeUdHB0pKSoZ97c9//nP87Gc/wyuvvIJ58+YNW7e2thYFBQVoaGgYso5er4fNZktbiIiIaGLKmbCk0+mwePFibNq0KVUmyzI2bdqEFStWDPm6Bx98EPfddx82btyIJUuWnPZ9Tpw4gZ6eHpSWlo5Ku4mIiCi35UxYAoC1a9fit7/9LZ555hkcOHAA3/jGN+D3+3HjjTcCAK677jrcddddqfoPPPAA7r77bjz55JOYMmUK2tvb0d7eDp/PBwDw+Xz4zne+g61bt+LYsWPYtGkTPvOZz2DatGlYs2ZNRj4jERERZZecGrN09dVXo6urC/fccw/a29uxYMECbNy4MTXou7m5GaLYn/8ef/xxRCIRfO5zn0s7z7333osf/vCHkCQJe/bswTPPPAOXy4WysjJcdNFFuO+++6DX68f1sxEREVF2yql5lrIV51kiIiLKPRNunqXJKHLsGGJdXZAjkUw3hYiIaNLKqctwk03jVZ+D7PcDAASjEZLdri42GySHHaLdDslm7y93qMdEux2S3QHJboNosUAQhAx/EiIiotzFsJSllHgcgkYDCAKgKFCCQcSCQcRGOAFniiSp4cpmg+hIhq0BAcueDFcnLTYbBJ1ubD4cERFRDmFYylKCJGH6tq1QZBmy14u4x4O4y4242w3Zo67VxZPalt0Dyj0eKKEQEI8j3teHeF8f0HSGbTCZ+nuyBvReiTZboucqWeaA5HSoa4cDIgfHExHRBMKwlOUEUUyFEpzhLOFyKJQIUy7InkSocqlBKu52JQKWJy1gqWHMo/ZmBQKIBQKItbWdWZuNxlRwkhz2AdsOaAZsD1xEqxWCyCF0RESUfRiWJjDRYIBoMEBbXHRGr1Piccg+30m9V65UkEr2cKnhy9W/uN1APN5/yfBMQlYyFA4RpgZdnA6IvFRIRERjjGGJTiFIUn9v1hlQZFkNWQMD1BBLLLXthhIIALLcf7nwTNpqMqV6r4bqtUoPWE4OeiciojPCsESjRhDF1GByVFWN+HVyJIJ43+kD1im9WLLcf6mwtQ3hkb6hRgPJ6YDG4YSUlwfJ6VT3nU5IDmdif0CZ0wnRaDyrnwkREeU+hiXKOFGng1hcdEaXC1MD34fttervvUpuK8EgEIsh3tWNeFf3iN9PMBgGBKi8QQNVf9BSe7l4NyER0cTAsEQ5KW3ge3X1iF8nh0JqcOrtRayvT+3RSlz+i7v6TimL9fUB0SiUUAixtjbE2kbegyVaLKlQlQpUg4Ws5GKzQZCks/uBEBHRmGFYoklFNBgglpRAW1IyovqKokD2BxB39fWHqr4+xHr7hg5ZLheQGL8l+3yIHj8+ssYJghoAk+EpzznspUGOvyIiGh8MS0TDEAQBksUMyWIGKipG9BpFliF7PP0BakDQSpX19qr7LnU/OV1D8nIhGhtH1sCB46+cyTFYHH9FRDSaGJaIRpkgiqm771Azstco0ag6FUPfyb1UfaeWJfaVQOCcxl8Ne2nQkV4maLVn98MgIpoAGJaIsoCg1UJTUABNQQFGOv95avzVyb1WfX2I9/Wqg91P6slSThp/NVKi1Tpg8Ppw46/Uni3Jbucko0Q0YTAsEeWoURt/dVLIGnT8ldcL2etFtLl5hI0T08dfDbw0OPBS4cDxV2Yzx18RUVZiWCKaJEZt/NVgdxKePP7qbCYZ1WrViUXPZIC7wXCWPw0iopFjWCKiIZ31+KuBlwCHHXvVi3hfYv6raBSxri7EurpG3j6jceRzXzmdkBwOjr8iojPGsEREo0rQaqEpLISmsHDEr5GDwdT4q6ECVXovlkud/yr5HMLWMxh/ZbONeOyVxumEaLNx/BXRJMewREQZJxqNEI1GaEtLR1RfHX/lH/ncV8nxV4oC2eOB7PEg2nQG469SlwdPnZZBk3dS6HI4OP6KaIJhWCKinKOOv7JAsliAysoRvUaJxxH3eIYZ3D6wTF1kn08df9Xbi3hv78gbqNGoA9wdjv71wO2BZY7+bY7BIspODEtENCkIkgRNYn6pkVIiEfVZg4ONveo9dXB7vLcXSjiszn/V04N4T8+ZtVGvHzpgnRSsBgYvjsMiGlsMS0REQxB0OmiLiqAtGvlDnlPzXw18iLO7/2HOg2673UAsBiUcRqyjA7GOjjNqp2g2jzhYpXqxOBaLaMQYloiIRtGZzn8FDBiDNTBguVyIuweELvfJYcudekyO7PdD9vsRbWkZeUMFQR3snlzsNog2+4BtmxqubPaT9m3qMwkZtGgSYVgiIsqwtDFYI5wDCxgwDsvlgux2I5ZYx12utO2Te7lkv18NWW43ZLcb0TNtsChCslohJsKTZLNBtNvUYJUMWwOClmRLBDG7jQ9/ppzEsERElKPOZhwWoI7FirvdatByexD3qL1UcbcnUe6G7Paoxz0eyB53op4HSiikDnp3uxE/l6DlsKeHK5utvxfLYoVks/avrVb1NVYrBL2eYYvGHcMSEdEkI+h0ZzwXVpIcDqu9U4kgldp2D9zvD1cDg5cSDp9b0II6j9fA8JQMVaLNCunktXVg0LJBslrUni1JOot3psmMYYmIiEZM1OshFhUBZzDoPSkVtIbq1UrMgRVPPIswbZ2YxkGJRs98KoeTP4PZrPZkWSzpa6ulP1RZbRAtZvVYcjFbIFnM6jxavANxUmFYIiKicXEuQUuRZciBgBqePF7IPq8arnw+de31Ie4dfq2EQgCQGhAfO4fPIuj1iRBlhmS2qAEsFawSISutfEDwMptTwUswmXhZMQcwLBERUdYTRDE1CH6kM72fTIlEEPf5Er1XPsjeIdbJ3i2/H7LPpwayxHYycCnhMOLhMOI9PWd1OTFFFBPhydwfvAaGqkRPlmgyJZbEdrLMbBpwzATBYGD4GgMMS0RENCkIOh00eXlAXt5Zn0OJxVIhKu7zQ/b7TgpUAwOWT933D1KWuKwIWYacuNw4KkQxLTylQpTZBMms9mSdEryGCV+i2cwABoYlIiKiEROSj7Kx23Euo5YURYESCiVCVzJUJYKX399f5vOplx+Ti9+fvp9YlEBAPbEsp4LZqBEEiEajGrSMRogGAwSTEaJBfaajaDJCSG4bDRCMRohGU/+2IVnHoJabkucwqescCGMMS0RERONMEIREqDCe1V2JJ1PiccjBEOSAH8rJYSoZsPxDlKct/XVTAUxRIAcCQCCA+Dm3dBDJn4XBoIYyYyJUJUNZYrvkh/dm7PmJDEtEREQ5TpAkSBYzJIt51M6pyDKUYLA/SAWDkINBtSwUghwIQg4GoARD6rFQEEog2L8dDKoBLvmaga8PBqFEIok3UqAEAoifJoyV3vfjUftsZ4phiYiIiE4hiCKExODzsZDsDVOCgVT4UkKnhio5GIISCmZ0ugaGJSIiIhp3yd4wjGJv2FjhkxCJiIiIhpFzYemxxx7DlClTYDAYsHz5crz33nvD1v/zn/+MmTNnwmAwYO7cuXj55ZfTjiuKgnvuuQelpaUwGo1YvXo1Dh8+PJYfgYiIiHJITl2Ge/7557F27Vo88cQTWL58OdavX481a9agvr4eRYPMCPvuu+/iC1/4AtatW4dPf/rTePbZZ3HFFVdg165dmDNnDgDgwQcfxMMPP4xnnnkGNTU1uPvuu7FmzRrs378fhgyNugfUEBeMBTP2/kRERNnEqDFmbIoBQVEUJSPvfBaWL1+OpUuX4tFHHwUAyLKMyspKfPOb38Sdd955Sv2rr74afr8f//jHP1Jl559/PhYsWIAnnngCiqKgrKwM//Iv/4Jvf/vbAAC3243i4mI8/fTTuOaaawZtRzgcRjgcTu17PB5UVlbC7XbDZrONymcNRANY/uzyUTkXERFRrtv2xW0waU2jek6PxwO73X7a7++cuQwXiUSwc+dOrF69OlUmiiJWr16NLVu2DPqaLVu2pNUHgDVr1qTqNzY2or29Pa2O3W7H8uXLhzwnAKxbtw52uz21VFZWnstHIyIiotOQ4+fyNL9zkzOX4bq7uxGPx1FcXJxWXlxcjIMHDw76mvb29kHrt7e3p44ny4aqM5i77roLa9euTe0ne5ZGk1FjxLYvbhvVcxIRjStFAZQ4EI8CchSQk9sxdYnH1PJ4VK0nxwYcjwLxeOJ1A+oqye3B6ibKlXjifIkleW5ZHrCdPCZDicegJF6nyDEochyIx6EocSAeS6xPPaeQeB8htR+DoMiAIkOQ4xAgZ/pPYEKRzmnO9HOTM2Epm+j1euj1+jF9D0EQRr27kYhynKKoASEeOWlJlMXCJx2PAvHwMHWS60FeFzvpdQPPJQ8MQCeFHlkNOEo8qoYJGlRcERCHiDgkxCFChggZQmJbgAIxtS0riTWERD31mJIoS9UbcB4ZIuTEeyiJ/fR6Q5xPGbwNCkQogghFEKAIEgARsiBCEJLlIpBcQwREERCkVPnARYAARUxspx0TUvuCIAKiAAgSIIgQBQHXiZmLLDkTlgoKCiBJEjo6OtLKOzo6UFJSMuhrSkpKhq2fXHd0dKB0wFOsOzo6sGDBglFsPRFNCLIMxIJqkIiFEkt4hOvhjkVGdg75nJ5vP66GGoYrKwKikBCHhBgkRBPrGCTElAHbqWMaxCEiqqjbMYiJtZQ6T3TA6+IQEYMEGSJiySCg9IeSZFmyTjy5KANfJyVCgUb9spYkdS1qIIgSIEqAoIEg9ZcJogYQJYiiBEHSAJIGYuKYKGnURZQgajQQRC00kgiNJEAjCpBEEVpJgCSq+xpJVNeiACmxLQnqcUkUIKb2AXFAuZQoFxPbWqG/rL9O/2uSa82Ac4pp5wE0oqhmlix/dttYy5mwpNPpsHjxYmzatAlXXHEFAHWA96ZNm3DrrbcO+poVK1Zg06ZNuOOOO1Jlr776KlasWAEAqKmpQUlJCTZt2pQKRx6PB9u2bcM3vvGNsfw4RDSaFEXt8Yj4gWgwsQTS17Hg0MdS69OUxUKZ/qRpZAiIQptYNIhCQgQaRBQNwooGEWjU8sR2aj+xRJRTy8LKSXUS9frLpP6woiTCCtJDTkyREE2GnER9QdJC0mih0Wig14jQSiJ0GjUI6BL7WklIrEe2nXy9VhKh1YgwJY5pJBG6AfU0kgBd6vUDzqFRtzVif3DRiGIqPBAl5UxYAoC1a9fi+uuvx5IlS7Bs2TKsX78efr8fN954IwDguuuuQ3l5OdatWwcAuP322/Hxj38cv/jFL3DZZZdhw4YN2LFjB37zm98AUJPyHXfcgZ/85Ceoq6tLTR1QVlaWCmRENMoURe0pifiBiC+xTm77BilP7IeHORbxqZeDxpEsaBAX9YiJOsREHaKCDhHoEBG0iChahKBFWNEiqGgQlLXwyxoEZA38cQ2CigZhRYcwtP2LokUYupP2+8v6A482FXDkEd6jo5UEGLRSYhFhTG5rJOi1IvQaNXjoNRJ0UnJbhClZpukvS67TXjNgXyeJ0GvTz6OTRIYPymk5FZauvvpqdHV14Z577kF7ezsWLFiAjRs3pgZoNzc3QxT7f3msXLkSzz77LH7wgx/g+9//Purq6vDiiy+m5lgCgO9+97vw+/342te+BpfLhY985CPYuHFjRudYIspachwIe4GwBwh5Ttr2nLTt7d8eeDzsHdtgI+mgaI1QNEbIGqMaaCQDoqIBEUGPsGBAEDqEFD0CihZ+RQefrIM3roU3poUnroErqi69UQ38shZB6BFSdAhCjyB0CEM34qAyHJ1GhEknwazTwKiTYNZJMOokGLUSHLpkmJES4UaEYcB2f/lJAeiUuhIkBhWic5JT8yxlq5HO00CUFeJRIOQGgn1A0KWuQ66h90Ou/rAT8Y1uWzRGQGdOLBZ1rVfXis6MqGRCWDAiIBgRgAF+RQ+PrIcnroMrZkBfTIueqBbdYS06wxp0hiS4wjLcwSii8dH71SYIgEWvgVWvgcWggUmngVkvwaTTwKTrX6thRz1m1Eow65MhKFkvUVcvwaSVoJFyZvYWoglppN/fOdWzREQnifgBfzcQ6FEXfzcQGLCfFn4SS8R77u8r6QGDDdBbAb0tsW07adua2o7rrPDDhD7ZiL64Hj1RHbpDEnpCMvr8EfT6o+gLRNDrj6CvLwJ3MApPMAp5xHlHARBNLP00ogCrQQ04Fr02FXYsibVV379t0WvUunrtSftq0JnsA1yJJjOGJaJsEvYBvg7A16mGnlOCUE+iPFF2Lo/E0dsBox0wOgGDQ10bHSftOwGDvT8AGeyA3gpZ1MEdjKLLF0a3N4wuXxhd3jC6fRE1/HQl1gF17QpGoCgRAK4zbqZOI8Ju1KYtNoOmf/vkYwO2GXKIaDQwLBGNtXgM8Hf1hyBf+4DtDsDb0b8f9Z/5+SU9YC4ATHmAqSCxna8uqcDjSA9Dehsgpf/zVxQFnlAMXd5k8Amjqzex9obR7etBl68V3d4Iun1hxEbe7ZNiM2iQZ9bBadYhz5RYm3VwmnTIM2vhTJQNDD8GrXTmPxMiolHEsER0LsI+wNMKeFoS6+R2Yt/XofYI4QyChc4CmAvVxZQPmPPVEGTKTwShggHl+Wr90/SexGUF3b4w2t0htLWF0OE5gXZPCO3uxJLYDkbPbBJBh0mLQoseBRY9Cq3qOt+SHn7yzDo4TDo4TFpoOUaHiHIQwxLRUGIRwH0ccDUB7hPpocidWIfdIzuXIALmIsBSBFhL1LWl+KQlUaa3nHFT3cEoTvQF0NIXxIm+IFpcQbS6gmhzh9DhCaHTG0Z8hD1BVoMmFXwKrXoUpoKQLq0836yHTsPwQ0QTH8MSTV7xGOA5Abiagb4mde1KrpvVMDSSHiG9DbCVDVjKE0tZIhgVqz1A4tlfTnIFIjjeG1QDkUsNROqi7ntDp78VXxSAIqsBJXYDSmyJdWK72GZAqV1dG3W87EVENBDDEk1ssTDQ2wj0HgF6jgA9DUDvUTUceVrUB2MOR2MEHFWAozI9BNnKAHsFYC1VBz+fI0VR0OuP4FiPH8e6A2jq8eNYT//aHTz9Yy7yzTqUO42ocBpR7jCizGFEqd2YCkQFFh1vVSciOgsMS5T75DjQd0wNQ72JQJTcdh3HsL1Dki4RhqoAR7W6dlb3b5sLTzse6Ez4wzEc6fKhodOHI12+VCBq6g7AGx6+d6jIqk+EIRPKHYlQ5DSi0qkGI5OO/5yJiMYCf7tS7ojHgL5GoOsg0HlQXXcdBLoPq09DH4rOCuTXAnlTgfxpQP5UwDlFDUSWYvXp2KOs1x9BQ6evf+ny4UinDy2uoW/1FwSgzG7ElAITqvPNmJKfXJtRlWfi5TEiogxhWKLsoyjqwOr2D4GOff3hqOew+rDUwWiMQF6tGoTypyaCUSIcjXLv0EDhWByHO3w40ObBgTYvDrZ7UN/uRY9/iHYCKLDoMLXQgmlFFtQUqGFoSoEJFU4Tb5MnIspCDEuUWbEI0F2vBqPUskd9HMdgtCagYDpQdB5QOAMoTKwd1WPSQ5SkKAq6vGHsb/PgYLs3EY48ONLlH/Ius3KHEXXFFkxLBKPk4jDpxqydREQ0+hiWaPzEY0DnfqBlp7q07VZ7jORBBi+LGjUIlcwBCmf2hyN71ZiGoiRXIIIPTrix57gLH5xw44MTLnR5B7/UZzdqcV6pFeeV2tSlxIapRWaOISIimiD425zGhqKot+G37ARadqnr1t2DP57DYAdK5gElc4HiOeq6cAag0Y9LU4OROD5scWPPiUQwOu5Cc2/glHqiAEwpMOO8UhtmldpSAanEZuAjNYiIJjCGJRod8SjQtgdofhdoehc4/p76DLOT6W1A2UKgfDFQvggonQ/YK8dsTNFgun1h7DjWhx3HerGjqQ97W9yDPrpjSr4J8yocmFdhx/xKB2aX2dhbREQ0CfE3P52dSAA4sR1o3qKGoxPbgehJvTGiVr2MVr44sSxRB1yPw2W0gZp7Ath6tAfbj/ViZ1Mfjnaf+vy1Iqse8ysdWFCphqN55Q7YTdpxbScREWUnhiUamXgMaN0FHHkdOPq6Go7kk+YFMjiAqhVA9Qp1XTIP0BrGvandvjDePdKDdxu68c6RbhzvPfXS34xiKxZPcWLpFCeWVOehwmnkpTQiIhoUwxINTlHUyR2PvqEGpGNvA2FPeh1b+YBwtFIdiD3OvUaAOuZo69EebG7oxjsN3TjY7k07rhEFLKh0YGlNHpZOcWJRlZN3pBER0YgxLFG/aAhofAs49N/A4VfVuY4GMjiA2o8DtZ8Eaj+hTuyYod6Y470BvF7fidcOduLdIz2IxOS04zNLrPjItAKsmlaApTV5sOj5V52IiM4Ov0EmO28HcPh/gPqN6uW1geOOJB1QuRyY+kk1IJXOP6eHwZ6LuKxgZ1MfNh3swOsHO3Gow5d2vNxhxEfr1HC0Ymo+CizjcycdERFNfAxLk5HrOLD/RWDfi0DLjvRjtnJg+hpg+iXAlI8AOlMmWggAiMVlvNfYi5f3tmHj3g50+/rnOZJEAYurnLjgvCJcMLMIdUUWjjkiIqIxwbA0WbhbgP1/Bfa9AJx4L/1Y2UI1HM24WB2UncHQEY3L2Hq0By9/2I5X9rWnPTbEZtDggplFuOC8Yny8rpB3qxER0bhgWJrIAr3A3v8CPvxP4PjWAQcEoHoVMPsKYOanAVtpploIQH2UyJ4Tbrzwfgv+9kEregcEJIdJizWzSnDpvFKsqM2HTjP+A8iJiGhyY1iaaOJRoGETsPuPQP1/pz9KpGoFMPtKYNb/AawlmWtjQosriBffb8Ffdp3Aka7+uY/yzTpcNLsEl80txfLaPGglBiQiIsochqWJousQsOsZYM+fAH9nf3npfGDeNWovkq0sY81LisZlbDrQgT9ua8bmhm4oiYmzDVoRF80qwZWLyvGRaQXQMCAREVGWYFjKZfEYUP8SsP3f1Vv+k8yFwLyrgflfUGfQzgItriA2vNeM57cfR+eAB9KuqM3HZxeV45I5JbAaOAaJiIiyD8NSLvJ2qL1IO54CvK1qmSAC0y8GFn4ZqPsUIGU+eCiKgneP9ODJzY14vb4TycevFVj0uHppBa5ZWoXKvMzdbUdERDQSDEu5pKseeOdhYM/z/WORTAXA4uuBxTcAjqqMNi8pGpfxjz2t+O1bjdjf1j/r96pp+fjismp8alYxB2oTEVHOYFjKBc3bgHfWA/Uv95dVLAOW3QzM+gygyY4JGD2hKJ7b1oyn3z2GNncIAGDUSvj8kgpcv3IKagstGW4hERHRmWNYylayDBx+RQ1JzVsShQIw8zJg1R1A5dIMNi6dOxjFk5sb8eQ7jfCG1IfrFlj0uHHVFFy7vIrPYSMiopzGsJSt3MeBDV8AFBkQtcD8a4BVtwMFdZluWcpgIamuyIKbP1aLzywog16TmUejEBERjSaGpWzlrFYHaxvswPn/nPGJIwcKRuJ48p1GPPHmkVRImlFsxe2r63Dx7BKIIh87QkREEwfDUjb7Pw9nugVp4rKCv+w6gV+8cgjtHnVMEkMSERFNdAxLNCJvH+7CT186gIPtXgBAucOI76yZgf8zv4whiYiIJjSGJRpWmzuI+/6xHy9/2A5AfZjtrRdMw3UrpsCg5ZgkIiKa+BiWaFDRuIyn3mnE+v89jEAkDkkU8OXzq3H7hXVwmnl3GxERTR45MzNgb28vrr32WthsNjgcDnzlK1+Bz+cbtv43v/lNzJgxA0ajEVVVVbjtttvgdrvT6gmCcMqyYcOGsf44WW1Xcx8ue/ht3P/yQQQicSyuduLvt34EP/w/sxmUiIho0smZnqVrr70WbW1tePXVVxGNRnHjjTfia1/7Gp599tlB67e2tqK1tRU///nPMWvWLDQ1NeHrX/86Wltb8Z//+Z9pdZ966ilcfPHFqX2HwzGWHyVrhWNx/L9XD+M3bx2BrAB5Zh3uvGQmPreoguOSiIho0hIUJfnc9+x14MABzJo1C9u3b8eSJUsAABs3bsSll16KEydOoKysbETn+fOf/4wvfelL8Pv90GjUnCgIAl544QVcccUVZ90+j8cDu90Ot9sNm8121ufJpA9PuPEvf96NQx1qb91nF5bjnk/PYk8SERFNWCP9/s6Jy3BbtmyBw+FIBSUAWL16NURRxLZt20Z8nuQPIxmUkm655RYUFBRg2bJlePLJJ3G6/BgOh+HxeNKWXBWLy1j/v4dwxa/ewaEOHwosOvz6y4vx/65ewKBERESEHLkM197ejqKiorQyjUaDvLw8tLe3j+gc3d3duO+++/C1r30trfzHP/4xLrjgAphMJrzyyiv453/+Z/h8Ptx2221DnmvdunX40Y9+dOYfJMu0u0O4bcP7eK+xFwBw2dxS3HfFHOQxJBEREaVkNCzdeeedeOCBB4atc+DAgXN+H4/Hg8suuwyzZs3CD3/4w7Rjd999d2p74cKF8Pv9eOihh4YNS3fddRfWrl2bdv7Kyspzbud4eqO+E2v/9AF6/RGYdRJ++tm5uGJheaabRURElHUyGpb+5V/+BTfccMOwdWpra1FSUoLOzs608lgsht7eXpSUlAz7eq/Xi4svvhhWqxUvvPACtFrtsPWXL1+O++67D+FwGHq9ftA6er1+yGPZTpYVPPzaYaz/38MAgFmlNjx27SLUFJgz3DIiIqLslNGwVFhYiMLCwtPWW7FiBVwuF3bu3InFixcDAF577TXIsozly5cP+TqPx4M1a9ZAr9fjb3/7GwwGw2nfa/fu3XA6nTkbhobjD8ew9k+78T/7OgAAXzq/Cj+4bBYnlyQiIhpGToxZOu+883DxxRfj5ptvxhNPPIFoNIpbb70V11xzTepOuJaWFlx44YX4/e9/j2XLlsHj8eCiiy5CIBDAH/7wh7SB2IWFhZAkCX//+9/R0dGB888/HwaDAa+++iruv/9+fPvb387kxx0Tx3sDuPn3O3Cw3QudJOInn52Dzy/JrUuHREREmZATYQkA/vjHP+LWW2/FhRdeCFEUcdVVV+Hhh/sfNBuNRlFfX49AIAAA2LVrV+pOuWnTpqWdq7GxEVOmTIFWq8Vjjz2Gb33rW1AUBdOmTcMvf/lL3HzzzeP3wcbBnhMu3PjUdvT4Iyi06vHElxZjcbUz080iIiLKCTkxz1K2y+Z5ll6v78Qtf9yFQCSO2WU2/O76pSixn/5yJBER0UQ30u/vnOlZojP35x3HcedfPkRcVvDRugI8/qXFsOj5R05ERHQm+M05QT39TiN++Pf9ANTZuB+4ah50mpyYg5SIiCirMCxNQL/b3Ij7/qEGpZs/WoPvX3oeBIHPdiMiIjobDEsTzG/fOoqfvqxO5PnPn5iK76yZwaBERER0DhiWJpB/f7s/KN12wTR861PTGZSIiIjOEcPSBPFfO0/gJy+pQen2C+vwrU9Nz3CLiIiIJgaO+J0AXjvYge/+1x4AwFc/UoM7VtdluEVEREQTB8NSjttxrBf//MddiMsKrlxYzsHcREREo4xhKYc19fjxlWd2IBSVccHMIjzwuXkQRQYlIiKi0cSwlKO8oSi++swOuINRzK904LEvLoJW4h8nERHRaOO3aw6Kywru2LAbhzt9KLbp8ZsvL4ZRJ2W6WURERBMSw1IO+sUr9dh0sBM6jYjffHkJim181hsREdFYYVjKMRv3tuFXbxwBADx41TzMr3RktkFEREQTHMNSDmlxBfHd/1SnCLj5ozW4YmF5hltEREQ08TEs5YhYXMbtz70PTyiG+ZUOfPfimZluEhER0aTAsJQjHt50GDua+mDVa/DINQt55xsREdE44TduDthypAePvN4AAPjplXNRlW/KcIuIiIgmD4alLOcLx/DtP38ARQE+v6QC/2d+WaabRERENKkwLGW5B/77IFpcQVQ4jbj38tmZbg4REdGkw7CUxbYc6cF/bG0CADxw1TyY9ZoMt4iIiGjyYVjKUoFIDN/7L3WagC8sq8KqaQUZbhEREdHkxK6KLOUKRJFn1iEal/H9SzlNABERUaYwLGWpMocR//WNlWjpC8Jq0Ga6OURERJMWL8NlMUkUOE0AERFRhjEsEREREQ2DYYmIiIhoGAxLRERERMNgWCIiIiIaBsMSERER0TAYloiIiIiGwbBERERENAyGJSIiIqJhMCwRERERDYNhiYiIiGgYDEtEREREw2BYIiIiIhoGwxIRERHRMDSZbsBEoCgKAMDj8WS4JURERDRSye/t5Pf4UBiWRoHX6wUAVFZWZrglREREdKa8Xi/sdvuQxwXldHGKTkuWZbS2tsJqtUIQhFE7r8fjQWVlJY4fPw6bzTZq56V0/DmPD/6cxw9/1uODP+fxMZY/Z0VR4PV6UVZWBlEcemQSe5ZGgSiKqKioGLPz22w2/kMcB/w5jw/+nMcPf9bjgz/n8TFWP+fhepSSOMCbiIiIaBgMS0RERETDYFjKYnq9Hvfeey/0en2mmzKh8ec8PvhzHj/8WY8P/pzHRzb8nDnAm4iIiGgY7FkiIiIiGgbDEhEREdEwGJaIiIiIhsGwRERERDQMhqUs9thjj2HKlCkwGAxYvnw53nvvvUw3aUJZt24dli5dCqvViqKiIlxxxRWor6/PdLMmvJ/97GcQBAF33HFHppsy4bS0tOBLX/oS8vPzYTQaMXfuXOzYsSPTzZpQ4vE47r77btTU1MBoNGLq1Km47777TvtsMTq9t956C5dffjnKysogCAJefPHFtOOKouCee+5BaWkpjEYjVq9ejcOHD49L2xiWstTzzz+PtWvX4t5778WuXbswf/58rFmzBp2dnZlu2oTx5ptv4pZbbsHWrVvx6quvIhqN4qKLLoLf78900yas7du349e//jXmzZuX6aZMOH19fVi1ahW0Wi3++7//G/v378cvfvELOJ3OTDdtQnnggQfw+OOP49FHH8WBAwfwwAMP4MEHH8QjjzyS6ablPL/fj/nz5+Oxxx4b9PiDDz6Ihx9+GE888QS2bdsGs9mMNWvWIBQKjX3jFMpKy5YtU2655ZbUfjweV8rKypR169ZlsFUTW2dnpwJAefPNNzPdlAnJ6/UqdXV1yquvvqp8/OMfV26//fZMN2lC+d73vqd85CMfyXQzJrzLLrtMuemmm9LKrrzySuXaa6/NUIsmJgDKCy+8kNqXZVkpKSlRHnrooVSZy+VS9Hq98txzz415e9izlIUikQh27tyJ1atXp8pEUcTq1auxZcuWDLZsYnO73QCAvLy8DLdkYrrllltw2WWXpf29ptHzt7/9DUuWLME//dM/oaioCAsXLsRvf/vbTDdrwlm5ciU2bdqEQ4cOAQA++OADbN68GZdcckmGWzaxNTY2or29Pe33h91ux/Lly8fle5EP0s1C3d3diMfjKC4uTisvLi7GwYMHM9SqiU2WZdxxxx1YtWoV5syZk+nmTDgbNmzArl27sH379kw3ZcI6evQoHn/8caxduxbf//73sX37dtx2223Q6XS4/vrrM928CePOO++Ex+PBzJkzIUkS4vE4fvrTn+Laa6/NdNMmtPb2dgAY9HsxeWwsMSwRQe312Lt3LzZv3pzppkw4x48fx+23345XX30VBoMh082ZsGRZxpIlS3D//fcDABYuXIi9e/fiiSeeYFgaRX/605/wxz/+Ec8++yxmz56N3bt344477kBZWRl/zhMYL8NloYKCAkiShI6OjrTyjo4OlJSUZKhVE9ett96Kf/zjH3j99ddRUVGR6eZMODt37kRnZycWLVoEjUYDjUaDN998Ew8//DA0Gg3i8XimmzghlJaWYtasWWll5513HpqbmzPUoonpO9/5Du68805cc801mDt3Lr785S/jW9/6FtatW5fppk1oye++TH0vMixlIZ1Oh8WLF2PTpk2pMlmWsWnTJqxYsSKDLZtYFEXBrbfeihdeeAGvvfYaampqMt2kCenCCy/Ehx9+iN27d6eWJUuW4Nprr8Xu3bshSVKmmzghrFq16pSpLw4dOoTq6uoMtWhiCgQCEMX0r05JkiDLcoZaNDnU1NSgpKQk7XvR4/Fg27Zt4/K9yMtwWWrt2rW4/vrrsWTJEixbtgzr16+H3+/HjTfemOmmTRi33HILnn32Wfz1r3+F1WpNXfe22+0wGo0Zbt3EYbVaTxkHZjabkZ+fz/Fho+hb3/oWVq5cifvvvx+f//zn8d577+E3v/kNfvOb32S6aRPK5Zdfjp/+9KeoqqrC7Nmz8f777+OXv/wlbrrppkw3Lef5fD40NDSk9hsbG7F7927k5eWhqqoKd9xxB37yk5+grq4ONTU1uPvuu1FWVoYrrrhi7Bs35vfb0Vl75JFHlKqqKkWn0ynLli1Ttm7dmukmTSgABl2eeuqpTDdtwuPUAWPj73//uzJnzhxFr9crM2fOVH7zm99kukkTjsfjUW6//XalqqpKMRgMSm1trfKv//qvSjgcznTTct7rr78+6O/k66+/XlEUdfqAu+++WykuLlb0er1y4YUXKvX19ePSNkFROO0oERER0VA4ZomIiIhoGAxLRERERMNgWCIiIiIaBsMSERER0TAYloiIiIiGwbBERERENAyGJSIiIqJhMCwRERERDYNhiYgmtfr6epSUlMDr9Y7Ze1xzzTX4xS9+MWbnJ6KxxRm8iWjC+cQnPoEFCxZg/fr1p6175ZVXYvHixfjXf/3XMWvP3r178bGPfQyNjY2w2+1j9j5ENDbYs0REk1ZzczP+8Y9/4IYbbhjT95kzZw6mTp2KP/zhD2P6PkQ0NhiWiGhCueGGG/Dmm2/i3/7t3yAIAgRBwLFjxwat+6c//Qnz589HeXl5quzpp5+Gw+HAP/7xD8yYMQMmkwmf+9znEAgE8Mwzz2DKlClwOp247bbbEI/HU6/71a9+hbq6OhgMBhQXF+Nzn/tc2ntdfvnl2LBhw5h8ZiIaW5pMN4CIaDT927/9Gw4dOoQ5c+bgxz/+MQCgsLBw0Lpvv/02lixZckp5IBDAww8/jA0bNsDr9eLKK6/EZz/7WTgcDrz88ss4evQorrrqKqxatQpXX301duzYgdtuuw3/8R//gZUrV6K3txdvv/122jmXLVuGn/70pwiHw9Dr9aP/wYlozDAsEdGEYrfbodPpYDKZUFJSMmzdpqamQcNSNBrF448/jqlTpwIAPve5z+E//uM/0NHRAYvFglmzZuGTn/wkXn/9dVx99dVobm6G2WzGpz/9aVitVlRXV2PhwoVp5ywrK0MkEkF7ezuqq6tH7wMT0ZjjZTgimrSCwSAMBsMp5SaTKRWUAKC4uBhTpkyBxWJJK+vs7AQAfOpTn0J1dTVqa2vx5S9/GX/84x8RCATSzmk0GgHglHIiyn4MS0Q0aRUUFKCvr++Ucq1Wm7YvCMKgZbIsAwCsVit27dqF5557DqWlpbjnnnswf/58uFyuVP3e3l4AQ18SJKLsxbBERBOOTqdLG3w9lIULF2L//v2j8p4ajQarV6/Ggw8+iD179uDYsWN47bXXUsf37t2LiooKFBQUjMr7EdH44ZglIppwpkyZgm3btuHYsWOwWCzIy8uDKJ76f8M1a9bgq1/9KuLxOCRJOuv3+8c//oGjR4/iYx/7GJxOJ15++WXIsowZM2ak6rz99tu46KKLzvo9iChz2LNERBPOt7/9bUiShFmzZqGwsBDNzc2D1rvkkkug0Wjwv//7v+f0fg6HA3/5y19wwQUX4LzzzsMTTzyB5557DrNnzwYAhEIhvPjii7j55pvP6X2IKDM4gzcRTWqPPfYY/va3v+F//ud/xuw9Hn/8cbzwwgt45ZVXxuw9iGjs8DIcEU1q//f//l+4XC54vV5YrdYxeQ+tVotHHnlkTM5NRGOPPUtEREREw+CYJSIiIqJhMCwRERERDYNhiYiIiGgYDEtEREREw2BYIiIiIhoGwxIRERHRMBiWiIiIiIbBsEREREQ0DIYlIiIiomH8f54XEg2fMSiJAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "for u0 in range(-4 * mM, 15 * mM, 2 * mM):\n", " plot_it(u0 * 0.1)\n", "\n", "plt.xlabel(\"t (ms)\")\n", "plt.ylabel(\"[u] (mM)\")\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here we see the stable fixed points at $0$ and $1$ and how all other trajectories (except for the one starting exactly at $\\alpha$, not shown) tend towards one of those values." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "
" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Alternatively, we might want to know how the kinetics change if we start at a fixed value (say 0.5) and vary $\\alpha$ instead. We can do this in the same way as above by setting the `value` property of $\\alpha$:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:10.508181Z", "iopub.status.busy": "2025-08-18T03:36:10.507722Z", "iopub.status.idle": "2025-08-18T03:36:10.514634Z", "shell.execute_reply": "2025-08-18T03:36:10.513113Z" } }, "outputs": [], "source": [ "def plot_α(α_init):\n", " u.initial = 0.5 * mM\n", " α.value = α_init\n", " n.finitialize(-65 * mV)\n", " n.continuerun(10 * ms)\n", " plt.gca().plot(t, y, label=f\"α = {α_init} mM\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here we have added a `label` argument to the plot and explicitly specify the `gca()` axes to allow displaying a legend; see the matplotlib.pyplot.legend documentation for more.\n", "\n", "Now we run it:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:10.517244Z", "iopub.status.busy": "2025-08-18T03:36:10.516938Z", "iopub.status.idle": "2025-08-18T03:36:13.547887Z", "shell.execute_reply": "2025-08-18T03:36:13.547471Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAtFJJREFUeJzs3XlcVPX+x/HXbMCw7yAIgiIqioooiFauiUuupZm7qJW37davMittMfOWZd5blmWuqbln5m7u+46KC4KAKMgm+z7MnN8fI6OkIiowgN/n4zGPmbPN+YDIvDnnu8gkSZIQBEEQBEGoI+TGLkAQBEEQBKEyiXAjCIIgCEKdIsKNIAiCIAh1igg3giAIgiDUKSLcCIIgCIJQp4hwIwiCIAhCnSLCjSAIgiAIdYrS2AVUN51OR2JiIlZWVshkMmOXIwiCIAhCBUiSRE5ODm5ubsjl5V+beeLCTWJiIh4eHsYuQxAEQRCER3Dt2jXq169f7j5PXLixsrIC9N8ca2trI1cjCIIgCEJFZGdn4+HhYfgcL88TF25Kb0VZW1uLcCMIgiAItUxFmpSIBsWCIAiCINQpItwIgiAIglCniHAjCIIgCEKdIsKNIAiCIAh1igg3giAIgiDUKSLcCIIgCIJQp4hwIwiCIAhCnSLCjSAIgiAIdYoIN4IgCIIg1Cki3AiCIAiCUKcYNdzs27ePvn374ubmhkwmY/369Q88Zs+ePbRp0wZTU1N8fHxYtGhRldcpCIIgCELtYdRwk5eXR6tWrZgzZ06F9o+NjaVPnz506dKF8PBw/v3vfzN+/Hi2bdtWxZUKgiAIglBbGHXizF69etGrV68K7z937ly8vb359ttvAWjWrBkHDhzgu+++IzQ0tKrKFARBEIQqI0kSknTr9a3l269L19/eh1vrJaQ7Xt/xXnfswz+Ov99+0u0dyx7P7eO55/qyX0cpE6UcZyuzB37tVaVWzQp++PBhunfvXmZdaGgo//73v+97TFFREUVFRYbl7OzsqipPEAShVpIkCa1OQqOVKNbq0Gh1lGglNFqdYVlTIqHR6dCU6NDcsa1EK1Gi06HVSZToJHSlz5JEifbWs07//uXvo0OrA61Od999dBKG59JAoF++1zr9su7WB7dOp19ful2i7HLp67vfs/R9bp+fO97jniHjHoHkfqHlXuToMEGDCSXI0aFEhwItSrQoZDr9M6XPWsN2Bbo7tmvvOK7s8bePvcd7yHTIkPTvxe3XciTk6JDfWi+/tU6BDrlMKruMjnwrb154/5eq+YGtgFoVbpKSknBxcSmzzsXFhezsbAoKClCr1XcdM2PGDD777LPqKlEQBOGxFZfoyC8uIa9YS15RCXlFJeTfel1YoqNIozU8F5XoKLzjuexrHUUlpc+3978dWHRodPqgUt6HbV0nQ4cZxZii0T/L9M+GdbLb20pfm96xzYQSTChBVfqQly5rUVGiDyqyO7Yb1pegumN96XsoZLX/H+OyJseo569V4eZRTJ48mXfeecewnJ2djYeHhxErEgShrioq0ZJdUEJ2oYbsAg3ZhSW3njVkF5SQdeu1PrBoywSY/CL96/ziEjRa43+4qRQyVAo5SrkME6UclaL0ISvzWnnrWS6ToZTLUJR56I8v3SaX32sfGQrZ3esMx916b4VchkqnwUSbg4m2AKU2HxOd/lmlLdA/SvJRaPNRagtQluSj1Obrn0vyUWgLUJTkoSjJR6HJR64tRK4t0j90xcb+dj+QJFeCXAlyhf5Zpl+WSpfv3PaPdTK5EkmuRPaP/W4vl65TgKz0WX7rtfz2a5m8gtvk+Fq6PPiLqkK1Kty4urqSnJxcZl1ycjLW1tb3vGoDYGpqiqmpaXWUJwhCHVKo0ZKeV0x6XjE384rJuPWcnldEep7m1nMxGfn6IJNVoKGoRFepNZgo5ViaKjE3UWBhokRtosBMJcdMpcBMqcBUJcdMqV9nqlJgptQ/myr1+5Q+m/1jnUohx0R5Z0iRY6KQozQEFxkymaxSvxYANIVQkA756frnomwozIbCrFuvs+54nX3HuluvSworv6Z/kqtAaQYqM/3zP18bltWgNAWVGhQm+tcKE1Cobj1Xwmu5Sv8sk9/336Oi/0pV8K9Zo9WqcBMSEsLmzZvLrNuxYwchISFGqkgQhNomt6iE5OxCkrMLSckuuvW6iJScW8s5haTmFJFfrH2k95fJwMpUibVahbWZCmu18tazChu1CiszJZam+oe5qRILEwXmJkosTMs+m5soUClq8FBk2hLIS4XcZMi/CQUZ+tCSf/N2gDG8ztC/1uRVzrlNLMHEAlTmt1+bmN96vt+20tcWoLLQh5J/BhWlGShq1ceicB9G/VfMzc0lOjrasBwbG0t4eDj29vZ4enoyefJkEhISWLJkCQCvvvoqP/zwA++//z5hYWHs2rWLVatWsWnTJmN9CYIg1CA6nURqbhHXM/K5nlHA9YwCEjJvPWfkk5RVSN5DhBalXIa9hUmZh4OFCXa3nu0tTLEzV2FjrjIEGCtTJXJ5Lf07WZL0ISU3RR9aDM/Jd6/Lv0nZvjIVJFOA2k7/MLO59bAGU+s7Xt96NrO5tf6O16ZW+tsfglAOo4abEydO0KVLF8NyaduY0aNHs2jRIm7cuEF8fLxhu7e3N5s2beLtt9/mv//9L/Xr1+fXX38V3cAF4Qmi1UkkZBQQk5ZLbFqe4XEtPZ/EzEKKtQ++NWRlqsTZ2hRnKzNcrE1xsTbD2doMZyv9aycrU+wtTLA2U1bN7RljkCT9LZ7sBMhKgOzrt54TIOu6/jk78eFu/cgUYOEIFk76sGJuD+YOoLa/x2t7/WtTa31bDUGoQjJJerLayGdnZ2NjY0NWVhbW1tbGLkcQhPvIKyrhcnIOkUk5XEnNJTYtn9i0XOLT88ttcKuQy3C1NqO+nRp3OzX17cypb6umvp0aVxszXKzNsDCto7ceinIh8ypkxEF6rP45Iw4y4/XhpTi3Yu+jtgNLF7B0vvV85+s71qntRVARqs3DfH7X0f/hgiDUFiVaHTFpeUQm6YPMpaQcIpOzuZZecN9jTJRyvB0s8Ha0wNvJAm8HCzwdzPUBxtoMZU1uq/K4CrMgLRrSLkN6zO0AkxGrbwPzIGp7sHEH6/pg7Xb7tY07WLvr1ylFJwyhdhPhRhCEaqPVScSk5nL2ehbnErI4ez2T84nZ9+1l5GRlShMXK3ycLWnkZIG3oyVejua42ahrb7uWipAk/ZWWtMuQFnXr+dbrnBvlH6u2AztvsPO6/bD1BBsPfXAxMa+GL0AQjEuEG0EQqkxKdiEnrmZw6moGZxOyOJ+Qdc8GvRYmCnxdrWjqakUTFyt8bz07WD4BVxDy0yE5ApLPQ1KE/nVaVPk9iyxdwbExODS6HWTsvcG2Aahtq6tyQaixRLgRBKFSSJLEldQ8TsSlczwugxNX07l6M/+u/dQqBS3crWlZ35aW9W3wd7fBy8Gibl+JAdDp4GY0JJ3VB5nSQJOdcO/95UqwbwiOvvog49jk1msffc8hQahh8jR5XM2+Snx2PGZKMzp7dDZaLSLcCILwyOJv5nMgOo2D0WkcjrlJel7ZkV5lMmjmak1gAztaeejDTCMnSxR1PchIEmRdg4RTkHASEk9DYjgU32dIetsG4OoPLs31D6dm+isxClW1li0ID6LVabmRd4PYrFjisuOIy4ozPKcUpBj2a+PcRoQbQRBqh4y8Yg5ducmB6DQORKfe1ejXVCmntYct7bzsaetlR5sGdlibPQEf0AWZcO3YrSBzSh9q8tPu3k9lDi4t9AHGtYX+tbOffhwXQahBirXFxGbFEp0ZzZXMK8RlxxGbFUt8djzF5UxXYW9mj6eVJ80dm1djtXcT4UYQhPuSJInolFz+vpjC3xeTORWfUWaCRaVcRhtPOzr6ONLRx4GW9W0xUdbhnkqlMq9B/BGIP6x/TrnAXQPayZX6EOMeCG5twL2N/taSGAFXqEG0Oi3Xcq4RnRlNVGYU0RnRRGdGczX7Klrp3gNequQqGlg3wMvaCy8brzLPNqY145ap+F8mCEIZGq2O43Hp/H0hhZ2Xku9qN9PExYqOPo481diBIG8HLOvqmDGlJAluXoHYvXD1kD7MZF+/ez/7hlA/SB9i3AP1V2VUZtVfryDcR1ZRFhfTL3Lx5kWiMqIMV2XudyXGSmVFY7vGNLRtiLe1N142Xnhbe+Nm6Yaiho8SXcd/KwmCUBElWh2HY26y8cwNtp5PIqtAY9hmopAT0siB7n4udGvqjJvtvSeprVOyb+jDTMxe/fM/G/3KFFCvFXiGgGd7/cPS2Ti1CsI/SJJEcn4yl9IvGcLMpfRL3Mi79zACZgozGtk2wsfWh8Z2jfGx9cHH1gdnc+daO0K3CDeC8ITS6iSOx6Wz8WwiW84lcfOOxsAOFiZ0aepM92YuPN3Yse6O6FuqKFcfYq7s1j+nXS67XWECHsHg9ZQ+0NRvq5+AURCMTJIkEvMSOZd2jos3bweZjKKMe+7vYeVBU/umNLVvSmNbfZBxt3JHLqtbt5Pr+G8sQRD+KTIph9UnrrHhTCIpOUWG9XbmKnq2qEffVvUI9nao+z2abl6BqO1weRtcPQjaOy/Ny8CtNXh3goadwKO9GPxOqBGyi7OJSIvgXOo5zqXpH+mF6Xftp5ApaGjbkGb2zWhq35Rm9s1oYt8EKxMrI1Rd/US4EYQnQFaBhg1nEllz4hpnrmcZ1luZKenZ3JXnWrnRoZEDqro8bUFJMcQfgsvbIWqbfsyZO9l5gU93aNhZf4VGbWeMKgXBoERXQmRGpCHInE09S1x23F37KWVKfO19ae7QnGYOzWhm3wwfWx/MlE9umy8RbgShjpIkiaOx6Sw/Gs/W80kU35riQCmX0a2ZM8+3qU+nJk6YKmt2w8DHoimA6L/hwga4vBWKsm9vkyuhQQdoHAqNe+gHyqul7QuEuiFfk8+Z1DOcTjnNqZRTnE09S0HJ3XOs1besj7+TPy0dW9LCsQXNHJphqngCRvN+CCLcCEIdk1tUwh+nE1h6+CqRybcHjWvqasULgfUZEOCOY12e1qAoR3+r6eIGiNoBmjt6e1k464OMbw9o2EWMLyMYVWp+KqdTThvCTGR65F3dr61MrGjp1BJ/R3/Dw85MXFV8EBFuBKGOiE7J4bfDV1l7KoHcohJAP9XBgAB3hgV50sLdutb2fHigoly4tAnO/wFXdoH2dlsibDzBrx806wf124G8Dt96E2q01PxUjiUd41jSMY4nHedazrW79nGzcCPAJYA2zm0IcA6gkW2jOtfYtzqIcCMItZgkSRyLTeeXfTHsvHR76POGThaMbN+A5wPr190RgrUafe+mc6v0webOKzQOPvow49cP6rUWt5sEo8gszOR48nGO3jjKsaRjxGbFltkuQ4avnS8BzgG0cdGHGVcLVyNVW7eIcCMItZBWJ7H9fBJz98Vw5lomoP/87t7MhdEhXnT0caibV2kkST/FwdlVELG27BQH9o3AfzD49QfnZiLQCNUuT5PHiaQTHE06yrEbx4jMiCyzXYaMpvZNCXINIqheEAHOAU9M76XqJsKNINQixSU61p66zs97rxB3a+RgE6WcFwLrM+Hphng71tGxV7ISIHw5nFkO6TG315s7gv8L4D9EPzKwCDRCNZIkiciMSA4kHOBQ4iFOJ5+mRCops4+PrQ/tXNsR7BpMW9e2NWZ6grpOhBtBqAVKQ80Pu6JJyNT3nrBRqxgV0oBRIV44WdXBBsJajb5h8KklEL0DJH1vL1Tm0LQPtHxR321bzJwtVKP0wnQOJx7mUOIhDiYc5GbhzTLb3S3daV+vPcH1gmnn2g5HtaORKn2yiXAjCDXYvUKNk5UprzzTkJeCPOvmyME3r8Dp3/RXanKTb69v0BECRkKzvmBqabz6hCeKVqflXNo59ifs51DCIc7fPI90xySpaqWaINcgOrp3pKNbRzytPY1YrVCqDv5mFITaT6eTWB+ewLfbL5cJNa92asTwYE/MVHVsbBqdVn+V5tgvELP79noLJ2g9TB9qHBsbrz7hiZKvyefwjcPsubaHfdf33TUCsK+dryHMBDgHYKIwMU6hwn2JcCMINcy+y6nM2HKJizf0A87V6VCTn66/7XR8PmTF31opg8bPQptR4NtT3HYSqkVyXjJ7r+9lz7U9HL1xtMxM2VYqKzq4d+Ap96fo4NYBZ3MxSWpNJ8KNINQQEQlZ/GfLJQ5E63sAWZkp+VdnH8Z08EJtUsdCTWI4HJsHEWugpFC/Tm0HbUZD2zCwa2DU8oS6r7Qx8O743ey+tpuL6RfLbHe3dKeLRxc6e3SmjUsbVHIRsmsTEW4EwciSswv5assl1p1OAMBEIWdkSANe7+KDnUUdutyt0+kbBh/6HuL2317v2hKCX4EWz4NKbbz6hDpPkiQi0iLYEb+DHXE7uJ573bBNhoyWTi3p7NGZzvU708i2Ud0cTuEJIcKNIBhJcYmOBQdj+X5nFHnF+iHX+7d2490eTfCwr0MzUJcU6celOfwDpF7Sr5Mr9ePRBL0CHkGiC7dQZXSSjjOpZ9get52d8Tu5kXfDsM1UYUoHtw508ejCM/WfwUHtYMRKhcokwo0gGMGeyBQ+/+sCMWl5AAR42vJZv+a0rG9r3MIqU0EGnFgAR3++3evJxAoCR0P7iWBT37j1CXVWia6EU8mn2HF1Bzvjd5JakGrYplaq6VS/E90bdOdp96cxV9WhPyQEAxFuBKEaXUvP5/ONF9hxQf9h72hpyge9mjIowB25vI5cvchL01+lOTYPinP166zcoP2rEDgGzMQgZkLlkySJM6ln2By7mW1x28r0cLJSWdHZozPdG3Sng1sHzJRmRqxUqA4i3AhCNdDqJBYejOXb7Zcp0GhRyGWM6eDFW90b1525n3KS9O1pTiy4Pc+Tc3Po+CY0HwTKOtR+SKgxLmdcZnPMZrbGbSUhN8Gw3sbUhq4eXXm2wbO0r9celeh190QR4UYQqtiFxGw+WHeWs9ezAAjytueLAS3wdakjc8pkXYcDs/Vduktn43YLgGfehya9RHsaodJdy7nGltgtbIndQnRmtGG9udKcbp7d6N2wN8H1gkUPpyeYCDeCUEUKNVr+tzOKn/fFoNVJWJkp+bB3M15s61E3bkFlJ8K+mXDqN9Bp9Os8gvWhxqebCDVCpUovTGdL7BY2x2zmbNpZw3qVXMUz9Z+hl3cvOtXvJG45CYAIN4JQJU7EpfPemrPE3mow3NvflU/7NsfZug784s1LgwPf6dvUlF6p8XoannkPvJ8RoUaoNBqthn3X9/HnlT/Zf32/YVJKuUxOkGsQvb17061BN6xNrI1cqVDTiHAjCJWoqETLdzui+GXfFXQSuFqb8Xn/5vRo7mrs0h5fYRYc+gGO/Hi7obBnCHSdAl4djVubUGdIksSF9Av8Gf0nW2K3kFmUadjW3KE5fRv1JdQrVExIKZRLhBtBqCQXErN5Z1U4l5JyAHghsD5T+/rV/gbDxflw7Gd9u5rCTP26eq31oUbcfhIqSWp+KhtjNrLhyoYy7Wic1E481/A5+jXqh4+djxErFGoTEW4E4TGVaHX8vC+G2X9fRqOVcLAwYcYg/9p/tUang7MrYOc0yEnUr3NsAl0/1s/MLUKN8Jg0Og17r+1lXdQ6DiYeRCfpADCRm9DNsxv9fPrRvl57lHLxUSU8HPETIwiP4XpGPv9eEc6JqxkA9PBz4ctB/jhamhq5sscUsxe2fwRJ5/TLNp7Q9SPwHwzyOjbPlVDt4rPjWRu1lj+j/+Rm4U3D+tZOrenn049Qr1DRjkZ4LCLcCMIj2hpxg/fXnCW7sAQrUyWf9mvOoDbutXs+mtRI2D4Forbpl01t4Jn/00+ToKoDjaEFoynSFrHz6k7WRq3lWNIxw3oHMwcG+AxggM8AvGy8jFegUKeIcCMID6lQo2X6pov8duQqAK09bPn+pYDaPR9UbgrsmQEnF4Ok1c/91HYcdJoEFmK+HeHRXcm8wprLa/gr5i+yivRjPcmQ0dG9Iy80foFnPJ4R49EIlU6EG0F4CNEpuby+/JSh0fCrnRrxfz18USnkRq7sEZUUw9GfYO9MKNZ/TTR9Drp/Bo6i8abwaApKCtgWt421l9cSnhpuWO9i7sKgxoMY6DOQepb1jFegUOeJcCMIFfTH6et8uC6CAo0WBwsTZr3Ymk6+TsYu69FF74Qtk+BmlH7ZLQB6TBfduoVHFp8dz8rIlayPXk92cTYACpmCTvU78bzv83R064hCtNkSqoEIN4LwAMUlOr7YdIElh/W3oTr6OPDdkNa1d0C+jKuw7UO4tFG/bOEEz34OLYeCvJZegRKMRqvTsj9hPysureBg4kHDendLd17wfYH+jfrjZF6L/wgQaiURbgShHElZhfxr2UlOxWcC8Ga3xrzVrTGK2jh9gqYADv4PDsyCkkKQKSD4Fej8gZipW3ho6YXp/BH1B6svrzZMWClDxlPuTzG06VCecn8KuUyEZcE4RLgRhPs4EnOT15efIi23GGszJd+92JpuzVyMXdbDkySI3AJbP4BM/dUnvJ6GXl+Di59xaxNqFUmSOJd2jhWXVrAtbhvFumJAPwP3QJ+BDPEdgoe1h5GrFAQRbgThLpIkMf9ALDO2XEKrk2jqasXPIwNp4GBh7NIeXtZ12PweRG7WL1u5QegX0HyQGIRPqLAibRFbYrfw+6XfuXDzgmF9c4fmDG06lJ5ePcWElUKNIsKNINyhUKPlwz/Ose6U/jL7wAB3vhzoj9qkljWC1Gn1E1vumqafB0quhJDX9ZNbmloauzqhlkgrSGNV5CpWRq4kvTAd0I8e3NO7J0ObDMXfyd/IFQrCvYlwIwi3pOYU8cpvJzgVn4lCLmNKn2aM7uBV+wblu3EW/noLEk/plz2C4bnZ4haUUGGX0i+x9MJSNsduRqPTAPpu3EObDuX5xs9jZ2Zn5AoFoXwi3AgCcD4xiwmLT5CYVYi1mZIfhwfyVONaNutwcZ5+IL7DP+oH4jO1ge6fQOBY0QtKeCCtTsve63tZenEpx5OOG9a3dGrJSL+RdPPsJgbbE2oNEW6EJ97WiCTeXhlOgUZLQycLfh3VloZOtezWTdQO2PgOZMXrl/0GQK+vwKqWT94pVLnc4lzWR69n2cVlXM+9DujHpunRoAfD/YbTyqmVkSsUhIcnwo3wxJIkiR/3XGHmtkgAnm7syA/D2mCjrkV/neanw9bJ+tm7AWw8oM+34Btq3LqEGu96znWWXVzGH9F/kKfJA8DaxJrBvoMZ2nQorhYiGAu1lwg3whOpRKtjyp8R/H7sGgBjOnjxcZ9mKGvTNAqXNsHGtyE3GWRyCJ4IXT4UDYaFcp1LPcfC8wvZGb8TnaQDwNvGmxHNRvBcw+cwV9XiOdIE4RYRboQnTl5RCa8vP8XuyFTkMvi0X3NGhXgZu6yKy0+HLe/DudX6ZUdf6P8jeLQzbl1CjaWTdOy/vp+F5xdyMvmkYX0Htw6M9BtJB7cOYsA9oU4R4UZ4oqTkFDJu0QnOJWRhppLzv6EB9Gheiy6/X9gAm96BvFT91ZoOb0LnyaASY4wIdyvWFrMpZhOLzy/mStYVAJRyJb29ezOm+Rga2zU2coWCUDVEuBGeGNEpuYxZeIzrGQXYW5gwf3RbAjxrSZfWvDTY/C6c/0O/7NRUf7WmfqBx6xJqpOzibFZHrmbZxWWkFqQCYKGyYLDvYIY3Gy7a0wh1ngg3whPheFw64xefIKtAg5eDOYvGBuHlWEtGHL64UT9uTX6afj6op96GTu+D0tTYlQk1TFJeEr9d+I01l9eQX5IPgLPamRF+I3jB9wWsTKyMXKEgVA8RboQ6b2vEDd5cEU5xiY4AT1t+HdUWB8taEAyKcvTzQZ1eql929oMBP4JbgHHrEmqcyPRIFp1fxNbYrZRIJQD42PowpvkYenv3RqWoRT0ABaESiHAj1Gkrj8czed05dBI86+fC/4YG1I6pFK4ehj9euTXRpQw6vgldPhJXa4QyTiWf4pdzv3Aw4aBhXZBrEGOaj+Ep96dq3+jaglBJRLgR6qxf9l3hy82XABjazoPpA/1RyGv4L/uSYtjzJRyYDUhg4wkD54JXR2NXJtQQkiRx5MYRfjn7CyeSTwAgl8l5tsGzjG0+luaOzY1coSAYn9H7/s2ZMwcvLy/MzMwIDg7m2LFj5e4/e/ZsmjRpglqtxsPDg7fffpvCwsJqqlaoDSRJ4uutlwzB5pVODZkxqBYEm+QLMK8rHPgOkKD1cJh4UAQbAdD/XO+9tpcRm0fw8o6XOZF8AqVcyfONn2fjgI180+kbEWwE4RajXrlZuXIl77zzDnPnziU4OJjZs2cTGhpKZGQkzs7Od+2/fPlyPvjgAxYsWECHDh24fPkyY8aMQSaTMWvWLCN8BUJNo9VJTPkzguVH9dMQTOrZlImdGxm5qgfQ6eDIj7DzM9AWg7kD9P0vNOtr7MqEGkAn6fj76t/8cvYXIjP0o2mbKkx5wfcFxjQfI3o+CcI9yCRJkox18uDgYNq1a8cPP/wAgE6nw8PDgzfeeIMPPvjgrv1ff/11Ll68yM6dOw3r/u///o+jR49y4MCBCp0zOzsbGxsbsrKysLa2rpwvRKgRikt0vLMqnI1nbyCTwZcD/XkpyNPYZZUvK0HftiZuv365cSj0+x6sXIxbl2B0JboStsRuYd65ecRmxQJgrjTnxaYvMspvFI7qWjaxqyA8pof5/DbalZvi4mJOnjzJ5MmTDevkcjndu3fn8OHD9zymQ4cOLF26lGPHjhEUFERMTAybN29m5MiR9z1PUVERRUVFhuXs7OzK+yKEGqOgWMurS0+y93IqKoWM715szXMt3YxdVvkuboQNr0NBBqgsIHQ6BI4B0Qj0iVasLWbDlQ3MPzffMJGllYkVw5sNZ3jT4dia2Rq3QEGoBYwWbtLS0tBqtbi4lP0L1cXFhUuXLt3zmGHDhpGWlsZTTz2FJEmUlJTw6quv8uGHH973PDNmzOCzzz6r1NqFmiWvqISwRcc5GpuOmUrO3BGBdG5y923NGkNTANs/huO/6pfrtYYXFoBDDb99JlSpwpJC1katZWHEQpLzkwGwM7VjVPNRDG0yFEsTMWeYIFRUreottWfPHr788kt+/PFHgoODiY6O5q233mLatGlMmTLlnsdMnjyZd955x7CcnZ2Nh4dHdZUsVLGcQg1jFx7nxNUMLE2VLBrbjrZe9sYu6/6SL8DacZByQb/c4U3oOgWUJsatSzCaPE0eKyNXsvj8YtIL0wH9wHtjWozh+cbPi4ksBeERGC3cODo6olAoSE5OLrM+OTkZV9d7N5CbMmUKI0eOZPz48QD4+/uTl5fHyy+/zEcffYRcfnfnL1NTU0xNxdggdVFWgYbRC44Rfi0TKzMlv40LprWHrbHLujdJghMLYNuHUFIIFs76Lt4+3YxdmWAkWUVZLL+0nKUXlpJdrL9d7m7pTliLMPr79MdUIX5vCcKjMlq4MTExITAwkJ07dzJgwABA36B4586dvP766/c8Jj8//64Ao1DoB2QzYrtowQgy84sZMf8oEQnZ2JqrWDoumBbuNsYu697y02HDG3Bpo37ZpzsM+Aksa/CtM6HK3Cy4yW8XfmNF5AryNHkAeFl7Md5/PL0b9kYlF6MJC8LjMuptqXfeeYfRo0fTtm1bgoKCmD17Nnl5eYwdOxaAUaNG4e7uzowZMwDo27cvs2bNIiAgwHBbasqUKfTt29cQcoS672ZuESPmH+PijWzsLUxYNj6YZvVqaM+3uAOw7mXITgC5Crp/Cu3/Bfe4yijUbcl5ySw6v4g1l9dQqNWPzdXYrjEv+7/Msw2eRSEXv8MEobIYNdy8+OKLpKamMnXqVJKSkmjdujVbt241NDKOj48vc6Xm448/RiaT8fHHH5OQkICTkxN9+/Zl+vTpxvoShGqWklPIiF+Pcjk5F0dLU5ZPCMbXpQZOBqgtgX1fw76ZIOnAvpG+0bBba2NXJlSzhNwEFpxbwB/Rf6DRaQBo4dCCl1u+TCePTshlIugKQmUz6jg3xiDGuam9UrILGTrvCDGpebhYm7J8QnsaOdXAHiTZN2DteLh6a+yl1iOg11dgWgNrFapMbFYsv577lU0xm9BKWgDaOLfhlZavEOIWIuZ9EoSHVCvGuRGEh5GaU8SwX48Sk5qHm40Zyye0x8vRwthl3e3KLlg7AfLTwMRSP9Kw/wvGrkqoRpHpkfx67le2xW1DQv+3Ywe3Dkzwn0Bb17ZGrk4Qngwi3Ag13s3cIkb8epTolFxcrc34/eX2NHCoYcFGp4U9/9HfhkIClxYweDE4+hi7MqGaRKRF8MvZX9h9bbdhXWePzrzs/zL+Tv5GrEwQnjwi3Ag1WkZeMcN/PUpkcg7OVqY1M9jkJOvHrimdQiFwDPT8D6jURi1LqB4nk0/yy9lfOJR4CAAZMkK9QhnvP54m9k2MXJ0gPJlEuBFqrKx8DSPmH+VSUg6Olvpg413TbkXF7NW3r8lL0U+h0Hc2tBxi7KqEKiZJEodvHOaXs79wMvkkAAqZgj4N+zDOfxwNbRoauUJBeLKJcCPUSFkFGkYuOMr5xGwcLEz4fUJwzWo8rNPqb0Ht+Q8ggbOf/jaUk6+xKxOqkCRJ7Lm2h3nn5nEu7RwAKrmKAT4DCGsRRn2r+sYtUBAEQIQboQbKKdSPPHz2ehZ25iqWT2hP45rU3Ts3BdZNgJg9+uWAkdDrazARw+TXVVqdlh3xO5h3dh6XMy4DYKYw4wXfFxjTfAwuFmIWd0GoSUS4EWqU3KISxiw8Tvi1TGzUKpaOD6aJaw0KNrH79e1rcpNBZQ59ZkHrl4xdlVBFNDoNW2K3MO/sPOKy4wCwUFkwtMlQRvqNxEHtYNwCBUG4JxFuhBqjoFhL2KLjnLyagbWZkmXjg2nuVkOmVNDp4OB3sOsL/aB8Tk31t6Gcmxq7MqEKFGuLWR+9ngURC0jITQDA2sSaEX4jGNZ0GDamNeTnUhCEexLhRqgRikq0vLL0JMdi07Ey1U+CWWPmiirIhD9ehctb9MuthkGfb8CkhjVuFh5bQUkBay+vZeH5haTkpwBgb2bP6OajebHJi1ioxL+5INQGItwIRlei1fHvFeHsu5yKWqVg4dh2tKops3snnYOVIyEjFhSm0HsmBI42dlVCJcstzmVl5EqWXFhCemE6AM7mzoS1CGNQ40GolaJbvyDUJiLcCEal00lMWnuOLRFJmCjk/DIqkLZe9sYuSy/8d9j4bygpBFtPGLIE3AKMXZVQibKKslh2cRlLLy4lpzgHAHdLd8b7j6dfo36YKEyMXKEgCI9ChBvBaCRJ4rO/zrP21HUUchnfDwvg6cZOxi4LSopg6wdwYoF+2ac7DJoH5jUkdAmP7WbBTZZcWMKKSyvIL8kHwNvGmwn+E+jl3QulXPxqFITaTPwPFozmm+2RLD58FZkMvhncktDmrsYuCTKvwapRkHgKkEHnD+CZ90EuZm6uC5Lyklh0fhFrL6+lUFsIQBO7JkxoOYHunt1RyBVGrlAQhMogwo1gFD/uiWbO7isAfDGgBQMDasDgZ9E79aMNF6SDmS08Px8adzd2VUIluJZzjQURC1gfvZ4SXQkALR1b8nLLl3mm/jNihm5BqGNEuBGq3ZLDcXy9NRKAD3s3ZXhwA+MWpNPB/m9h93RAgnqt9e1r7Ixcl/DYYrJimH9uPptiNqGVtAC0c23HBP8JtK/XXoQaQaijRLgRqtWak9eZ+ud5AN7s6sPLzzQybkEFGbe6eW/VL7cZBb1mgsrMuHUJjyUyPZJ55+axPW47EhIAHd078rL/y7RxaWPk6gRBqGoi3AjVZsu5G7y/5gwAYzt68fazRp6H6cZZWDUSMuL03bz7fAttRhq3JuGxhKeEM+/cPPZd32dY19WjKy+3fJnmjs2NWJkgCNVJhBuhWhyMTuOtFeHoJBjStj5T+vgZ95bA6WWw6Z1b3bwbwIu/Qb1WxqtHeGSSJHE06Sjzzs7jWNIxAOQyOaENQhnfcjy+dmIyU0F40ohwI1S5c9ezeHnJCYq1Onr7uzJjUEvkciMFG00hbJ0EJxfplxv3gEG/gNrOOPUIj0wn6dh7bW+ZGbqVciX9GvUjrEUYDaxFmylBeFKJcCNUqdi0PMYsPEZesZaOPg5892JrFMYKNpnxt7p5nwZk0OVDePpd0c27ltHqtGy/up155+YRlREFgKnClOcbP8/YFmNxtagBQwoIgmBUItwIVSY5u5CR849yM68Yf3cbfh7ZFlOlkcYRif77VjfvDP1Vmud/1Q/OJ9QaGq2Gv2L+Yv65+cTnxANihm5BEO5NhBuhSmQVaBi94BjXMwrwdrRg4dh2WJoa4cdNp4N9M2HPDEDST58wZIl+OgWhVigsKWRt1FoWnV9EUl4SADamNoxoNoKXmr4kZugWBOEuItwIla5Qo2X84uNcSsrB2cqUJWFBOFqaVn8h+emw7mWI3qFfDhwLvb4CpRFqER5abnEuKyJX8NuF3wyTWTqpnRjdfDSDfQdjrjI3coWCINRUItwIlapEq+P15ac4HpeBlZmSxWFBeNgb4UMoMVzfzTszHpRm0GcWBAyv/jqEh5ZZmMnSi0tZfml5mcksw1qE0d+nP6YKEU4FQSifCDdCpZEkicnrzvH3xRRMlXLmj25Hs3rW1V/Iqd9g0/+BtgjsvGDIb1CvZfXXITyUpLwkllxYwprLaygoKQD0k1mO9x9PL+9eqOQqI1coCEJtIcKNUGm+2hrJ6pP6Gb5/GNaGIO9qnkVbUwhb3oNTS/TLvj1h4FzRzbuGu5J5hQURC9gcs5kSST/vUzP7ZkxoOYFunt2Qy0RvNkEQHo4IN0Kl+HV/DHP36ifCnDHQn2f9XKq3gIyr+m7eN8IBGXT9CJ76P9HNuwY7nXKaBecWsOf6HsO6INcgxrYYS0e3jmLeJ0EQHpkIN8JjW3fqOl9sugjApJ5NGdLOo3oLiPob1pV287aHF+ZDo67VW4NQITpJx/7r+1kQsYBTKacAkCGjm2c3wlqE4e/kb+QKBUGoC0S4ER7L7kspvL/mLADjnvLm1U4Nq+/kd3XzbnOrm3c1hyvhgTQ6DVtit7AwYiHRmdEAqOQq+jXqx+jmo/G28TZyhUJNodVq0Wg0xi5DMBITExPklXDFXYQb4ZGdvJrBxGUnKdFJDAxw56PezarvVkJBhr6bd9R2/bLo5l0j5WvyWRu1liUXlhjGqLFQWTCkyRBGNBuBs7mzkSsUagpJkkhKSiIzM9PYpQhGJJfL8fb2xsTE5LHeR4Qb4ZFcTs4hbNFxCjU6Ojdx4usXqnG+qBtnYOVIyLwqunnXUBmFGSy/tJzfL/1OVlEWAA5mDozwG8GQJkOwNjFCLzqhRisNNs7Ozpibm4s2V08gnU5HYmIiN27cwNPT87F+BkS4ER5aQmYBo+YfI6tAQ4CnLT8Ob4NKUU0Nd08v1XfzNszmvVR0865BruVc47cLv/FH1B8UagsB8LTyZEyLMfRr1E+MUSPck1arNQQbBwcxjcaTzMnJicTEREpKSlCpHn34BxFuhIeSnlfMqPlHScouxMfZkgWj22FuUg0/RnfN5h0Kg34W3bxriPCUcBafX8yua7vQSToA/Bz8CGsRRnfP7ijkRppTTKgVStvYmJuLUaefdKW3o7RarQg3QvXIKyph7KLjXEnNw83GjCVhQdhZPN590Qq5azbvj+Bp0c3b2LQ6LTvjd7L4wmLOpp41rO/o3pExzccQ7Bosbi0ID0X8vAiV9TMgwo1QIcUlOiYuO8WZa5nYmqtYMi4IN1t11Z84euet2bzTxWzeNUS+Jp8/ov/gtwu/kZCbAOh7Pj3X8DlG+Y3Cx87HyBUKgvCkE+FGeCCdTuLd1WfYdzkVtUrBwjHt8HG2quqTwv5vYfd0QIJ6rfXdvO0aVO15hftKzktm+aXlrL682jDnk42pDS82eZGXmr6Eo9rRyBUKgiDoiev6QrkkSeLzjRfYcCYRpVzG3JGBBHhWcTuXggxY8RLs/gKQoM1oCNsmgo2RRKZH8uH+D+m5ricLIhaQU5xDA+sGfBz8MTte2MEbAW+IYCMI1Sg+Pp4+ffpgbm6Os7Mz7733HiUlJZV+nk8//RSZTEbPnj3v2jZz5kxkMhmdO3eu9PNWBnHlRijXj3uusOhQHADfDG5FJ1+nqj3hjbP62bwz4kBhCs/NgoARVXtO4S46SceBhAMsubCEozeOGta3cW7D6Oaj6ezRWcz5JAhGoNVq6dOnD66urhw6dIgbN24watQoVCoVX375ZaWfr169euzevZvr169Tv359w/oFCxbg6elZ6eerLOK3k3BfK47FM3NbJABTn/NjQIB71Z4w/HeY/6w+2Nh6wrjtIthUszxNHssuLqPf+n68tvM1jt44ikKmoKdXT37v8zuLey2mq2dXEWwEAcjNzWXs2LFYWVnh4uLCzJkzSUhIwNzcnNzc3Co55/bt27lw4QJLly6ldevW9OrVi2nTpjFnzhyKi4vveUxcXBwymYxVq1bx9NNPo1aradeuHZcvX+b48eO0bdsWS0tLevXqRWpqapljnZ2d6dGjB4sXLzasO3ToEGlpafTp06dKvsbKIK7cCPe0NSKJD/84B8C/Ojci7KkqHB5fUwhb3odTt/7z+DwLg34B82qeVfwJFp8dz++XfueP6D/I0+QBYKWyYmDjgQxvNhw3SzcjVyg8aSRJokCjrfbzqlWKCvfYGTNmDOfOnWPPnj0kJyczaNAgIiIi6N69O5aWlvc9rrxtACNGjGDu3Ln33Hb48GH8/f1xcbk9OXFoaCgTJ07k/PnzBAQE3Pd9P/nkE2bPno2npydhYWEMGzYMKysr/vvf/2Jubs6QIUOYOnUqP/30U5njwsLCeP/99/noo48A/VWb4cNr9sCpItwIdzkSc5M3V5xGJ8GLbT14L7RJ1Z0sPUbfzTvpHCCDzh/AM++Lbt7VQJIkDt84zPKLy9l3fR8SEgBe1l4Mbzacfo36Ya4S444IxlGg0eI3dVu1n/fC56EVGrsrLS2NdevWsWzZMgIDAwEYOHAgS5YsYf78+eUeGx4eXu52a+v7j+CdlJRUJtgAhuWkpKRy3/fdd98lNDQUgLfeeouXXnqJnTt30rFjRwDGjRvHokWL7jruueee49VXX2Xfvn0EBgayatUqDhw4wIIFC8o9nzGJcCOUcT4xiwmLT1BcoqOHnwvTB7aourEnLmyAP1+Domwwd4Tn54nZvKtBviafjTEbWX5xOVeyrhjWP+3+NMObDSfELUTcdhKEB4iOjkaSJEJCQgzrgoKCWL16Nf369Sv3WB8f4wyX0LLl7dHcSwORv79/mXUpKSl3HadSqRgxYgQLFy4kJiYGX1/fMu9VE4lwIxjE38xn9ILj5BSVEORtz/9eCkBZFdMqaDWw4xM4Mke/7NEeBi8Ea3Hroyol5Caw4tIK1katNXTlNleaM8BnAC81fQkvGy/jFigId1CrFFz4PNQo560IU1P9VCJ3TvDo5OSEr68vjo7l9x58nNtSrq6uHDt2rMy65ORkw7by3Dnib+kfrf9cp9Pp7nlsWFgYwcHBREREEBYWVu55agIRbgQAUnOKGLngKGm5RTR1tWLeqLaYVfA/+UPJug6rx8L1W/85O7wB3T4BxaMPsy3cnyRJHE06yopLK9h9bbdhagQPKw+GNR1Gf5/+WJlU8ZhFgvAIZDJZ9Uzt8oi8vb2Ry+VERUXh5qb/w2zDhg3Ex8cjSVK5V7wf57ZUSEgI06dPJyUlBWdnZwB27NiBtbU1fn5+D/+FVFDz5s1p3rw5Z8+eZdiwYVV2nspSc39yhGqTXahh9IJjXL2Zj4e9miVhQdioqyBsRP8NayfoRxs2tYGBP0HTmtvavjbLKspiw5UNrIpcRVx2nGF9SL0QhjcbztP1nxa3ngThMdja2jJo0CCmT59OUFAQly9fZuvWrajVanbt2kW3bt3ue+zj3Jbq0aMHfn5+jBw5kq+//pqkpCQ+/vhjXnvtNcPVpKqya9cuNBoNtra2VXqeyiDCzROuUKNlwuITXLiRjaOlCb+FBeNsbVa5J9FpYc9/YN9M9KMNt4LBi8G+CntgPaHO3zzPyksr2RK7xTArt4XKgucaPsdLTV+ikW0jI1coCHXHnDlzGD9+PO7u7igUCmbPno2pqSnDhw9n+vTpjBs3rtLPqVAo2LhxIxMnTiQkJAQLCwtGjx7N559/Xunn+icLC4sqP0dlkUmSJBm7iOqUnZ2NjY0NWVlZ5V76exKUaHX8a9kptl9IxtJUyYqX29PC3aZyT5KbCmvHQexe/XLbMAidAapKDlBPsMKSQrbGbWVV5CrOpZ0zrG9s15ihTYbSp2EfLFS155eS8OQpLCwkNjYWb29vzMzE74YnWXk/Cw/z+S2u3DyhJEnioz8i2H4hGROFnHmj2lZ+sLl6SN++JjcJVBbQ97/QcnDlnuMJdjX7KqsiV7E+ej3ZxdmAfgLLZxs8y9CmQ2nt1FrMsiwIwhNJhJsn1Mxtkaw8cQ25DP73UgAhjRwq7811Ojg4G3Z9AZIWnJrqJ710qsLxcp4QGp2Gvdf2svryag4lHjKsd7NwY3CTwQz0GYiDuhL/LQVBEGohEW6eQL/uj+HHPfrxTb4c6E/PFuV3H3wouSnwxytwZZd+ueWL8Nx3YCJuizyO+Ox41kat5c/oP7lZeBMAGTKecn+KoU2H0tGtIwp5FfRuEwRBqIVEuHnCrDt1nS82XQTgvdAmDA2qxInPYvbCugmQmwxKNfSeqZ8bStwaeSRF2iL+vvo3a6PWcjzpuGG9g5kD/X36M9h3MPWt6pfzDoIgCE8mEW6eILsuJfPemrMAhHX05l+dK6nnjLYE9v4H9n0DSODUTD8on3Ozynn/J0xURhTrotbxV8xfZBVlAfqrNB3dO/JC4xd4xuMZVHIxLpAgCML9iHDzhDh5NZ1/LTuFVicxMMCdj/s0q5zGplkJsHY8xN9q/9FmNPT8D5iIOYkeRr4mn21x21gTtYazqWcN610tXBnoM5CBPgOpZ1nPiBUKgiDUHiLcPAEik3IYu/A4hRodnZs48fULLZHLKyHYRG6F9RP1g/KZWEHf2eD/wuO/7xNCkiTOpp1lffR6tsRuMczGrZQp6eTRiecbP08Htw6iLY0gCMJDEuGmjruWns+oBUfJLiyhjactPw5vg+px54sqKYa/P709N1S91vDCAnAQA8RVRFJeEhtjNvJn9J9lRg/2tPJkUONB9Pfpj6O6/LlpBEEQhPsT4aYOS8stYtSCYyRnF+HrYsmCMe0ef66W9BhYEwaJp/XL7f8F3T8FZdUO+13bFZQUsCt+FxuubOBw4mEk9GNnqpVqunt2Z4DPANq6thVTIgiCIFQCEW7qqJxCDWMWHiM2LQ93WzVLwoKxNTd58IHliVgLf/0birLBzBYG/ARNe1dGuXWSJEmcST3D+uj1bIvbRq4m17At0CWQ/o3608Orhxg9WBAEoZIZ/c/EOXPm4OXlhZmZGcHBwXdN5f5PmZmZvPbaa9SrVw9TU1N8fX3ZvHlzNVVbOxSVaHnlt5NEJGRjb2HCb+OCcLV5jCHNi3Lgj4n6KzZF2eDRHiYeFMHmPm7k3uCXs7/Qd31fRm4ZydqoteRqcnG3dOfVVq+yeeBmFvVcxMDGA0WwEQThocXHx9OnTx/Mzc1xdnbmvffeo6Sk5L77x8XFMW7cOLy9vVGr1TRq1IhPPvmE4uLiSq9t0aJFyGQymjW7u7fs6tWrkclkeHl5Vfp5/8moV25WrlzJO++8w9y5cwkODmb27NmEhoYSGRlpmMr9TsXFxTz77LM4OzuzZs0a3N3duXr1aq2YobS6lGh1vLH8NIeu3MTCRMGise1o6GT56G94/YS+N1RGLMjk8PT/QacPQCEu+t0puzibHXE72By7meNJx8vcdnq2wbMM8BlAoEuguO0kCMJj0Wq19OnTB1dXVw4dOsSNGzcYNWoUKpWKL7/88p7HXLp0CZ1Ox88//4yPjw8RERFMmDCBvLw8vvnmm0qv0cLCgpSUFA4fPkxISIhh/fz58/H0rMSx1cojGVFQUJD02muvGZa1Wq3k5uYmzZgx4577//TTT1LDhg2l4uLiCp+jsLBQysrKMjyuXbsmAVJWVtZj11/TaLU66e0Vp6UGkzZKjT/aLB2MSn2MNyuRpL1fS9KndpL0ibUkzWouSXEHK6/YOqCwpFDaHrddemvXW1LAkgCpxaIWhsfYrWOl9VHrpbziPGOXKQg1XkFBgXThwgWpoKDA2KU8lJycHGnMmDGSpaWl5OzsLH399dfS9evXJbVaLeXk5FTJOTdv3izJ5XIpKSnJsO6nn36SrK2tpaKiogq/z9dffy15e3uXuw8gzZ07V+rTp4+kVqulpk2bSocOHZKioqKkTp06Sebm5lJISIgUHR1tOGbhwoWSjY2N9Prrr0vjx483rL927ZpkamoqffDBB1KDBg3ue87yfhaysrIq/PlttD+/i4uLOXnyJJMnTzask8vldO/encOHD9/zmA0bNhASEsJrr73Gn3/+iZOTE8OGDWPSpEkoFPfuLjtjxgw+++yzKvkaahJJkvjsr/OsO52AQi5jzrA2dPB5xB43mddg3cu3x65p8Tz0mQVq20qrt7bS6rQcTz7OpphN/H317zLtaHxsfejTsA+9vHvhbuluxCoFoQ6QJNDkV/95VeYVHlV9zJgxnDt3jj179pCcnMygQYOIiIige/fuWFre/4p5edsARowYwdy5c++57fDhw/j7++Pi4mJYFxoaysSJEzl//jwBAQEVqj0rKwt7e/sH7jdt2jRmzZrFrFmzmDRpEsOGDaNhw4ZMnjwZT09PwsLCeP3119myZUuZ48LCwujcuTP//e9/MTc3Z9GiRfTs2bNM3VXJaOEmLS0NrVZ71xfq4uLCpUuX7nlMTEwMu3btYvjw4WzevJno6Gj+9a9/odFo+OSTT+55zOTJk3nnnXcMy9nZ2Xh4eFTeF1JDzNpxmcWHryKTwbeDW/Gs3yP+AEWshb/ehqIsMLGE3t9Aq6FP9BQKkiRxIf0Cm2I2sTV2K6kFqYZtrhau9PbuTW/v3jSxFxODCkKl0eTDl27Vf94PEys0F15aWhrr1q1j2bJlBAYGAjBw4ECWLFnC/Pnzyz02PDy83O3W1tb33ZaUlHTPz83SbRURHR3N999/X6FbUmPHjmXIkCEATJo0iZCQEKZMmUJoaCgAb731FmPHjr3ruICAABo2bMiaNWsYOXIkixYtYtasWcTExFSoxsdVqxpO6HQ6nJ2d+eWXX1AoFAQGBpKQkMDMmTPvG25MTU0xNa3b3ZR/2XeF73dFA/B5/xYMCHiEqwZFObD5fTizXL/s3haenwf2DSux0tolJjOGbVe3sTlmc5nxaKxNrOnh1YM+3n1o49JGtKMRhCdQdHQ0kiSVaVMSFBTE6tWr6devX7nH+vj4VHV595WQkEDPnj0ZPHgwEyZMeOD+LVu2NLwuDVH+/v5l1hUWFpKdnX1XKAsLC2PhwoV4enqSl5dH7969+eGHHyrpKymf0cKNo6MjCoWC5OTkMuuTk5Nxdb33LNX16tVDpVKVuQXVrFkzkpKSKC4uxsTkMbs610K/H4vny836K13v92zCyPYNHv5N7tloeBIonrz5i2KyYtget51tcduIzow2rDdVmNLZozN9vPvwlPtTqJ7A740gVCuVuf4qijHOWwGlfzTf+bnj5OSEr68vjo7lNwl4nNtSrq6ud/UqLv0cvd9nZ6nExES6dOlChw4d+OWXX8rdt5RKdft3XemUPfdap9Pp7jp2+PDhvP/++3z66aeMHDkSpbL6IofRwo2JiQmBgYHs3LmTAQMGAPpvzs6dO3n99dfveUzHjh1Zvnw5Op0OuVz/1/Lly5epV6/eExlsNpxJ5MM/zgEwsXMj/tX5If8a0Gpg79ew/1uQtGDjAYN+gQYdqqDamisuK45tcdvYdnUbURlRhvVKuZKQeiGEeoXSzbMbliaP0etMEISHI5NV6PaQsXh7eyOXy4mKisLNTX/7bMOGDcTHxyNJUrlz9z3ObamQkBCmT59OSkqKoVfxjh07sLa2xs/P777HJSQk0KVLFwIDA1m4cKHhM7Qq2dvb069fP1atWnXfsFZVjHpb6p133mH06NG0bduWoKAgZs+eTV5enuH+3ahRo3B3d2fGjBkATJw4kR9++IG33nqLN954g6ioKL788kvefPNNY34ZRrHrUjLvrAxHkmB4sCfvhz5ke4+Ui/DHK3DjjH65xQvQ59snptFwXFYc269uZ3vcdiIzIg3rlTIl7d3aE+oVShePLtiY2hixSkEQaipbW1sGDRrE9OnTCQoK4vLly2zduhW1Ws2uXbvo1q3bfY99nNtSPXr0wM/Pj5EjR/L111+TlJTExx9/zGuvvWa4mnTs2DFGjRrFzp07cXd3JyEhgc6dO9OgQQO++eYbUlPvaDf4gKs9j2vRokX8+OOPODg4VOl5/umRwk18fDxXr14lPz8fJycnmjdv/kjtWl588UVSU1OZOnUqSUlJtG7dmq1btxru68XHx5dJlx4eHmzbto23336bli1b4u7uzltvvcWkSZMe5cuotY7E3GTi0lOU6CT6t3ZjWv8WFZ/hW6eFIz/CzmmgLQK1nb4nVItBVVt0DRCTFcPOqzvZfnU7l9JvN1pXypQEuwUT2iCUrp5dRaARBKFC5syZw/jx43F3d0ehUDB79mxMTU0ZPnw406dPZ9y4cZV+ToVCwcaNG5k4cSIhISFYWFgwevRoPv/8c8M++fn5REZGotFoAP2VnejoaKKjo6lfv36Z99P3+K46arUatVpdpee4F5lUwa8sLi6On376iRUrVnD9+vUy3xATExOefvppXn75ZZ5//vlqudz1qLKzs7GxsSErK6vcS3811ZlrmQybd4S8Yi3dmznz04jAik+EmREH6/8FVw/qlxv3gH7fg1XVJndj0Uk6zqedZ2f8TnbG7yzTKFghU9C+nv4KjQg0gmBchYWFxMbG4u3tjZnZY4ymLtR65f0sPMznd4Wu3Lz55pssXryY0NBQvvjiC4KCgnBzc0OtVpOenk5ERAT79+9n6tSpfPbZZyxcuJB27do9+lcn3NPl5BxGLzxGXrGWkIYO/DCsgjN8SxKc/g22TobiXFBZQM8voc3oOtfFW6PTcCLpBLvid7Hr2i5S8lMM25RyJcGuwTzb4Fm6enbFzszOiJUKgiAIVaVC4cbCwoKYmJh73jNzdnama9eudO3alU8++YStW7dy7do1EW4q2ZXUXIbNO0pmvoZWHrbMG90WM9W9By4sIycZNrwBUdv0y54h+gkv7b2rtuBqVFBSwKGEQ+yM38ne63vJLs42bFMr1Tzt/jTdPLvxdP2nsTKxMmKlgiAIQnWoULgpbdBbET179nzkYoR7i0vLY9i8I6TlFuFXz5rFY9thaVqBf7rz62Hj21CQDgoT6DoFQl4DeQVCUQ2XVpDGgYQD7I7fzaHEQxRqCw3b7Ezt6OLZha4eXWnv1h5TRd0e50gQBEEoq1YN4vckup6Rz/Bfj5KcXYSviyVLxwdja/6Abu+5qbD5/+DCn/plV38Y+Au43L+bYE0nSRKX0i+x9/pe9l/fz7m0c4bJKQHcLNzo6tmVbp7dCHAOQFEHApwgCILwaCocbrp27Vqh/Xbt2vXIxQhl3cgq4KV5R0jILKChkwXLxrfH3qKcYCNJ+ukTNr+nv1ojU8DT78Az74Oy9o0DVFBSwNEbR9l7fS/7ru8r034GoJl9Mzp5dKKbZzea2DWpeI8xQRAEoU6rcLjZs2cPDRo0oE+fPmVGJxSqRkp2IcPmHeVaegENHMxZPr49Tlbl3F7JSYKN70DkJv2yiz8MmAP1WlVPwZXkRu4N9l3fx97rezmWdIwibZFhm1qppn299nSq34mn6z+Ns7mzESsVBEEQaqoKh5uvvvqKhQsXsnr1aoYPH05YWBgtWrSoytqeWGm5RQz79SixaXm426pZPqE9rjb36R4pSRC+HLZNhsIskKug0/vw1Nu1YvqEYm0xp1JOcTDhIAcTD5YZIRj0t5ueqf8MnTw60c61nWg/IwiCIDxQhcPNe++9x3vvvcfhw4dZsGABHTt2pEmTJoSFhTFs2LBaOWZMTZSRV8yIX48SnZKLq7UZv09oj7vtfQZAyroOf70F0X/rl90CoP8ccGlefQU/JEmSuJp9lYOJBzmYcJATyScoKCkwbJfL5LRyaqUPNPU74WPrI243CYIgCA+lwoP4/VN+fj6rV69mzpw5XLhwgcTExFoRcGryIH5ZBRqG/3qEiIRsHC1NWfVKexo63WM+I50OTi2G7VOgOAcUptBlMoS8AYqa10Y8tziXo0lHOZhwkEOJh0jITSiz3VHtSAe3DnR060iIW4gYf0YQnjBiED+hVLUO4ncvp06dYu/evVy8eJEWLVqIdjiPKadQw+gFx4hIyMbewoTlE4LvHWxSL+uv1sQf0i/XD9JfrXHyrd6Cy6HVablw8wJHbhzhQMIBzqaepUQqMWxXypUEOgfSwV0faHztfMXVGUEQBKHSPFS4SUxMZNGiRSxatIjs7GxGjBjB0aNHy52JVHiw0mATfi0TG7WKpeOC8XX5x2BzJUWwfxYcmAXaYlCZQ9ePIfhVo49bI0kSVzKvcDTpKEdvHOVE0glyNDll9mlg3cBwdaadazvMVeZGqlYQBKF2i4+PZ+LEiezevRtLS0tGjx7NjBkzUCof/JFeVFREcHAwZ86c4fTp07Ru3bpSa1u0aBFjx46ladOmXLx4scy21atXM2TIEBo0aEBcXFylnvefKhxuevfuze7du+nRowczZ86kT58+FfpGCuXLvhVsTsdnYm2mZOm4YPzc/nG5Le4gbPw3pF3WLzcOhT7fgK1ntddbKjE3kaM3jnLkxhGOJR0jrSCtzHYrlRXtXNvRwa0DHdw74GHlYaRKBUEQ6g6tVkufPn1wdXXl0KFD3Lhxg1GjRqFSqfjyyy8fePz777+Pm5sbZ86cqbIaLSwsSElJ4fDhw4SEhBjWz58/H0/P6vncqnA62bp1K/Xq1SM+Pp7PPvuMzz777J77nTp1qtKKq+uyCzWMmn/7is2y8cG0cL9jAseCDNjxib59DYCFM/T6CpoPrPY5oW4W3OR40nHD1ZlrOdfKbDdVmBLgHEBwvWDa12tPM/tmYiA9QRDqtNzcXN544w3WrFmDubk57777LsOGDaNx48akpKRgaXmPpgWPafv27Vy4cIG///4bFxcXWrduzbRp05g0aRKffvopJib3H9Nsy5YtbN++nbVr17Jly5YHnksmkzF37lz++usvdu3aRYMGDViwYAFOTk6MHz+e48eP06pVK3777TcaNWpkOE6pVDJs2DAWLFhgCDfXr19nz549vP322/z++++P/414gAqHm08++aQq63jiZBVoGLXgGGeuZWJrrr8VZQg2kgTn18GWDyDv1sB1gWOg+6egrp7Gtsl5yZxIPsHJ5JOcSD5BbFZsme0KmYIWji0MYaalU0vRTVsQhEojSVKZnpTVRa1UV7gN4JgxYzh37hx79uwhOTmZQYMGERERQffu3csNNg8KPSNGjGDu3Ln33Hb48GH8/f1xcXExrAsNDWXixImcP3+egICAex6XnJzMhAkTWL9+PebmFW8WMG3aNGbNmsWsWbOYNGkSw4YNo2HDhkyePBlPT0/CwsJ4/fXX7wpLYWFhdO7cmf/+97+Ym5uzaNEievbsWabuqiTCjRFk5WsYueAoZ69nYWuuv2LT3O1WsEmPhS3vQ9R2/bJjE+j7X2gQcv83fEySJHE99zonk0/qw0zSCa7nXr9rv8Z2jQl21YeZQJdALE0q/68SQRAE0I9QHrw8uNrPe3TY0Qq1CUxLS2PdunUsW7aMwMBAAAYOHMiSJUuYP39+uceGh4eXu728nkBJSUl3BYTS5aSkpHseI0kSY8aM4dVXX6Vt27YP1d5l7NixDBkyBIBJkyYREhLClClTCA0NBeCtt95i7Nixdx0XEBBAw4YNWbNmDSNHjmTRokXMmjWLmJiYCp/7cYhGM9UsK1/DiPlHOZeQhZ25imXj2+vb2GgK4OB/9Y2GtUX6iS6ffhee+jcoK/eKiCRJxGbHGoLMyeSTJOcnl9lHLpPT1L4pgS6BtHVpSxvnNtia2VZqHYIgCLVVdHQ0kiSVaVMSFBTE6tWr6devX7nH+vj4VHV5ZXz//ffk5OQwefLkhz62ZcuWhtelIcrf37/MusLCQrKzs+8KZWFhYSxcuBBPT0/y8vLo3bs3P/zwwyN+FQ/nocPNzZs3mTp1Krt37yYlJQWdTldme3p6eqUVV9dk5hczYv5RQ3fvZeODaVbPGi5v01+tyYjT79iwM/T+BhwbV8p5i7RFXLh5gfCUcP0jNZz0wrL/TkqZkuaOzWnr0pZAl0ACnAPElRlBEIxGrVRzdNhRo5y3IkxN9X903tnGxcnJCV9fXxwdHcs99nFuS7m6unLs2LEy65KTkw3b7mXXrl0cPnzYUHOptm3bMnz4cBYvXnzfWu4c5qX0dt291v0zCwAMHz6c999/n08//ZSRI0dWayekhz7TyJEjiY6OZty4cbi4uIjxSSooI08fbM4n3h7HpqlpBvz+6u35oKzcoOeX4DfgsRoMp+anEp56O8hcvHkRjU5TZh8TuQktnVrS1lUfZlo5tarwf2pBEISqJpPJavSQEd7e3sjlcqKionBzcwNgw4YNxMfHI0lSuZ+Nj3NbKiQkhOnTp5OSkoKzs35+vR07dmBtbX3fYVn+97//8cUXXxiWExMTCQ0NZeXKlQQHV92tP3t7e/r168eqVavuG9aqykOHm/3793PgwAFatapdEzIaU0pOISN/PUZkcg4OFib8HhaA7+WfYd+3UFIAciW0/xd0mgSmD3e1pERXQnRmtCHIhKeE3zUCMIC9mT2tnVrT2rk1rZxa0cKxBSaK2jdTuCAIQk1ga2vLoEGDmD59OkFBQVy+fJmtW7eiVqvZtWsX3bp1u++xj3NbqkePHvj5+TFy5Ei+/vprkpKS+Pjjj3nttdcMV2aOHTvGqFGj2LlzJ+7u7nd1vy69ctSoUSPq16//yLVUxKJFi/jxxx9xcHCo0vP800OHm6ZNm1JQUP0t2GurxMwCht+aBNPZypQ/ehTgvvZZSL/VqMrraeg9E5ybPfC9JEkiOT+ZiLQIzqWdIyItgoi0CPJL8svsJ0NGY7vGhjDT2qk19a3qi6tsgiAIlWjOnDmMHz8ed3d3FAoFs2fPxtTUlOHDhzN9+nTGjRtX6edUKBRs3LiRiRMnEhISgoWFBaNHj+bzzz837JOfn09kZCQajaacd6oearUatbr67wo89NxSx48f54MPPmDq1Kn3nHahps3X9E/VObfU1Zt5DJt3lITMAtpbp7PAbT3mcbcmubR0gR7Twf+F+96Cyi7O5nza+TJhJrUg9a79LFWWtHRqSWun1rRybkVLx5aivYwgCLWGmFtKKGW0uaVsbW3Jzs6ma9euZdaX3mPUarUP+5Z1UlRyDsN/PUpBTgZfWW1kSMkmZHEa/S2ooJeh82Qwu/2PU6QtIjI9sswVmbjsuLveVyFT0NiuMc0dmuPv6I+/kz+NbBqJAfMEQRAE4ZaHDjfDhw9HpVKxfPly0aD4Ps4nZjH618M8W7SdSeo12Gqy9Bt8noXQL8m3rc/ljMtcjLvIpfRLXLx5kajMKEp0JXe9V33L+vg7+tPCsQX+Tv40tW8qGv4KgiAIQjkeOtxERERw+vRpmjRpUhX11Hqn4jP4fsEilugW4qe6Sg4yjrv4cNGvFxflWi4efJ/Y7Fh00t3d5uxM7fQh5laYaeHYAjuz6hmRWBAEQRDqiocON23btuXatWsi3PyDpNWy9us53Cj6g4420cw3UXHRzJ1rSgVQDHF/ltnfwcwBPwc/mjk0o5l9M5o5NMPNwk1cCRMEQRCEx/TQ4eaNN97grbfe4r333sPf3/+uBsV3jmb4JFn8n3e5eTYGhaw+VjnONDQ/h9xVhrmrRJFXPRrW0wcZPwc/mtk3w8ncydglC4IgCEKd9NDh5sUXXwT0wyqXkslkT3yDYgeTRqRzDa2UQYK1PY66CQw5tBHLvESQJWDipcLMT4WZnwlmzU3QNjNBYWPz4DcWBEEQBOGhPHS4iY2NffBOT6Dn3n2N43sbcnThUooL00mTn+BIm9fwzIrG6+JqiI2lODaW7E2bDMeo6tfHzM8PM79m+md/f5R2oo2NIAiCIDyOhx7nprar6nFuslKSWfnZh+SkJYPMEhOrwZhZONDSD7xKItFcukDhhQtort896zboA4+6pT9m/i31z35+yI0wAJIgCEJ1EePcCKUqa5ybCoWbI0eO0L59+woVlp+fT2xsLM2bN6/Q/tWtOgbxy0lPY820j0lPvI5caYHSfBByhRO2LuZ0fN6HBv4O6LKzKbx4kcILFym8cIHCiAiK7zUNvUKBaePGqP39MWvpj7plS0wbNUJWjROQCYIgVCURboRS1RpuGjduTMOGDRk/fjy9e/fGwsLirn0uXLjA0qVLWbhwIV999RWjRo16yC+pelTXCMX52VmsmT6F1LgYVGYWmNk8j6ZYP1OsZ3N7nh7ii61L2UnhtNnZFEZEUHD2HAXnzlF49iwlqXePSCxTqzHz80Pt74+6VUvUAQGo7jMbrCAIQk0nwo1QqlrDjUaj4aeffmLOnDnExMTg6+uLm5sbZmZmZGRkcOnSJXJzcxk4cCAffvgh/v7+j/fVVaHqnH6hMDeXdTM+4UZ0JCozNb4dJhB7VoVOKyFXymjTowFtejZAZXL/0YU1yckUnD1LYWngOXcOXV7eXfsp3eph3joAdYD+YdbEF9k/erIJgiDURCLcPJz4+HgmTpzI7t27sbS0ZPTo0cyYMQPlA67ob9q0ic8//5yzZ89iZmZGp06dWL9+faXWtmfPHrp06YKtrS03btwo8+95/PhxgoKCAP2sBvdSreHmTidOnODAgQNcvXqVgoICHB0dCQgIoEuXLtjb2z/MWxlFdYYbgOKCfP74+nOuX4hAaWpK9/HvcSXcjGsX0gGwcjDj6Rd98W7pWKH3k3Q6imNjKTh7jsJzZykIP0PhpUugKzsooEyt1l/ZCQhAHdAadatWorGyIAg1kgg3FafVamndujWurq7MnDmTGzduMGrUKCZMmMCXX3553+PWrl1r2Kdr166UlJQQERHBkCFDKrW+0nDj4eHBV199xUsvvWTY9uqrr7Jlyxbi4+NrXrip7ao73ABoigrZMGsGceEnUapM6PfuR+h0HhxYHUVuRhEAXv4OPP2iL9aOD994WJeXR8G5CArCT5N/+jQFp8PRZWfftZ9Jw4aoA1pjHhCAOjAQEy8vMWigIAhGV1vDTW5uLm+88QZr1qzB3Nycd999l2HDhtG4cWNSUlKwtKz8CYy3bNnCc889R2JiIi4uLgDMnTuXSZMmkZqaiomJyV3HlJSU4OXlxWefffZQM5V7eXkxfvx4Ll++zLp163BwcOD7778nJCSE8ePHs3PnTho2bMiCBQto27YtcDvcfPzxxxw5coQdO3YAUFBQQL169XjzzTeZNm1alYcbeYW/SuGRqUzN6P/uxzRqG0yJppg/v/kChfwawz5tT5tQT+RyGXHnbrL8s6Mc3xRLiebhxgqSW1hg0T4Yx1dfxfPnn/E9cpiGmzbiOu1zbAYNwsTbG4DimBiy1q7jxsdTiOnVm6hnnuH622+TvmwZhZGXkXR3TwkhCIJgDJIkocvPr/bHw/y9P2bMGA4dOsSePXtYuHAhU6ZM4cMPP6R79+7lBhtLS8tyH6+++up9jz18+DD+/v6GYAMQGhpKdnY258+fv+cxp06dIiEhAblcTkBAAPXq1aNXr15EREQ88Gv87rvv6NixI6dPn6ZPnz6MHDmSUaNGMWLECE6dOkWjRo0YNWrUXd+3kSNHsn//fuLj4wH9lSMvLy/atGnzwHNWBtHlppooVSr6vv0Bf333FVdOHGH9N1/Q//8+ImRgW5q0r8e+FZdJiMzg2F+xXDqSxDNDfWnQ3OGRziWTyzFt1AjTRo2wGzwYgJKMDArCwyk4HU7+qZMUnj2HNjWNnC1bydmyFQCFjQ3qtm0xb9cW87btMGvaRPTKEgTBKKSCAiLbBFb7eZucOonM3PyB+6WlpbFu3TqWLVtGYKC+zoEDB7JkyRLmz59f7rHh4eHlbi/vqkRSUlKZYAMYlpOSku55TExMDACffvops2bNwsvLi2+//ZbOnTtz+fLlcpuU9O7dm1deeQWAqVOn8tNPP9GuXTsG3/psmTRpEiEhISQnJ+N6R8cWZ2dnevXqxaJFi5g6dSoLFiwoM/hvVROfXNVIoVTR9+1JbJz9NdHHD/PnN1/Q792PaBjQjv7/bk30yRQOro4iO7WAjd+foXFbZ54a4ou59d2XGR+W0s4Oqy5dsOrSBQBdURGFZ8+Sf+IE+cdPkH/6NNqsLHJ37iR3505Af0VIHdgG87btMG/bFnWL5sjucclTEAThSRMdHY0kSYSEhBjWBQUFsXr1avr161fusT4+PlVdXhm6W1flP/roI55//nkAFi5cSP369Vm9erUhvNzLnVMqlYaoOzsNla5LSUkpE25AP5PBW2+9xYgRIzh8+DCrV69m//79lfNFPYAIN9VMoVTx3L8nsem/XxN17BAbvplO33c+pFFgEI3butCghQPH/orl7K5rRJ1IIf5COh2e96FZh3qV2j5GbmqKebt2mLdrBxNB0mgovHBBH3aOHSf/1Cl0OTnk7dtP3j79D6NMrcY8MBCLkPaYt2+PWbNmyOTizqYgCJVPplbT5NRJo5y3IkxNTQHKtHFxcnLC19cXR8fyO4g8qC3OiBEjmDt37j23ubq6cuzYsTLrkpOTDdvupV69egD4+fmVqb9hw4aG20b3c+f8kaWfQfdap7tHs4ZevXrx8ssvM27cOPr27YuDw6PdjXgUItwYgUKppM9b77Ppf18TdfQQG779kn7/N5lGgcGYmCl5anBjfINc2L30EmnXctn92yUijyTReXgT7FzvHmOoMshUKtStWqFu1QqHceOQtFqKIiNvX9k5cQJtRgZ5Bw6Qd+CA/uuwscE8OFgfdoLbY+ItGigLglA5ZDJZhW4PGYu3tzdyuZyoqCjc3NwA2LBhg6EnUHm/Cx/ntlRISAjTp08nJSUFZ2dnAHbs2IG1tXWZ8HKnwMBATE1NiYyM5KmnngL0Q7zExcXRoEGDcmt5HEqlklGjRvH111+zZcuWKjvPPc9dkZ3eeeedh37jjz/+uFZ0DTcWhVJJnzffZ/P333D5yAE2fDuDvm9/gE87/UjQzg2sGfxBW87uvs7RDTEkRmWy4otjBPb0IjC0AQpV1V4xkSkUt+a98sN+1CgknY6iqCjyDh8m/8hR8o8fR5uVRc727eRs3w6A0tUVi+BgzEPaYxESguof94UFQRDqCltbWwYNGsT06dMJCgri8uXLbN26FbVaza5du+jWrdt9j32c21I9evTAz8+PkSNH8vXXX5OUlMTHH3/Ma6+9ZriadOzYMUaNGsXOnTtxd3fH2tqaV199lU8++QQPDw8aNGjAzJkzAQxtZ6rKtGnTeO+996r1qg1UMNzMnj2bkJCQe3Yxu5cDBw7w+uuvi3DzAPqA8x4ymYzIw/v567v/lAk4coWc1t09adjaib2/Xyb+/E2Ob4wl+kQynYc3xa2xbbXVKpPLMWvSBLMmTXAYMwZJo6EgIoL8I0fIO3yEgtOnKUlKIuvPP8n6808ATLy9sejYEYunOmIRFIS8Bv8VJgiC8LDmzJnD+PHjcXd3R6FQMHv2bExNTRk+fDjTp09/qG7XFaVQKNi4cSMTJ04kJCQECwsLRo8ezeeff27YJz8/n8jISDQajWHdzJkzUSqVjBw5koKCAoKDg9m1axd2VTz+mYmJyQNv01WFCo1zI5fLSUpKMlwCexArKyvOnDlDw4YNH7vAymaMcW4eRKfVsvmHb4k8tA+5Qkn/9/SNjO8kSRLRJ1PYvyqKguxiAPyedqPDIB9M1ca/u6grKKDg9GnyDh8h78gRCs+fLzOwoEylQt02EMunnsLiqacx9W0sbmEJggDU3nFuhMpXrYP4LV68mKFDhxoueT3I8uXL6d+//z3noDK2mhhuQB9wNv1vJpePHEChUjHw/U9o0LL1XfsV5mk4/McVLhxIBMDSzpTOw5vSoEX1XvJ7EG12NnlHj5J34CB5+/ejSUwss13p7Hz7qk6HDmL0ZEF4golwI5QSIxQ/opoabgC0JSX89d1/uHLiCEoTUwZN/hQPv3vP05VwOYNdv10iO7UAgKbtXek4uDFmFjVvPilJkiiOjSPvwAFyD+wn/9hxpMLC2zvIZJj5+2P5VEcsnn4adcuWyBT3n29LEIS6RYQboZQIN4+oJocbgBKNhg3fTif29AlUpmY8/9E03Js0u+e+mmItR/+M4cyuayCBubUJnYY1oWFrp2qu+uHoioooOHmS3P36nldFUVFltivs7LB85hksu3TGomNHFFZWxilUEIRqIcKNUMpo4UYul5fbVkKrfbipA6pbTQ83ACXFxfzx9efEnwvHRG3O4I+/wNXH9777J8VksWvJRTKS8gFo3M6Fp19sjNqydgy4p0lKIu/gQX3YOXgQXU7O7Y1KJeZt22LVpTOWnTtjUoXdFgVBMA4RboRSRgs3f97qCVNKo9Fw+vRpFi9e/NCTchlDbQg3oJ9sc92MT7l+MQIzC0sGT/0SZ6/7N9Au0Wg5vjGW09vjkSRQW6l4ZmgTfAIr1gi8ppA0GvJPnSZ3zx5y9+yhODa2zHYTb28sO+uDjnmbAGSqmncbThCEhyPCjVCqxt2WWr58OStXrrwr/NQ0tSXcABQX5LPmy6ncuHwJtZU1Qz6ZgaNH+VcukuOy2bXkIumJeQA0auNEp5eaoLaqHVdx/qk4Lo7cvXvJ2bOH/OMnoKTEsE1uZYXl009h2a0blp06oaiCGXgFQah6ItwIpWpcuImJiaFly5bk5uZWxttVmdoUbgCK8vNYPe1jkmOiMLex5cVP/4O9W/1yj9FqdJzYGsepLVfR6STUViq6jGiKd6ua3RbnQbQ5OeQdPKS/qrN3L9qMjNsbVSos2rfHqnt3rLp1RWmEcRUEQXg0ItwIpWpUuCkoKGDy5Mls2bKFyMjIx327KlXbwg1AQW4Oq6d9RGpcDJZ29rz46VfYutZ74HGp8Tn8veiC4SpOsw71eGpwY0xqwLg4j0vSaik4e5bcXbvJ+fvvsrevZDLUAQH6oNO9GyaensYrVBCEBxLhRihltHBjZ2dXpkGxJEnk5ORgbm7O0qVLHzgbqrHVxnADkJ+dxarPJnPzejzWTi4M/ewrrBwefHWiRKPl6IZYwv+OBwms7M3oNqYZ7r51a1yZoitXyPl7Jzl//03huXNltpn6+uqDzrPdMW3aVAweKAg1jAg3QimjhZvFixeXWZbL5Tg5OREcHFzlwzhXhtoabgDyMjNY+ekkMm4kYu9Wnxc/+wpza5sKHZsYlcHOxRfJTisEGbTq5kH7/g1RqureeDKaGzfI2bmLnJ1/k3/sONzRg0/l7o5Vjx5Y9+qJmb+/CDqCUAOIcCOUqlG3pWqT2hxuALLTUlgxdRI5N1Nx9m7EkKlfYmpesZGgiwtLOLgm2jC6sZ2rOd3H+uHcoPZ9HypKm5lJzp495Pz9N3kHDpYZPFDl5oZVz576oNOihQg6gmAkItw8nPj4eCZOnMju3buxtLRk9OjRzJgxA6Xy/k0OLl++zHvvvcfBgwcpLi6mZcuWTJs2jS5dulRqbXv27KFLly7Y2tpy48aNMv+ex48fJygoCNDf9bmXygo3FZpa+uzZs+jumCfoQc6fP0/JHb1ahMpj7ejMCx9PQ21tQ0rsFf746jM0RYUPPhAwMVPSZURT+rzWEnNrEzKS8ln71UmOb4pFq634v29torC1xXbAADx++AHfw4dw//5/WPfpg8zcHE1iIukLFhA3eAhXuj9L8syZFJyLuO9/OkEQBGPTarX06dOH4uJiDh06xOLFi1m0aBFTp04t97jnnnuOkpISdu3axcmTJ2nVqhXPPfccSUlJVVKnlZUVf/zxR5l18+fPx7Oa2kBWKNwEBARw8+bNCr9pSEgI8fHxj1yUUD57t/q88NE0TM0tSLh0gQ2zZqAt0Tz4wFu8/B15aWowjdo4o9NJHPsrlnVfnyQzOb8KqzY+uVqN9bPP4v7tN/geOoj7//6Lde9e+qCTkED6/AXEDR7MlWd7kPLNNxREnBdBRxCMRJIkNEXaan88zP/53Nxcxo4di5WVFS4uLsycOZOEhATMzc2rrOfw9u3buXDhAkuXLqV169b06tWLadOmMWfOHIqLi+95TFpaGlFRUXzwwQe0bNmSxo0b85///If8/HwiIiLuey4vLy+++OILRo0ahaWlJQ0aNGDDhg2kpqbSv39/LC0tadmyJSdOnLjr2NGjR7NgwQLDckFBAStWrGD06NGP/02ogAp1m5EkiSlTpmBubl6hN73fN1ioPM5eDRk46RPWfDmFuPCTbP7+W/q89R5yecXa0JhZqgid0JyoE47s+/0yKVdzWPnlcZ4e0phmHerV+Vs0cjMzrHv0wLpHD3QFBeTu20/21i3k7tmL5vp1bv46n5u/zkfl4YF1z55YP/ccZk3uP0q0IAiVq6RYxy9v7a328778306oTCv2e3TMmDGcO3eOPXv2kJyczKBBg4iIiKB79+5YljPuVnnbAEaMGMHcuXPvue3w4cP4+/vj4uJiWBcaGsrEiRM5f/48AQEBdx3j4OBAkyZNWLJkCW3atMHU1JSff/4ZZ2dnAgMDy63lu+++48svv2TKlCl89913jBw5kg4dOhAWFsbMmTOZNGkSo0aN4vz582U+N0aOHMnMmTOJj4/H09OTtWvX4uXlRZs2bco9X2WpULh55plnHqqLd0hICGq1usL7z5kzh5kzZ5KUlESrVq34/vvvDfflyrNixQpeeukl+vfvz/r16yt8vrrCvakf/f/vI9Z//TmXjxzARK2mxytvVjiYyGQyfNu54uZjy9+LLpAQmcnu3y4Rf/4mnYc3rZGTcFYFuVqNdWgPrENvBZ29+8jeupXcvXvRXLvGzXnzuDlvHqa+vlj3fQ6bPn1QubkZu2xBEIwoLS2NdevWsWzZMkNAGDhwIEuWLGH+/PnlHhseHl7u9vLakyQlJZUJNoBh+X63mGQyGX///TcDBgzAysoKuVyOs7MzW7dufWBHoN69e/PKK68AMHXqVH766SfatWvH4MGDAZg0aRIhISEkJyfj6upqOM7Z2ZlevXoZbpktWLCAsLCwcs9VmSoUbvbs2VNlBaxcuZJ33nmHuXPnEhwczOzZswkNDSUyMhJn5/tPHRAXF8e7777L008/XWW11QZerdrQ5833+eu7/xCxewem5uZ0Gjn+oa68WNqZ0e+tAMJ3xHP0zxiunEolOTab7mP8cG9S83vAVSa5Wo11z1Cse4aiy88nd+9esjZtIm/vPoouXyb121mkfjsLddtAbJ7ri3XPUBS2tsYuWxDqHKWJnJf/28ko562I6OhoJEkiJCTEsC4oKIjVq1c/cEgUHx+fx6rxYUmSxGuvvYazszP79+9HrVbz66+/0rdvX44fP069evcfN61ly5aG16Uhyt/f/651KSkpZcINQFhYGG+99RYjRozg8OHDrF69mv3791fml3ZfFftXrEKzZs1iwoQJjB07Fj8/P+bOnYu5uXmZe3X/pNVqGT58OJ999hkNG95/vqUnRePgDoROfAuAk5v+5MjaFQ/9HnK5jDahDXh+UiA2zmpyM4pYP/s0h/+4UmcbGz+I3Nwc61698PjhBxof2I/r559hHhQEMhkFJ06S9OmnXH76Ga796zWyN29GV1Bg7JIFoc6QyWSoTBXV/qjoH4ampqYAmJjcntrGyckJX19fHB8wQrqlpWW5j1dfffW+x7q6upKcnFxmXenyP8NFqV27drFx40ZWrFhBx44dadOmDT/++CNqtfqu4V3+SXXH/H2l35t7rbtXp6NevXpRUFDAuHHj6Nu3Lw4ODuWeqzIZdaja4uJiTp48yeTJkw3r5HI53bt35/Dhw/c97vPPP8fZ2Zlx48Y9MAUWFRVRVFRkWM7Ozn78wmug5p26UZSfx+5Fv3Bo9TJMzc1p07v/Q7+PcwNrhnzYjoOro7hw8Aantl3l+qV0ng1rjq1Lxdpc1UUKGxvshgzBbsgQNDdukL15M1l/baTo0iVyd+0id9cu5ObmWD37LNbPPYdFhxBkiro3hpAgCHre3t7I5XKioqJwu3WbesOGDcTHxyNJUrkh6XFuS4WEhDB9+nRSUlIMdzd27NiBtbU1fn5+9zwmP1/fWUQuL3s9Qy6XP1RP6IelVCoZNWoUX3/9NVu2bKmy89yLUa/cpKWlodVq73n/8H73Dg8cOMD8+fOZN29ehc4xY8YMbGxsDA8PD4/HrrumatOrHx2GDAdg9+J5ROze8UjvY2KmpMvIZvR8uQWm5kpDY+MLBxNF7yFAVa8eDuPG0XD9HzT8awMOr7yCyt0dXX4+WX/+ybUJE4ju0pWUb7+l6MoVY5crCEIVsLW1ZdCgQUyfPp2CggLOnDnD1q1bUavV7Nq1q9xjfXx8yn2U1ySjR48e+Pn5MXLkSM6cOcO2bdv4+OOPee211wxXk44dO0bTpk1JSEgA9IHIzs6O0aNHc+bMGcOYN7GxsfTp06fyvin3MG3aNFJTUwkNDa3S8/yT0W9LPYycnBxGjhzJvHnzHnjZr9TkyZPJysoyPK5du1bFVRpX+0FDCewzAIDtP39P9PEjj/xejdo4M3RKEO5NbCkp0rL7t0tsmxdBYV7Fu53XdaaNG+P89r9p9PcOGixfht2wl1DY2FCSksLNeb8S0+c5Yoe8SMbvv6PNyjJ2uYIgVKI5c+ZgZmaGu7s73bt3Z/bs2cyePZvhw4c/sFHxo1IoFGzcuBGFQkFISAgjRoxg1KhRfP7554Z98vPziYyMRKPR/652dHRk69at5Obm0rVrV9q2bcuBAwf4888/adWqVZXUWcrExARHR8dq74Fr1BGKi4uLMTc3Z82aNQwYMMCwfvTo0WRmZvLnn3+W2T88PJyAgAAUd1zuL72kJpfLiYyMpFGjRuWes7aPUFwRkiSxbe5/Ob/nbxQqFS98OI36fi0e+f10OsnQ2Fink7CyN6PH+Oa4NqzY1A9PGl1xMbl79pD1x3py9+0zTP8gU6mw7NYN24EDsOjYEVk5o4kKwpNEjFAslKrWEYqriomJCYGBgezcudOwTqfTsXPnzjIt0Es1bdqUc+fOER4ebnj069ePLl26EB4eXqdvOT0MmUxGj5ffoFHbYLQaDX98/TkpcTGP/H53Nja2djQjJ72QP745xaltV5F04jbVP8lNTLDu0QOPn36k8d49OH8wCVNfXySNhpytW7n2yqtEdelC8syZFEVFGbtcQRCEOsfoc0utXLmS0aNH8/PPPxMUFMTs2bNZtWoVly5dwsXFhVGjRuHu7s6MGTPuefyYMWPIzMys8Dg3T8KVm1Ka4iLWTp9KwqXzmNvY8tLnM7F1vX+Xv4ooKihhz7JLRJ9IAcCzuT3dx/ihtjJ5wJFPNkmSKLp4kcz168n+ayPajAzDNrMWLbAZNBCbvn1RWFkZsUpBMA5x5UYoVSeu3AC8+OKLfPPNN0ydOpXWrVsTHh7O1q1bDY2M4+PjuXHjhpGrrJ1UJqYMeH8KTg28yc/KZM2XU8jLzHjwgeUwVSvpMa45nYc3QaGSE38+nRVfHCMh8vHet66TyWSY+fnh+uGHNN67h/o/fI9lt26gVFIYEUHy59OIevoZEj+YTP6pU6LhtiAIwmMw+pWb6vYkXbkplZeZwe9T3yMrOQmnBt4M+WQGZhblD/9dETcTctk2L4KMpHyQQdveXrTr441cXrenbqhMJTdvkvXXX2SuWUNx9O2eVSaNGmH7wgvYDOiP8gEjiApCbSeu3AilKuvKjQg3T4jMpBv8PvU98rMyqd+sBYM+/AyVieljv6+mSMv+lZe5eEh/dc2tsS3PhjXH0u7x3/tJIkkSBeHhZK5eQ/aWLUi3BgSUqVRYPdsd28GDMQ8ORiY3+sVWQah0ItwIpUS4eURPargBSImLYeWnH1BckE+jtu3p985k5JU00NzlY0nsWRaJpkiLmaWKbqOb4eVfse76Qlna3FyyN24ic/VqCs+fN6xXeXjor+YMHICqnHEwBKG2EeFGKCXCzSN6ksMNwLUL51j75VS0Gg0tujz7UBNtPkhmcj7bfo0g7VouAK27e9B+QCMUSnG14VEVXrhAxurVZP+1EV2u/vuKQoFll87YDX1JPxKyuJoj1HIi3AilRLh5RE96uAGIOn6Yv76dgSTpaNf/BZ4ZNqbS3lur0XFwXTTndl8HwLWhNT3Gt8DKXvzCehy6ggKyt24jc/VqCk6dMqxXeXpiN3QoNgMHiLY5Qq0lwo1QSoSbRyTCjd653dvZPvd/AHQaEUbbvoMq9f1jTqeyc8lFigtKMLNQ0T3MjwbNq2/StLqsKDqajFWryPpjPbqcHABkJiZY9+6N3UtDMWvZstpHAxWExyHCjVBKhJtHJMLNbcf+XMP+5YsA6Pmvt2neqVulvn9WagHb5kWQGp+j703Vy4t2z4neVJVFl59P1qZNZPz+O0UXLhrWm/n5YfvSUGz69EFu/uROdirUHiLcCKXqzDg3gvG06/e8YR6qbXP/y5WTxyr1/W2c1Ax6rw3Nn3EHCU5sjmPDf8PJzy6u1PM8qeTm5tgNHoz32rV4rVyBTf/+yExMKLxwgaQpU4nq1Jmk6V9SFPPoo1MLglDzvPnmmwQGBmJqakrr1q2r7DyLFi1CJpPRrFmzu7atXr0amUyGl5dXlZ3/cYhw8wSTyWR0GhGG3zNdkXQ6Nn73H65fOv/gAx+CUqWg87AmPBvmh9JUQUJkBiunHyMxSgz6V1lkMhnqVq1w++o/+Ozdg/N776Hy9ESXk0PGb78R07sPV0ePIXv7dqSSEmOXKwhCJQgLC+PFF1+s8vNYWFiQkpLC4cOHy6yfP38+np6eVX7+RyXCzRNOJpfT45U3adimHSWaYtZ//Tmp8XGVfh7fIFcGf9AWu3oW5GcVs/67cDE3VRVQ2tnhMC6MRlu34DFvHpZdu4JcTv7RoyS8+RZXeoRyc/58MUO5UONJkoSmsLDaHw/TUiM3N5exY8diZWWFi4sLM2fOJCEhAXNzc3JLezdWgf/973+89tprNGzYsMLHyGQyfv75Z5577jnMzc1p1qwZhw8fJjo6ms6dO2NhYUGHDh24cuVKmeOUSiXDhg1jwYIFhnXXr19nz549DBs2rNK+psom2twIAGiKClnzxRQSL1/E0s6eoZ/PxMbZpQrOo2XP8ktcPpoMgJe/A93G+GFmoar0cwl6msREMlasJHPVKrSZmQDI1Gps+vXDfsRwTBs3Nm6BwhPvXu0sNIWF/G/0C9Vey5uL16CqYLufF154gXPnzrF8+XKSk5MZNGgQL774IhkZGWzYsOG+x1lalj9C/IgRI5g7d+4Dz//pp5+yfv16wsPDH7ivTCbD3d2dWbNm0bp1ayZNmkR4eDgNGzbk/fffx9PTk7CwMGxtbdmyZQugvy3173//m927d9O5c2du3LiBubk5X3zxBSdOnKBz587Mnj2buLi4B56/oiqrzY2y0ioSajWVqRkDJ33Cik/e5+b1eNZ+OZWhn3+NubVNJZ9HQfcxfrj52LJ/ZRRx526yavpxQl9ugYuXCJtVQeXmhvM7b+P4r4lkb9pE+pLfKIqMJHPlSjJXrsSiQwh2I0Zi2bmTGDNHECooLS2NdevWsWzZMgIDAwEYOHAgS5YsYf78+eUe+6AwUlV/eI8dO5YhQ4YAMGnSJEJCQpgyZQqhoaEAvPXWW4wdO/au4wICAmjYsCFr1qxh5MiRLFq0iFmzZhFTg9vziXAjGJhZWvL8R5/z+5T3yLiRwLoZnzJk6nRM1JXb40Ymk9H8aXecG1iz9ZdzZKcVsu6bkzz1QmNadHIX3ZiriNzMDNvnn8dm0CDyjx8n47el5OzcSd6hw+QdOozK0xP74cOwGTRIzE4uGJ3S1JQ3F68xynkrIjo6GkmSCAkJMawLCgpi9erV9OvXr9xjfXx8HqvGR9WyZUvD69LJqf39/cusKywsJDs7+66AFRYWxsKFC/H09CQvL4/evXvzww8/VE/hj0D8mSaUYWXvyAsfTUNtZU1yTBR/fvslJRpNlZzLydOKIR+2w7uVI7oSiX0rLrNj/nmKC0Wj16okk8mwCAqi/vf/o9H27diPC0NubY0mPp7kGf8hulNnkqZ9QVFsrLFLFZ5gMpkMlZlZtT8q+seV6a0QZGJiYljn5OSEr68vjo7lTz1jaWlZ7uPVV1999G9cOVSq27f/S7/Oe63T6XR3HTt8+HCOHDnCp59+ysiRI1Eqa/a1kZpdnWAU9m71GfTBp6z6/EPiz4Wzdc4ser/5LnJ55cxDdSdTcxW9XvXnzM5rHF53hagTKaQl5NHrlRbYuVpU+vmEskzqu+Py3ns4vfYaWRv+In3pbxRHXyFj2TIyli3D4pmncRgzBvOQEHFFTRDu4O3tjVwuJyoqCjc3NwA2bNhAfHw8kiSV+//FWLelHoe9vT39+vVj1apVFWoPZGziyo1wT64+vvR79yPkCiWRh/eze9EvD9WL4GHIZDJad/ek/zsBmNuYkHEjj9X/OcGVUylVcj7hbnJzc+yGvkjDv/7Cc8F8LLt0AZmMvH37iQ8bR+zAQWSuX49ULMYoEgQAW1tbBg0axPTp0ykoKODMmTNs3boVtVrNrl27yj3Wx8en3IfzAybGjY6OJjw8nKSkJAoKCggPDyc8PJziKv7/uWjRItLS0mjatGmVnqcyiHAj3JdXywB6vf4OyGSEb9vEkXUrqvR8bj62DPmwHW6NbdEUatn6SwQH10aj0959iVSoGjKZDIsOHfD46Ucabd2C3fDhyNRqii5d4sYHk4nu/ixpv8wTXckFAZgzZw5mZma4u7vTvXt3Zs+ezezZsxk+fPgDGxU/jvHjxxMQEMDPP//M5cuXCQgIICAggMTExCo7J4BarcbBoXZMoyO6ggsPdHrrX+xa+DMA3cf/i1bP9q7S82m1Oo78cYXwv68B4O5rS4/xLTC3NnnAkUJV0GZmkrFyFRlLl1KSmgqAzNwc20GDsB89ChMPDyNXKNR2YvoFoZSYfkGoNgE9+9L++aEA/D3/Jy4fPVil51Mo5HR8oTGhE1qgMlWQcDmTVdOPceOKuFpgDApbWxxfeRmfnX9Tb8YMTH19kfLzyVi6lCuhPbn+5lvknz5t7DIFQRAMRLgRKqTD4OG07NYTJInN/5tJfMTZKj+nT6AzL3zQFjtXc/Kyiln/7SnO7r5WZW1/hPLJTEywHTgA7z/X4zH/Vyyeegp0OnK2b+fqS8OIG/oS2du2I2m1xi5VEIQnnAg3QoXIZDK6jZ9I46AOaEtK+PObaSTHXnnwgY/Jvp4FL3zQlkZtnNHpJPavjGLHggtoisQHqLHIZDIsO3bE89d5eG/4E5tBg5CpVBSEh5Pw1ltcCe1J+tJl6AoKjF2qIAhPKBFuhAqTyxX0fuNdPPz8KS4oYN2MT8hIqtoGbAAmZkpCJzSn4ws+yOQyoo4ns+arE2Qm51f5uYXymfn64vbldHx27cTh1VdQ2NiguX6d5C++ILprN1LnzKEkQ0ySKghC9RLhRngoShMT+r83BWevRuRnZbJ2+hRyM9Kr/Lyl3cUHvB2AubUJ6Yl5rJ5xnJjTqVV+buHBlE5OOP/73/js2Y3L1Cmo6tdHm5FB2vc/EN2tO8kzZqC5ccPYZQo1nLjlLFTWz4AIN8JDMzU3Z9DkT7F1qUdWSjLrZnxCYV7VzYB7J7fGtgz5qB31fGwoLtSy5edzHP5DdBevKeRqNfbDhtFo6xbcvv0G06ZNkfLzSV+8hOhne5A46QOKoqKMXaZQw5SOkpufL67GPulKx+pRKB5v0FjRFVx4ZJnJSayY+h55mRnUb9aCQR9+hsqkYvOyPC6tVsfhdVc4s/NWd/EmtvQYJ7qL1zSSJJF34CA3f/2V/KNHDestu3TBYcJ4zNu0MWJ1Qk1y48YNMjMzcXZ2xtzcXIyI/QTS6XQkJiaiUqnw9PS862fgYT6/RbgRHktKXAwrP/2A4oJ8fNq1p+/bk5E/ZuJ+GFEnktn12yVKirRY2JrS8+UWuDas3JnMhcpRcO4cN+f9Ss6OHXDr1466TRscxo8XM5ILSJJEUlISmZmZxi5FMCK5XI63t3eZObtKiXBTDhFuKt+1C+dY++VUtBoNLbr0oMcrb1TrX13piXls/eUcGUn5yBUynhosZhevyYpiY0lfsJCs9euRbk3KatrYB/tx47Dp0wfZHRP5CU8erVaLpoom6xVqPhMTE+T3+UNHhJtyiHBTNaKOH+avb2cgSTqCBw7hqaGjqvX8xYUl7FpykSun9A2MfYNd6DysKSrT6ruKJDwcTUoKGUuWkLFiJbpcfZstZb16OIwZje3gwcjNzY1coSAINYkIN+UQ4abqnN25jR2/fA9Al9ETaNO7f7WeX5Ikzuy8xqF1V5B0Eg7uFvR82R9bF/EhWZNpc3LIWLGC9CVL0KamAfpRke1Hj8Ju+HAU4v+pIAiIcFMuEW6q1tE/VnFgxRIAer/xLs2e6lztNSRGZbB13nkKsosxMVPQbYwfDVs7VXsdwsPRFRWR9eef3Px1Ppr4eADkFhbYDRuG/ZjRKGvJhH2CIFQNEW7KIcJN1ZIkiT2L53FqywbkCgUD3p+Kd+vAaq8jL7OIbb9GcCNaPx9Vm1BPgvs1RK4QjVZrOqmkhOyt27j588+GbuMyMzNsBw/GYVwYKldXI1coCIIxiHBTDhFuqp6k07H5h2+5dHAvSlNThkz5knqNm1R7HaK7eO0m6XTk7tlD2tyfKTx7ay4zlQrbAf1xGD8ekwYNjFugIAjVSoSbcohwUz20JRrWfz2NuDOnMLOyZuinX+FQ38MotZTpLm5jQs9X/EV38VpEkiTyDx8mbe7P5B87pl8pl2PdqxcOr7yMma+vcQsUBKFaiHBTDhFuqk9xYQGrp31EUvRlrBycGPr511g7GqftS/qNPLb+fKu7uFxGx8E++HeuL7qL1zL5p06T9vNc8vbuM6yz7NYNx1dfQe3vb8TKBEGoaiLclEOEm+qVn53Fyk8mkZ54HXt3D4Z+9hVqK+N83/XdxS9x5VQKAI3budBlhOguXhsVXrxI2s+/kLNtm2FAQIsOHXB45RXMg9qJ0CoIdZAIN+UQ4ab6Zael8PvU98m9mUa9xk144eMvMDFTG6UWSZI4u+s6h9ZGo9NJ2LtZ0OsV0V28tiqKieHmL/PI+usv0GoBUAcE4PjqK1g884wIOYJQh4hwUw4Rbozj5vV4VnwyicLcHDxbtGLgpE9Q3mN47eqSGJ3Jtl8iyM8uRmWmoNvoZjQKcDZaPcLjKb6eQPqC+WSuWYt0a+I9U79mOL78ClY9nhVTOwhCHSDCTTlEuDGeG9GRrJ72MZrCAhoGBtHvnQ9RKJVGqycvq4jtv54nMSoTgIBnPWk/QHQXr800KSmkL1pMxooVSLdmmDZp1AjHV17GundvZEb8eRME4fGIcFMOEW6M69r5s6yb8SklmmKaduxEr9ffQS43XpsXrVbHkT+uEP63vru4W2NbeoxvjoVN9cxuLlSNkowMMn5bSvrSpeiyswFQeXjg8PIEbPv3R2bEq4aCIDwaEW7KIcKN8cWcPs6fM79Ap9Xi3y2UZye8bvS2EdEnU9i15CKaIi3mNib0nNCCej62Rq1JeHza3Fwylv9O+sKFaDMyAFC6uuIwbhy2g19AbmZm5AoFQagoEW7KIcJNzRB5+ACb/vs1kqQjsM8AOo0cZ/SAk5GUx5a5t7uLd3jeh5ZdRXfxukCXn0/m6tXcnL+AkhR9bzmFoyMOY8dg++JQFJYWRq5QEIQHEeGmHCLc1Bzndm9n+9z/AdBh8HBCXnjJyBXpu4vvXnqJ6BP6D0Cfts50GdEUEzPRVqMu0BUVkfXHH9yc9yuahAQAFDY22I0ehf2IEWKSTkGowUS4KYcINzXLqc1/snvxPAA6jxpPYJ8Bxi2IW93Fd1/n0Bp9d3G7ehb0eqUFdq7ir/u6QtJoyNq4iZs//0xxXBxwa5LO4cP1k3Ta2xu3QEEQ7iLCTTlEuKl5jqxdwcFVSwF49uU3aNkt1MgV6d2IzmTbvAjysopRmSroOqoZPoGiu3hdImm15GzbRtrcnym6fBnQT9Jp9+KL2IeFoXIR/96CUFOIcFMOEW5qHkmS2L98Ecc3rAWZjD5vvEvTjp2MXRYA+dnFbJsXYegu3rq7B+0HNkIhuovXKYZJOn/8icKICABkKhU2zw/CYfwETOq7G7lCQRBEuCmHCDc1kyRJ7Jz/I2d2bEEml9P33x/QOLiDscsCQKfVcWR9DKd3xANQr5ENPcY3x9JO9LSpayRJIu/gIdLm/kTBiZP6lUolNn374vDyBEy9vY1boCA8wUS4KYcINzWXpNOx9cfvuLB/N3KFgr7vfIhP22Bjl2Vw5XQKuxZfpLhQi5mFiu5j/WjQwsHYZQlVJP/4cdLm/kzewYP6FTIZ1r164vDKq5g1ETORC0J1E+GmHCLc1Gw6nZYtP8zi0sG9yBVK+r/3EQ0D2hm7LIOs1Hy2zTtPanwOAG16NiC4r7cY1bgOKzh7lrS5P5O7a5dhnZiJXBCqnwg35RDhpubTabVs+t9MLh85gEKlYsC7H+PVOtDYZRmUaLQcXBNNxF59V+L/b+/O46Oq7/2Pv2afSSZ7IAtJSIBAgLCHsIqgCFWkRVzQYqXa9ra9ti7c3p/YVm3rVcSlFxWua2vdEMQKKLghIIKyh0VkJwlhS0LWSWaS2c75/TFhIBLGBZKTDJ/n43EeJ/M935n55DySzDvnfM/3pPSIYcIvcrHHyazG4axx3z4qXniBuo/OuhP56NEk/ubXROTlaVydEOFPwk0IEm46Br/Px/K5czi0ZQNGk5kp9z1I134DtS6rmYNby1jzxj68jX5sUYHTVBl95DRVuHMXFlL5wovULl8evBN5RF4eCb/9DZEjR8qkj0K0Egk3IUi46Tj8Pi/v/X02hds2YzRbmDrrIdL79te6rGZqylx8/PJuKo7Wgw7yrs5k6LVZ6PXyARfuPMeOUfnSy9S++y6q1wuAtX9/En/zG+zjxkrIEeIik3ATgoSbjsXn9fLeU49QtH0rJouVqff/hbTeuVqX1YzP62f92wf5et0JALr0jOWqX8jNNy8V3tJSKv/5T2reXoza2AiAJSeHxF//B1ETJqAzaHdjWCHCiYSbECTcdDw+j4elTzzMkV3bMVltXP/Hv9GlV2+tyzrHgS2lfPbGfrzuwGmqq37Rl/Qcmen2UuGrrKTqX69S/eabKC4XAOasLBJ+/R/ETJqEzmTSuEIhOjYJNyFIuOmYvB43S+f8jZLdOzHbbFw36y+k5fTVuqxzVJc6+fil3VQed545TTUpU66muoT4a2qoeuNNql5/HaW2FgBTWhoJv/wlMVOvQ282a1yhEB2ThJsQJNx0XF53I0sf/xslu3dhsli57r4H290YHACfx8+6RQfY88VJIDDp3/g7+hCdYNO4MtGW/PX1VL/1FlX/ehV/ZSUAxqQkEn5xB7E33ojeJj8PQnwfEm5CkHDTsXndjSx78hGO7NqO0WzhJ//9ZzL7D9K6rBYd2FLK2jf342n0Y7YZGXdrjtyb6hKkNDRQs/gdKv/xD3xlZQAY4uOJv/3nxN1yCwa7XeMKhegYJNyEIOGm4/N5PLz/v7MpLNiCwWTiJ//1J7IGtc95RhwVDXzyj68pK3IA0Gd0KqNvzMZkkUGmlxrF46F2yVIqX3oJ77FjAOhjYoi/9Vbif3YrhthYbQsUop37Pp/f7WIgwPz588nMzMRqtTJs2DA2b9583r4vvfQSl112GXFxccTFxTF+/PiQ/UX4MZrN/Pi//kiPocPxe70se/J/OLR1k9ZltSg60cZ1fxjMkB91BR3sWX+CxbO3UHGsTuvSRBvTm83ETbuJ7h99SOqcxzBnZaHU1lIxfz6HrriS8qeewtd0+koIcWE0DzeLFi1i5syZPPTQQxQUFDBgwAAmTpxIeXl5i/0/++wzbrnlFtasWcOGDRtIT09nwoQJHD9+vI0rF1oyGE1ce88seg4fjd/n4/2/P8qBTV9oXVaLDAY9w6d05yd3DyQyxkx1qYvFj21l15qjXGIHTgWgMxqJ+clP6Lb8fbrM/V8svXqhuFxUvvQyh64cT+mjj+ItLdW6TCE6NM1PSw0bNoyhQ4cyb948ABRFIT09nd///vfMmjXrW5/v9/uJi4tj3rx53Hbbbd/aX05LhRfF7+fD+YF7Uen0eq6+cya9R4/Vuqzzaqj3sPq1fRTvqgAgs18CV8zojc0uV9BcqlRVpX7NZ1Q8/zyNu3YBoDOZiLnuOhJ+9UvM6ekaVyhE+9BhTkt5PB62bdvG+PHjg216vZ7x48ezYcOG7/QaLpcLr9dLfHzL84m43W4cDkezRYQPvcHA1b+bSd/Lr0RVFD6Y9xQ7PvlA67LOy2Y3c81v+3HZtJ4YjHqKv6pk4cObObqvSuvShEZ0Oh1RV4wjc9FC0v/xMhFDh6J6vdS8/TaHf3Q1J+6bhbuwUOsyhehQNA03FRUV+P1+kpKSmrUnJSVR+h0Py953332kpqY2C0hnmz17NjExMcElXf4LCjt6vYGJv7mbgRMngaqy6h//x8Z3F7XbUz46nY7+49K4YVYecckRuGo9vDd3B+sXH8Tn9WtdntCITqfDPmoUXV9/ja5vvE7k6NHg91O7bBmFk67l2L330rhvn9ZlCtEhaD7m5kI89thjLFy4kCVLlmC1Wlvsc//991NbWxtcjh492sZVirag0+u54vbfMPz6mwH4YtHrrH39H+024AAkptm58Y9D6TumCwA7Vx3l7Ue3cqpEBhtf6iLy8sh4+SUyF7+N/corQVWp+/AjiqZcx9Hf/icNO3dqXaIQ7Zqm4SYxMRGDwUBZ09wPp5WVlZGcnBzyuU8++SSPPfYYn3zyCf37n38iN4vFQnR0dLNFhCedTseom25l7G2/AmDbiqV8/PzTKP72ezTEZDYw9qe9mHRnf2zRZqpPOnlnzla2fVSMorTfYCbahq1fP9LnzyNr2TKir7kGdDrq16yheNrNlNzxC1xbtmhdohDtkqbhxmw2M2TIEFatWhVsUxSFVatWMWLEiPM+7/HHH+fhhx/mo48+Ii+vfc5vIrQzZNJPmPjbe9Dp9Hz92ae8/7+P4fN4tC4rpMx+idzyYD7dBnZC8atsXFrIkicLqD3VoHVpoh2w9upJl78/RbcVK4i57jowGHB++SVHfnYbxdNvpX7dunZ9lFKItqb51VKLFi1ixowZvPDCC+Tn5zN37lzefvtt9u3bR1JSErfddhtdunRh9uzZAMyZM4cHH3yQBQsWMGrUqODr2O127N9hpk+5WurScXDLBlbMnYPf5yMjtz8/+cOfMdsitC4rJFVV2b+xlM8XHcDb6MdoMXDZjdn0HpWCTqfTujzRTniOHafy5Zeo/fe7qF4vELgTecIdtxN99dVyk04RljrcDMXz5s3jiSeeoLS0lIEDB/LMM88wbNgwAMaOHUtmZib/+te/AMjMzOTIkSPnvMZDDz3EX/7yl299Lwk3l5YjX+1g2ZOP4G1soHNmd66b9RD2uPZ/p25HRQOrXt3LiYM1AGT2T2TcrTlERMsl4+IMb1kZVf98herFi1Gb7kRuTEkh/rbbiL3xRgz2SI0rFOLi6XDhpi1JuLn0lB46wLtz/kqDo5boTp2ZOuuvJKS1/6vmFEVl56dH2fjeYRSfitVuYszNPekxpLMcxRHN+GtqqF64iKo33sBfEZhDSR8VRdzN04i79WeYkuSeZqLjk3ATgoSbS1NN6Unefewhqk+ewBIZyZQ/PEBan1yty/pOKo7V8+kre6g8Xg9At0GduPyWXnIUR5xDcbupfe89qv75Cp6iokCjyUTM5Mkk3HE7lh49tC1QiAsg4SYECTeXLpejlqVPPMzJA/swGI386M6Z5Iwco3VZ34nfp7Dtw2K2fXgERVGxRBoZM60n2UOT5CiOOIeqKNSvWUPlP/5JQ0FBsN0+diwJv7gDW16e/NyIDkfCTQgSbi5tXo+bD555kkNbAjNgj7n1DvKuva7D/KE/dbSO1a/tpeJo4ChOZv9Exk7vRWSMRePKRHvl2r6dqn/+k7pPV0HTn3tr//4k3HEHUVeNR2eQO9SLjkHCTQgSboSi+PnstZfZ/uH7AAyceC3jZvwKfQf5I+/3K2z/+AhbVhSj+FUsEUYuuymbnsOSO0xIE23PXVRE1b9epXbJEtSmqRFMaWnE3Tqd2OuvxxAVpXGFQoQm4SYECTcCApdcb1uxlLWv/wOAzIFDuPbu/4clouNcXVJ5vJ5Vr+4Nzmic0TeBy2/pSXSiTePKRHvmq6yk+s03qX5zAf7aWgD0kZHEXD+V+FtvxZyRoXGFQrRMwk0IEm7E2Q5sXM+H8/8Xn8dNfGoaU+57kLjkVK3L+s4Uv8L2lSVsXl6E4lMxmvXkT+7GgCvS0Bs69N1VRCtTGhqoff99ql57Dc+hw4FGnQ77FVcQf9ttROQPlSOBol2RcBOChBvxTWWFh1j6xMPUV1VitUcx+d77ycg9/y092qPqUiefvbk/OC9OYrqdcbfm0Lmr/IyL0FRVxfnFl1S99irOz9cF2y29exN/221ET7oGvVmuzBPak3ATgoQb0ZL66iree/IRTh7aj95g4Irbf82Aq67RuqzvRVVV9n55ki//fQi3y4dOB/3GpTHsx90wW41alyc6AHdhIVWvv07t0mWoDYFbfxgSEoi75Rbibp6GMTFR4wrFpUzCTQgSbsT5+DwePnnhGfau/wyAARMmMW7GLzEYO9ZU9i6Hhy/eOciBzYEb0trjLIy5uSdZAzppXJnoKPw1NdS88w5Vb7yJr7QUAJ3JRPSkScRN/ym2fv00rlBciiTchCDhRoSiqiqbly5m/cLXAEjpmcPke2cRFd/x/mMt+bqStW/tx1HRCEDWgERG35gtA47Fd6Z6vdStXEnVq6/RsHNnsN3arx9xP/0p0ddcjd4i0xCItiHhJgQJN+K7KCzYwgfznsTtdBIRE8u199xHep+O99+q1+Nn64oidqw8iqKoGEx6Bk/syuAJGRjNHePSd9E+NOzcSdWbb1L34UfBm3UaYmOJuX4qcbfcgjktTeMKRbiTcBOChBvxXdWUnuS9px7hVEkxOr2eMdNvZ8ikKR3yCpLKE/WsW3SQ4/urAYhKsDL6xmyyBiR2yO9HaMdXWUnNO/+metFCfCdOBhp1OuxjxhA3/adEjh6NTi9X6omLT8JNCBJuxPfhdTey8qX57F23BoCew0cz8Td3YbZFaFzZ96eqKocLTvHFOwepr3YDkNEnntE3ZROX3HHm9xHtg+r3U792LdVvLsD5xRfBdlNGBnE330zs1OswxMZqV6AIOxJuQpBwI74vVVXZ8ckKPnv1JRS/n7iUVCbdfR9JWd21Lu0H8br9bPuwmO2flqD4VPQGHQPHpzPk6ky5qkr8IO6iImoWLqTm3SUodYFJJXUWC9E/mkjsjTdiGzJEjhCKCybhJgQJN+KHOr5/Lyuefpy6ylMYjEbG3HoHg340ucP+0a4pc7Hu7YOUfF0JgC3azLDJWfQemSITAIofRHG5qF2xguoFb+HeuzfYbu7WjdgbbiBmyk8wxsdrWKHoyCTchCDhRlyIhvo6Pn7uaQ5v3QhAtyH5/Oi392CL6pg/S6qqUvxVJesXH8RxKjCvSXxqJCOn9iCjb3yHDW5CW6qq0rhrF9WLF+P44ENUlyuwwWQiavyVxN14IxHDh8vYHPG9SLgJQcKNuFCqqrL9o+V8/sY/8Pt82BMSmfT7P5DWO1fr0n4wv09h99rjbPmgCLfTB0BaThyjbuhBYprcUFH8cP76ehwrPqBm8WIad+8OtpvS04m9/npipl6HqXNnDSsUHYWEmxAk3IiLpazoMCuenkP1yRPodHqG/ngqI26cjtHUsSb9O1uj08u2D4vZ9dkxFJ8KOsgZkcLwH3cjMlbmMxEXpnHvXmoWL6b2vfdR6usDjQYD9rFjib1+KvbLLkPXgX9/ROuScBOChBtxMXkaXKx+5QW+XrsKgE4ZmVz9u/+iU9csjSu7MLWnGti47DCHtpYDYDTr6X9FOoOuysAaKR8+4sIoDQ04PvqYmsWLaSgoCLYb4uOJmXwtMVOmYO3dW8MKRXsk4SYECTeiNRzc/CUrX5pPg6MWvcHIyJumM/THU9HrO/ZEeaWFtXzxziFKC2sBMNuMDLoqg/5XpMmVVeKicB86RM07/6Z2+XL8FRXBdktODjFTfkLM5MkYExI0rFC0FxJuQpBwI1qLq7aGT16cFxxsnNqzNz/6z3uIS+micWUXRlVVindVsHFZIVUnnADYokwMntiV3Mu7YDR17AAn2gfV56N+/XpqlyylfvXq4CzIGAzYx4whZsoU7OPGyh3KL2ESbkKQcCNak6qq7Pl8NatfeQFPgwujyczwG24h79rrMBg79pEOVVE5uK2Mze8VUdt0ZVVkrIW8azLpPTIFg1GufBEXh7+mhtoPPqB26TIad+0KthtiYoi6+kfEXHsttsGD5WqrS4yEmxAk3Ii24Kgo55MXnuXIru0AdOqaxYRf30Vy92yNK7twfr/C/g2lbFlRFJzp2B5nYfDErvQelSJHcsRF5T58mNqly6hdtgxfeXmw3ZiSQvQ1VxMzaRKW3r1l2oJLgISbECTciLZy+ijOZ6+9TGN9HTqdnsHX/JhRN92KyWrVurwL5vP6+frzExR8cgRXrQeAiBgzg67KoO+YLpjkxpziIlL9fpwbN+JYvoK6lSvPXG1FYJLA6EnXEDNpEubMTO2KFK1Kwk0IEm5EW3PV1rDm1ZfY98VaAKI7JTHu5/9B9yH5YfHfps/rZ+8XJyn4+EjwSI4tysTA8RnkXt5FBh6Li05xu6lfuxbHig+oX7MG1eMJbrPm5hI9aRLR11yNKSlJwyrFxSbhJgQJN0IrRdu3svLl+dRVnAIga+AQxs74D+JTO/aA49P8PoX9G0vZ9lExjopGACwRRvqO6UL/cWlExsg8OeLi89fXU/fppziWr8C5YQP4/cFttoEDiZowgagJEzCnhcfv2aVMwk0IEm6EljyNDWx6dxFbly9F8fvQG4wMuXYKw6dOw2y1aV3eRaH4FQ5sKWPbh0eoKQtMu6836uiZn8yg8RnEp8odyEXr8FVW4vj4YxzLVzSbPwfA2rcvURMnEj3hKjl11UFJuAlBwo1oD6pOHOezV1+kaMc2AOzxCYyZfjs5I8eEzRUgqqJStKuCHStLOHm4NtjeNTeBQVdlkNozNixOy4n2yVtWTt2nK6n7+BNcW7eCogS3WXr1ImrCVURPnIilRw8NqxTfh4SbECTciPZCVVUKCzaz5tWXqC0rBSCpWw/GTL+djNwBGld3cZUW1rJ9ZQmFO05B01+cxHQ7/cam0XNoEkYZfCxaka+ykrpPV1H3ySc4N25sdurKnJWF/YpxRF1xBbaBA9EZ5GexvZJwE4KEG9He+Dwetq1YyuZli/E0BOaPyRo4hMum306njExti7vIaspc7Fx1lL0bTuL3Bv6TtkQa6TMyldzLuxCdGB6n5kT75a+poW71Guo+/hjnl1+emSwQMMTGYr/8cuzjxhE5ejQGu5xCbU8k3IQg4Ua0Vy5HLRv/vZCdKz9A8ftBp6PvmCsZedNPiU4Mr7smN9Z72fPlCXavPU5dZWDwMTrI7JdI/7FppOXEodPLKSvRuvx1dTjXr6du9RrqP/8cpfbM6VOdyUREfn7gqM64cZhSUzWsVICEm5Ak3Ij2rrr0BOsXvs6BDesA0BuM9LviKvKn3Bh2IUdRVI7sruSrz45xdE9VsD060UrvUan0HpEidyMXbUL1+XAVFFC/5jPqV6/Gc+RIs+2WXr2wXzaayNGXETF4EDq5DUSbk3ATgoQb0VGcPLSf9W+9RsnunUB4hxyA6lInu9ceZ9+Gk3gaA2MidDro2i+RPqNS6JqbgN4QHoOtRfvnLiyifs1q6tasoaFge7MBybqICCKHDSPystHYR4/GnJGhYaWXDgk3IUi4ER3Nsb272fDOW81CTu648QydfD2xySkaV3fxeT1+DheUs2f9CU4eOnOaICLGTM7wFHoNTyY+RcZCiLbjq67Guf4LnOvXU//FF83uXg5g6pqBffRlRI4eRWR+PvpI+flsDRJuQpBwIzqqY3t2s+HfCyjZHbiRoE6nJzt/BHk/nkpKj14aV9c6qkud7PniJPs3nqSh7szAz04ZUfTMTyJ7aJJMDijalKoouPfto379FzjXrcO1fTv4fGc6mEzYBvQncthwIoblYxs4UO5kfpFIuAlBwo3o6I7t2c3mZYuDc+QApPXJZejk68kaOCRs5sk5m9+nULSzgv2bSinZXYmiBP5s6XSQlhNHz/xkug3qJLd6EG3OX1+Pa9Mm6tevx/n5OrzHjzfbrrNaiRg8iIhhw4kclo81NxedUX5OfwgJNyFIuBHhoqKkmK3Ll7B3/VoUf+A/x/jUNAZMmETfy6/EEhGhcYWto6Hew6Gt5RzYXEZp4ZnTVkaTnsz+iXQf3JmMvvESdESbU1UVb0kJzk2bcG3chHPz5nNOYekjI4nIyyNi2DAihuZhzclBZzJpVHHHIuEmBAk3ItzUVVZQ8OF77Pr0w+A8OSaLlT5jxjFwwiQSw2yunLPVnnJxYHMZBzaXBW/1AGAw6cnoE0/3wZ3J7J+IxSZBR7Q9VVXxHD6Mc+MmXJsCYefsy80BdDYbtgEDiBg8GNuQwdgGDJT5dc5Dwk0IEm5EuHK7XOxZt5odH6+g6vjRYHta71z6XTmR7PwRmCxWDStsPaqqcqqkjkNbyzm8vTx4404AvUFHWk483Qd3ImtAIja7jH8Q2jg9Xse5cROuzZtxbd9+TthBr8eak4NtyBAihgzGNngwps7hd3XkDyHhJgQJNyLcqarK0a+/Yscnyzm0ZSNq0yWsZlsEvUZeRu7Yq0jJ7hW293VSVZWKY/UUbj/F4YJyqkvPHNFBB8lZ0XTNTaBrbiKJ6faw3Q+i/VMVBc/hw7i2FeAq2EbDtoJzxuwAmFJTsfbvj61/f2wD+mPt0we97dKbzVvCTQgSbsSlpK6ygt1rVvL12k+pLS8LtsenptF37HhyRo4hulN4/1dYddJJ4fZyDm8/RcXR+mbbImPMwaCT1jtOxukIzXlLS2koKGgKPAW49+2Db35MGwxYevYMhJ3+/bD174+5W7ewvy+WhJsQJNyIS5GqKBzds5uvP1vJgU1f4vO4g9tSsnuRM3IMPYePxh6foGGVra++upEjuysp/qqSY/uq8HnOTMymN+hI6RFDWk48aTlxdM6IkkkDheb89fU07v6ahl27aNi1k8adu/CdOnVOP31kJNbevbH27YOld2+svftg6d4trK7MknATgoQbcalzu1zs37COfes/4+je3Wf+K9TpSMvpS6+RY8jOH0FkbJy2hbYyn9fPiYM1HPmqkuLdlThONTTbbrYaSO0ZR1pOYIlPiZRTWEJzqqriKyujYWdT2Nn1FQ1ff43qcp3TV2exYOnZMxB6+vTG2qcPlp490Vs75tg7CTchSLgR4oz66ioObPyC/V9+zokDe89s0OlI6dGT7kOG0W1IPonpXcP6g11VVWrKXBzbV82xfdUcP1CN2+Vr1ici2kyXnrGk9IglpUcM8al29HJzT9EOqD4f7sOFNO7ZE1j27sG9dx+K03luZ4MBS7csLDm9sfTMxpKdjTU7G2Nqarv/HZdwE4KEGyFa5qgo58CG9ezfsI7SwwebbYvulET3Ifl0G5JPWu9cjGE+L4eiqFQcrWsKO1WcOFSL36s062O2GkjuHkNK90DYScqMxmgO7zEPouNQFQVvSQmNe/fSuGdvU+jZi7+qqsX++shILD16BAPP6cWQkNBuQo+EmxAk3Ajx7eqrKiks2MLhbZso+WonPq8nuM1otpDWuy9d+w0ko99AOmVkhuWsyGfzexVKC2s5ebiGE4dqKT1ci9ftb9ZHb9CRmGanc9doOmdG0zkzirjkSDm6I9oNVVXxlZfT+PUe3AcO4D54MLAUFYHX2+JzDHFxWHr0wJyV1bRkYsnKwtSlS5uP55FwE4KEGyG+H6+7kSNf7aRw2yYKC7bgrKlutt0WHUNG7gC69htIWp9cYpNS2s1/eq1F8StUHndy8nANJw/VcuJQDa5azzn9TBYDnTKiAmGnaxRJmdFEJVjDfv+IjkX1evEcOXIm7Bw8iPvAQTwlJedeqXWayYQ5IyMYdsyZZ8KPMa51xutJuAlBwo0QP5yqqlQePcKRr3Zy5KvtHNuzG6+7sVmfiJhYuuT0oUuvvnTJ6UPnzG7ow/wSVVVVqatspKzYQXmxg/IjdZSX1OH7xtEdAEuEkYQudhLS7CQ2reNTIzHJKS3RziiNjXgKC3EfOoS7qAhPUTGeoiI8xcWobvd5n2eIicE2NI/0efMuaj0SbkKQcCPExeP3eTl5cD9HvtpJye6dlB0+gN/XfCCuyWIlJbsnSd17ktytB0ndsonu1Dnsj14oikr1SSflRxyUF9dRfsRBxbF6FP+5f3J1OojpHEFCFzuJaZHEp9iJS4kgupMNg1yOLtoZVVHwnjgZCDpFRXiKi4Lhx1daCkDEiOF0feWVi/q+Em5CkHAjROvxeTyUFR7i+P49HN/3Ncf378HdwhUbtqhokrpnB8NO56zuRCUkhn3g8fsUqkudVB6rp+JYPZXHA+uGupbHO+j1OqI72YhLjiAuOZK45Ahim76W+2WJ9khxufAcPQp+P9Y+fS7qa0u4CUHCjRBtR1UUKo+VcOLgfsoKD1JWeIhTR4qDdzE/myUyksT0riSmZwbWGYGvrXa7BpW3LWetOxh0Ko/XU33SRXWZq8XTWqdFRJuJTYogOtFKdKKNmE42ohMDiy3KFPZBUVx6JNyEIOFGCG35PB4qSoopLTxEWeFBSg8fpPJYSfAeWN9kj08gIS2DuJQuxKWkNq27EN2pE3p9+I5TUVUVZ42b6lIX1aXOprWLmlInzhYGL5/NZDE0BR0r0Z1sRCfYsMdZiIq3Yo+zYLVL+BEdj4SbECTcCNH++Lxeqk8co6KkmFNHj1BRUkzF0SPUVZw7zfxpBqORmKSUYOiJTUohulPnwJLYKWzvgA7gbvBRU+qitsKF41QDtRWNOE414KhooL7GDd/yV91g0mOPtWCPt2CPDQQee1PwscdZiYg2Y7Ob0Mll7KIdkXATgoQbIToOt8tJRckRqk4co7r0BNUnjlN98jg1ZSfxn2dejtMiYmKJTuxEdKekYOiJSuiEPS4ee1w8ETGxYXkVl8/rp66yEUdFI46KBmpPNVBX1Uh9VSP11W5cdZ5vDT8AOr0OW5SJiGgzEdEWImLMTV9/Y4mxYLYa5EiQaHUSbkKQcCNEx6cofuoqKqg+eTwYemrLS3GcKsdRUY6noeFbX0On0xMRG4s9Lp7IuHjssU3r+HgiY+OwRcUQER2DLToGs80WNh/efp+Cs8ZNfXUjdVWBdX21u2lpxFnjpqHe+50C0Gl6gw5rpAmr3dTyOtKEzd78sSXCKEeGxPci4SYECTdChDdVVXE7ndSeKsNRUY6jPBB4HKfKqK+qpL66CmdN9XnH+LTEYDRii4oOLE2BxxYVHQg/UdFY7HasEZFYIiOxRNix2u1YIiIxms2t+J22HsWv0FDvxVXrweXw4HK4A+vg46al1o2n8fyDnkPSgdliwBxhxGIzYradWZu/8dhiM57pZzVithkwWQwYLQa5VP4S0uHCzfz583niiScoLS1lwIABPPvss+Tn55+3/+LFi3nggQcoLi4mOzubOXPmcM0113yn95JwI4RQFD8NDseZsFNddda6EmdNNQ11DlyOWnwhJiv7NgaTCUtEJJbIpvDTFHrMNhtmqw2T1YbZaj1rbW2hLdDXYGqfg4B9Hj8N9V4anU1LfdPibGHt9NJQ78X7QwNRC/RGHSaLAZM5EHhMFgMma/PHRstZ2yxGTBY9RrMBg1GP0azHaDJgMOkxmgKPDUZDYG3SYzTq5QhTO/F9Pr81nyhh0aJFzJw5k+eff55hw4Yxd+5cJk6cyP79++ncufM5/b/88ktuueUWZs+ezbXXXsuCBQuYMmUKBQUF5ObmavAdCCE6Gr3eQGRsHJGxcSR9S1+vu5GGOgcNDgcNjlpcp7+uq8XlqKXB4cDtcuJ2OnG76nE7nTS6nKCq+L1eXLU1uGprLrhmnV6PyWLFaDZjNFsCa5O56bE5+NjwjcfB/iYTRrMFg8kUWAxG9EYjhqZFbzRhMBia2kxNbcZma4PRiN5gbBayjGYDUfEGouK/+wBuv0/B7fLhafCdWTe0vA72aTzztbfRj6IE/i9XfCpunw+389zpBS4Wg1EfDD9nQtCZcGQw6TEY9OiNuvOvjXr0htNrPYbTbcE+egwGXfO1UYdOr0Ov16E3nP468Dp6vQ6dQYdOx5nHel27DMBa0PzIzbBhwxg6dCjzmqZpVhSF9PR0fv/73zNr1qxz+k+bNg2n08ny5cuDbcOHD2fgwIE8//zz3/p+cuRGCNHaVEXB09gQCDrO+mD4aXQGwo+n0YW3sRFPYyPexoamrwNrb2MDnsaGpm2N+Dw//MhRa9EbjOiNhqbAY0Kv16MzGJo+eA3odE1r/Vlr/ZnHZ9r06PQG9Iam9em2pj7B7U1tOp0OnV6PTqdDVUH1g6IEZoNW/KD4VVSF4NeKH/x+FcWv4vcHgpDfFzjtpvibnudreo2mfoHncNaYI903lqa2ZiGipa/Pt/1Mm+6b28/7mt/9fXT6QBAOhiLdmRAU2Oc69AY9NG0L7Hua9jVn7ePAoPLTgen083VN37pOp2/axjfWgedEJ0Yy7Mf9WvgefrgOc+TG4/Gwbds27r///mCbXq9n/PjxbNiwocXnbNiwgZkzZzZrmzhxIkuXLm2xv9vtxn3WYWWHw3HhhQshRAg6vT5wOioikuhO5x6B/j4UxY+30d0UegJhx+/14vO48Xk8gcXbtPa48Z21zR9sb95P8fnwNy2Kv2n9jceBNu85t9MAUPyBfhdyyk6EN0tkGsN+/O0HHFqLpuGmoqICv99PUlLzA8NJSUns27evxeeUlpa22L+06X4W3zR79mz++te/XpyChRCijen1BiwREVgiIjR5f1VVURUFf1PQUfx+/D7vNwKSH1VRzqyVcx+f26epn7/p6/P1VQNrxR8Yp6MqSqAmVQEVVFVBVQKPA7WqgPqNtsAaVQ0+N7D99Ouc/vrsvkrg6NA3+yoqKmpwQHqgb9O+ChR0esedOfjT9Npn79NmfYLbzuqnNr1e8Pm00OfM16ra/PWa9zn9vOALN6v13PYz79Vse/PmM9/zWd/naRGx2vy8nqb5mJvWdv/99zc70uNwOEhPT9ewIiGE6Dh0Ol3wNJHJonU1Qnw3moabxMREDAYDZWVlzdrLyspITk5u8TnJycnfq7/FYsFikd9IIYQQ4lKh6QQBZrOZIUOGsGrVqmCboiisWrWKESNGtPicESNGNOsPsHLlyvP2F0IIIcSlRfPTUjNnzmTGjBnk5eWRn5/P3LlzcTqd3H777QDcdtttdOnShdmzZwNw9913c/nll/PUU08xadIkFi5cyNatW3nxxRe1/DaEEEII0U5oHm6mTZvGqVOnePDBByktLWXgwIF89NFHwUHDJSUl6PVnDjCNHDmSBQsW8Oc//5k//vGPZGdns3TpUpnjRgghhBBAO5jnpq3JPDdCCCFEx/N9Pr/lphxCCCGECCsSboQQQggRViTcCCGEECKsSLgRQgghRFiRcCOEEEKIsCLhRgghhBBhRcKNEEIIIcKKhBshhBBChBUJN0IIIYQIK5rffqGtnZ6Q2eFwaFyJEEIIIb6r05/b3+XGCpdcuKmrqwMgPT1d40qEEEII8X3V1dURExMTss8ld28pRVE4ceIEUVFR6HS6i/raDoeD9PR0jh49KvetakWyn9uG7Oe2Ifu57ci+bhuttZ9VVaWuro7U1NRmN9RuySV35Eav15OWltaq7xEdHS2/OG1A9nPbkP3cNmQ/tx3Z122jNfbztx2xOU0GFAshhBAirEi4EUIIIURYkXBzEVksFh566CEsFovWpYQ12c9tQ/Zz25D93HZkX7eN9rCfL7kBxUIIIYQIb3LkRgghhBBhRcKNEEIIIcKKhBshhBBChBUJN0IIIYQIKxJuLpL58+eTmZmJ1Wpl2LBhbN68WeuSws7s2bMZOnQoUVFRdO7cmSlTprB//36tywprjz32GDqdjnvuuUfrUsLS8ePHufXWW0lISMBms9GvXz+2bt2qdVlhxe/388ADD5CVlYXNZqN79+48/PDD3+n+ROL8Pv/8cyZPnkxqaio6nY6lS5c2266qKg8++CApKSnYbDbGjx/PwYMH26w+CTcXwaJFi5g5cyYPPfQQBQUFDBgwgIkTJ1JeXq51aWFl7dq13HnnnWzcuJGVK1fi9XqZMGECTqdT69LC0pYtW3jhhRfo37+/1qWEperqakaNGoXJZOLDDz9kz549PPXUU8TFxWldWliZM2cOzz33HPPmzWPv3r3MmTOHxx9/nGeffVbr0jo0p9PJgAEDmD9/fovbH3/8cZ555hmef/55Nm3aRGRkJBMnTqSxsbFtClTFBcvPz1fvvPPO4GO/36+mpqaqs2fP1rCq8FdeXq4C6tq1a7UuJezU1dWp2dnZ6sqVK9XLL79cvfvuu7UuKezcd9996ujRo7UuI+xNmjRJveOOO5q1TZ06VZ0+fbpGFYUfQF2yZEnwsaIoanJysvrEE08E22pqalSLxaK+9dZbbVKTHLm5QB6Ph23btjF+/Phgm16vZ/z48WzYsEHDysJfbW0tAPHx8RpXEn7uvPNOJk2a1OznWlxc7733Hnl5edx444107tyZQYMG8dJLL2ldVtgZOXIkq1at4sCBAwDs3LmT9evXc/XVV2tcWfgqKiqitLS02d+PmJgYhg0b1mafi5fcjTMvtoqKCvx+P0lJSc3ak5KS2Ldvn0ZVhT9FUbjnnnsYNWoUubm5WpcTVhYuXEhBQQFbtmzRupSwVlhYyHPPPcfMmTP54x//yJYtW7jrrrswm83MmDFD6/LCxqxZs3A4HOTk5GAwGPD7/TzyyCNMnz5d69LCVmlpKUCLn4unt7U2CTeiQ7rzzjvZvXs369ev17qUsHL06FHuvvtuVq5cidVq1bqcsKYoCnl5eTz66KMADBo0iN27d/P8889LuLmI3n77bd58800WLFhA37592bFjB/fccw+pqamyn8OYnJa6QImJiRgMBsrKypq1l5WVkZycrFFV4e13v/sdy5cvZ82aNaSlpWldTljZtm0b5eXlDB48GKPRiNFoZO3atTzzzDMYjUb8fr/WJYaNlJQU+vTp06ytd+/elJSUaFRRePrv//5vZs2axc0330y/fv342c9+xr333svs2bO1Li1snf7s0/JzUcLNBTKbzQwZMoRVq1YF2xRFYdWqVYwYMULDysKPqqr87ne/Y8mSJaxevZqsrCytSwo7V155JV999RU7duwILnl5eUyfPp0dO3ZgMBi0LjFsjBo16pypDA4cOEDXrl01qig8uVwu9PrmH3UGgwFFUTSqKPxlZWWRnJzc7HPR4XCwadOmNvtclNNSF8HMmTOZMWMGeXl55OfnM3fuXJxOJ7fffrvWpYWVO++8kwULFrBs2TKioqKC525jYmKw2WwaVxceoqKizhnDFBkZSUJCgoxtusjuvfdeRo4cyaOPPspNN93E5s2befHFF3nxxRe1Li2sTJ48mUceeYSMjAz69u3L9u3b+fvf/84dd9yhdWkdWn19PYcOHQo+LioqYseOHcTHx5ORkcE999zD//zP/5CdnU1WVhYPPPAAqampTJkypW0KbJNrsi4Bzz77rJqRkaGazWY1Pz9f3bhxo9YlhR2gxeWVV17RurSwJpeCt573339fzc3NVS0Wi5qTk6O++OKLWpcUdhwOh3r33XerGRkZqtVqVbt166b+6U9/Ut1ut9aldWhr1qxp8e/xjBkzVFUNXA7+wAMPqElJSarFYlGvvPJKdf/+/W1Wn05VZZpGIYQQQoQPGXMjhBBCiLAi4UYIIYQQYUXCjRBCCCHCioQbIYQQQoQVCTdCCCGECCsSboQQQggRViTcCCGEECKsSLgRQgghRFiRcCOE6FD2799PcnIydXV1rfYeN998M0899VSrvb4QonXJDMVCCM2NHTuWgQMHMnfu3G/tO3XqVIYMGcKf/vSnVqtn9+7djBkzhqKiImJiYlrtfYQQrUOO3AghOoySkhKWL1/Oz3/+81Z9n9zcXLp3784bb7zRqu8jhGgdEm6EEJr6+c9/ztq1a3n66afR6XTodDqKi4tb7Pv2228zYMAAunTpEmz717/+RWxsLMuXL6dXr15ERERwww034HK5ePXVV8nMzCQuLo677roLv98ffN7//d//kZ2djdVqJSkpiRtuuKHZe02ePJmFCxe2yvcshGhdRq0LEEJc2p5++mkOHDhAbm4uf/vb3wDo1KlTi33XrVtHXl7eOe0ul4tnnnmGhQsXUldXx9SpU7nuuuuIjY3lgw8+oLCwkOuvv55Ro0Yxbdo0tm7dyl133cXrr7/OyJEjqaqqYt26dc1eMz8/n0ceeQS3243FYrn437gQotVIuBFCaComJgaz2UxERATJyckh+x45cqTFcOP1ennuuefo3r07ADfccAOvv/46ZWVl2O12+vTpw7hx41izZg3Tpk2jpKSEyMhIrr32WqKioujatSuDBg1q9pqpqal4PB5KS0vp2rXrxfuGhRCtTk5LCSE6jIaGBqxW6zntERERwWADkJSURGZmJna7vVlbeXk5AFdddRVdu3alW7du/OxnP+PNN9/E5XI1e02bzQZwTrsQov2TcCOE6DASExOprq4+p91kMjV7rNPpWmxTFAWAqKgoCgoKeOutt0hJSeHBBx9kwIAB1NTUBPtXVVUB5z9FJoRovyTcCCE0Zzabmw32PZ9BgwaxZ8+ei/KeRqOR8ePH8/jjj7Nr1y6Ki4tZvXp1cPvu3btJS0sjMTHxoryfEKLtyJgbIYTmMjMz2bRpE8XFxdjtduLj49Hrz/3fa+LEifzyl7/E7/djMBh+8PstX76cwsJCxowZQ1xcHB988AGKotCrV69gn3Xr1jFhwoQf/B5CCO3IkRshhOb+8Ic/YDAY6NOnD506daKkpKTFfldffTVGo5FPP/30gt4vNjaWd999lyuuuILevXvz/PPP89Zbb9G3b18AGhsbWbp0Kb/61a8u6H2EENqQGYqFEB3K/Pnzee+99/j4449b7T2ee+45lixZwieffNJq7yGEaD1yWkoI0aH8+te/pqamhrq6OqKiolrlPUwmE88++2yrvLYQovXJkRshhBBChBUZcyOEEEKIsCLhRgghhBBhRcKNEEIIIcKKhBshhBBChBUJN0IIIYQIKxJuhBBCCBFWJNwIIYQQIqxIuBFCCCFEWJFwI4QQQoiw8v8By0rZgHiy9ccAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "for my_α in [0 * mM, 0.2 * mM, 0.4 * mM, 0.6 * mM, 0.8 * mM, 1 * mM]:\n", " plot_α(my_α)\n", "\n", "plt.gca().legend()\n", "plt.xlabel(\"t (ms)\")\n", "plt.ylabel(\"[u] (mM)\")\n", "plt.show()" ] } ], "metadata": { "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.10" } }, "nbformat": 4, "nbformat_minor": 1 }